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Zur Kinematik des Rollgleitens 
Von Gerhard Grüß { in Freiberg 


Zur Klärung der kinematischen Verhältnisse beim Rollgleiten werden die beiden Polkurvenscharen, die 
raumfesten Spurkurven und die für den sich bewegenden Körper, körperfesten Rollkurven für alle Bauas NgS- 
zustände vom reinen Gleiten bis zum reinen Rollen bestimmt. 


To investigale the kinematic relations of roll-gliding, the sets of the two polecurves — the fived polecurves 
as well as the movable polecurves of Ihe moving solid — are evaluated jr all states of motion, from absolute 
gliding to absolute rolling. ; 


Pour Eclairer les conditions cineligues du glisser roulant les deux faisceaux de courbes polaires, les traces 
fixes aussi que les courbes roulantes du corps mouvant, sont determ Ws pour tout Etat de mouwvement depuis 
le rouler pur jusqu’au glisser pur. 


C UHeJAIbIO HCCJIENOBAHUA KUHEMATUKH KATAHUA CO CKONbIKEHHEM ONPENENAMTCA MIA BCEX 
CAIYyYaeB ABHIKEHHA OT YUCTOTO CKOJIBIKEHNA AO YUCTOTO KATAHAUA 06a ceMeücTBa Ha HIPaBJIAIOIIHX 
KPHBbIXJIHHHM CJIene Ha INIOCKBCTU KATAHNA WU CBASAHHBIe C ABH3KYIIMHMCA TEeJIOM JIHHHH 
KAaTaHHag. 


1. Problemstellung 


Viele Bewegungen, denen in der Technik eine gewisse Bedeutung zukommt, sind als ein 
aus Rollen und Gleiten gemischter Vorgang zu verstehen, der hier als ‚‚Rollgleiten‘ bezeichnet 
? werden soll; man denke etwa an die Anfangsbewegung der Triebräder einer Lokomotive auf 
E glatter Schiene oder an die gegenseitige Bewegung der Zähne zweier Zahnräder. Das ebene 
. Rollgleiten eines starren Körpers an einem zweiten, der als raumfest angesehen werden soll, 
R kommt dadurch zustande, daß eine mit dem beweglichen starren Körper festverbundene Kurve € 
E an einer mit dem raumfesten Körper verbundenen Schiene C teils abrollt und teils gleitet. Im 
folgenden soll versucht werden, zur Kenntnis der kinematischen Verhältnisse dieses Bewegungs- 
 typus dadurch beizutragen, daß die beiden Polkurven, die raumfeste Spurkurve und die körper- 
feste Rollkurve, des im Rollgleiten befindlichen Körpers bestimmt werden. Dazu ist zunächst 
F die Fragestellung mathematisch zu präzisieren. 


Die glatte, reelle Kurve Ü sei in der festen x,y-Ebene mit Hilte ihrer Bogenlänge s als 
; Parameter durch die Gleichungen 


226), y=yßa) e +) = NEE U NN 


bestimmt. Entsprechend sei in einer zweiten Ebene, die gegenüber der x,y-Ebene beweglich 
und durch das kartesische Koordinatensystem z, ı) orientiert zu denken ist, die glatte Kurve € 
mit der Bogenlänge 3 gegeben: 


R ger (9) i 
wre ri), y=y($), ee lgs Era nn 5. (2), 
he das Koordinatensystem so gewählt sei, daß der Punkt na —=(0 von & Koordinatenanfangs- 
punkt und die g-Achse Tangente im Punkt $ = 0 sein möge, eine Voraussetzung, durch die ja 
die Bueoeinbell der Fragestellung nicht eingeschränkt wird (Bild 1): 


d d 
x(0) = 90) 0, EEE ao BEL E  ea 


Br Die Kurve € möge sich nun längs Ü een Voraussetzungen gemäß bewegen: 


Kar 
ie eL; bei Beginn der Bewegung, also zur Zeit &—=0, mögen sich die Anfangspunkte s = 0 
n Bu 80 der beiden Kurven decken; 


1 
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2. in jedem Zeitpunkt ti der Bewegung mögen die beiden Kurven dieselbe Gerade zur 


Tangente haben; 
3. die längs C und & gemessenen Bogenlängen, die von dem Anfangspunkt s = 0, 83—= 0 
£ und dem zur Zeit gemeinsamen Berührungs- 
punkt s(t) und 3(t) begrenzt werden, mögen 
ein zeitlich konstantes Verhältnis bilden: 


3 — As mit konstanter Rollgleitzahl A. 


Durch diese Bedingung sei der Vorgang 
des Rollgleitens definiert. Als Grenzfälle sind 
im Rollgleiten enthalten: 

erstens für A=0 der Fall, daß die 
Kurve E an der Schiene Ü so entlang gleitet, 
daß der Anfangspunkt = 9 —= 0 von Estets 
Berührungspunkt bleibt; das ist der Fall des 
reinen Gleitens von € auf O; 

zweitens der durch A=1 gekennzeich- 
nete Fall reinen Rollens. 


Das Rollgleiten im engeren Sinn ist also 
durch einen von 0 und 1 verschiedenen Wert 
der Rollgleitzahl A gekennzeichnet, die auch 
negativ oder größer als 1 seinkann. Natürlich 
wird bei wirklichen Bewegungen die Voraus- 
setzung, daß A zeitlich konstant sei, nicht oder 

X wenigstens nicht streng erfüllt sei; doch soll 
/ hier mit konstantem 4 gerechnet werden, um die 
Bild 1 Verhältnisse nicht allzusehr zu komplizieren. 


X,(5) X(S) 


2. Translationspunkt und Winkelgeschwindigkeit 
Zur Bestimmung der raumfesten Polkurve oder Spurkurve dienen die Gleichungen 


ee meyt+n We RE RE TE 
darin bedeuten &,, %, die Koordinaten des sogenannten Translationspunktes, x, und %, seine 
Geschwindigkeitskomponenten und & die Winkelgeschwindigkeit des starren Körpers, der hier 
im Falle des Rollgleitens mit der Kurve & festverbunden zu denken ist oder auch durch sie 
repräsentiert wird; das körperfeste Koordinatensystem ist also das System r, ), dessen r-Achse 
Anfangstangente von € ist; sein Anfangspunkt x = 4 = 0 werde als Translationspunkt gewählt. 


Schließlich ist, wenn man den Winkel zwischen den beiden Abszissenachsen y„=-Oundy=0 


mit @ bezeichnet, die Winkelgeschwindigkeit ® = . 


Die raumfeste Polkurve der rollgleitenden Kurve & läßt sich also berechnen, wenn nur 
die drei Funktionen &,(t), Y.(t) und &(t) den Bedingungen des Rollgleitens gemäß bestimmt sind. 


Ebenso kann dann die körperfeste Polkurve oder Rollkurve nach den bekannten Formeln 
IRESSR 
Em = ar (2 sin 9 — %0 605 9) 


IyR® eher. 
1m =; (io p + jasin 9) 
leicht berechnet werden. 


Zur Bestimmung von &, %, und & stehen nun folgende Gleichungen zur Verfügung: 
erstens fällt der Punkt ı@). y(8) der Kurve € gemäß der a a. 


BE Asa 1 (6) 
mit dem Punkt x(s), %(s) der Kurve C zusammen, so daß die Formeln der Koordinateriie Zeiss 


formation gelten: 
s()= © + x(As) cos P—Y(As)sinp,. | | 
y(s)= Yy+ t(As)sin@ + Y(As) cos p 


(6) ° 
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und zweitens haben die Kurven O’und Ein ihrem gemeinsamen Punkt P— ®dieselbe Tangenten- 
richtung, die mit der x-Achse den Winkel d und mit der g-Achse den Winkel 9 — p einschließt: 
dr 


de 
ds = 059 cos (d —p) cos p — sin (d — p) sin 9= 7, 


. cos) 
dy _ sin d — sin(d — 9) cos p + cos(d 2 
Fake i Y p ( pin g—zz not ggsinn 


sin, 


Im folgenden soll nun die Ableitung nach der jeweils maßgebenden Bogenlänge immer 
durch einen Akzent bezeichnet werden, so daß 


‚_.de ‚_dy 

ande ale FE 
aber 

‚ad dy 


und infolgedessen — der Punkt bedeute die Ableitung nach der Zeit t — 
=), eye) 
| era), Hey 
gilt. Mit Rücksicht auf 2? +9”? —=1 
folgt aus (8) 


und 


csp— er + y'y x | ET N NT Te 9); 
sinp=—ay+yr] 


und indem man diese Werte i in (7) einsetzt und nach den Koordinaten des Fronsle Hons punkte 


auflöst, 
= as) — As) (Hy) + HA) ey + Yr), | 
Yy=yl)—rA) (ey +YE) ya) tYyYy) | 


Die Winkelgeschwindigkeit = 9 erhält man aus Formel au und den durch Ableitung dieser 
Gleichungen entstehenden Formeln 


— sn 99 = (ar + y”’Yy)s ne (Er yy)As, 
csgP=(-EyY Hy) + ey HE), 
und zwar ergibt sich nach einfacher Zwischenrechnung 
ol Ya NEN): 2:0: 11). 
Zu demselben Ergebnis gelangt man fast noch einfacher durch folgende en. Betrach- 


tung: o ist die Differenz der Tangentenwinkel 9 von C und d — go von G,alog = 9 — (d —p), 


und daher 
I TE ae 


Dam ist nun m Em Ez — y'x' = K, die Krimmung der Kurve C, und entsprechend stellt 


ds 
dd. IE 9) _ r' Pr A 
die Krümmung von € dar, so _ daß man wieder (11) oder in anderer Form 
@=älK-AS8) nu ne: TER EUR (12) 


erhält. £ 


Durch die Formeln (10) und (11) oder Se sind nun &, Y, und & mit Hilfe der Gleichungen 
der Kurve ( und & dare et 


3. Spurkurve und Rollkurve 


Um die beiden Pölkurven des Rollgleitens zu bestimmen, hat man 2, %o und & aus (10) 
und (12) in (4) und (5) einzusetzen. Wenn man bedenkt, daß in diesen beiden Formeln die Ab- 
Kir 


= 
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leitungen von % und %, nach der Zeit t vorkommen, wird man zunächst or eine Außerst, lang- 
wierige Rechnung und ein unübersichtliches Ergebnis gefaßt sein; aber überraschenderweise ist. 
das Endresultat erstaunlich einfach. 


Man berechne zunächst 


en H=a— ai tyyIHHleytyT)o. e 


Damit ergibt sich mit Rücksicht auf (12) 
ED. U Bye) SERIEN 

-E@E HN) Ey ty) tre@d + YH Ky TEE Hy} 

—=20—sy'(l —A) 


3% 2 

ö | sy ara HE) Alert + YyY) Ä 

e n — 1 Dun I y'r 3 N. ER +Ay' 2.) 2 & 

3 N. | —y@ "+ Y'y +Au ri 9} jr 
: und ebenso | 

mMo—yotan=y—e—y ty Het) 

| + fa! ar (a Hy) Ay a Hy Er) 4 

3 rear tr yy + Anz 7 +Ay' 0) 2 k 

of + y(— xy’ + y"x , Er D +Ay v”)} : 

und eliminiere rechter Hand die zweite Ableitung der Koordinaten 2, Y, E y, mit ul der Glei- : 

chungen 3 

"=—y'K, y'=+w«K, E 

"’=—yß8, er a 

die sofort aus i 

ER Bel) ER a U Yy% ER art yo B = 

Eye, very "= E40 | 

folgen. . | 


und / Ar N, 
yo—=yo+sa(1l—A)+s(K—AR)- EV OR 5 
EYE) HREHF N) HE HN HH HN 
=yo+s«(1—J), | 
schließlich unter der Voraussetzung 
0=s(K—AR) +0 


als Parameterdarstellung der raumfesten Polkurve oder Spurkurve BIRNEN k = g 


r 
1 
vn Pe \ r 
Fe ; . b 
1 A Se le A ar ten 5 ae a de 


vr ol ee 

En = TREE RUE AIR NT NT N 
ee en 

Bi EROELIE AR 2 RENTEN, 2, 


wie man sieht, sind diese überraschend einfachen Formeln den allgemeinen. ee 3 
raumfesten Polkurve sehr ähnlich. Um nun auch die körperfeste Polkurve zu bestimme 
&, und %, aus 1) berechnet; das ist SEE als eg zu Keane 


und IndeHR man diese Werte in (5) einsetzt, 


et 
» 


n=,@ sin —y Co 


EI NET Re 
n= " (@ 005 P— y'sit 
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Daraus folgt mit Rücksicht auf (8) und (12) und o=g9=g's 


on uk) 3 


4 eo REN HEMER MAN, 
| 


> 
k 
. 
das ist die Parameterdarstellung der körperfesten Polkurve oder Rollkurve; sie stimmt mutatis 
\ mutandis wörtlich mit der raumfesten der Polkurve überein. 
j 
1 4. Beispiele 

Zunächst die Grenzfälle reinen Gleitens (A= 0) und reinen Rollens (A=1). Für A=0 
liefert (13) als raumfeste Polkurve die Evolute E von Ü: 


y'(s) EN NS) 
K(s) ’ Ym = IT Re 
Bei der Berechnung der körperfesten Polkurve ist zu beachten, daß das Argument von 


% 4, 2, und ft die Bogenlänge 3 = Asist, also hier 8— 0, und daher folgt aus (14) mit Rück- 
sicht auf (3) 


0 u 
>m 9 m K (s) 


Im—% (s) Se 


Die Rollkurve ist also, da ja s veränderlich ist, ein Stück der Normalen von CE im Punkt 3=0. 
Die Normale von €, die ja zugleich immner Normale von ( ist, rollt ohne zu gleiten an der Evolute 


von Ü ab — dieser aus der Differentialgeometrie bekannte Satz erscheint hier kinematisch 
| gedeutet. 
Im Falle des reinen Rollens (A= 1) wird nach (13) und (14) 
m = %(8), Ym = Y(8) 
und 
Em nn £(s) ’ Nm = y6s), 


d.h. die raumfeste Polkurve ist mit C’ identisch und die körperfeste Polkurve mit &, was ja von 
vornherein klar ist. | | 
Die zu dem allgemeinen Fall des Rollgleitens (A + 0 und 1) gehörige raumfeste Polkurve (13) 
ist Element einer Kurvenschar, der einerseits die Kurve (C und andererseits ihre Evolute Z als 
Grenzkurven angehören ; ebenso erscheint die körperfeste Polkurve in eine Kurvenschar (14) ein- 
 gebettet, der die Kurve € und ihre Normale im Anfangspunkt 3 — 0 angehören. Dieser irgend- 
wie überraschende Sachverhalt: daß die Kurve C und ihre Evolute als Kurven ein und der- 
selben Schar erscheinen (und ebenso die Kurve & und ihre Normale im Anfangspunkt), mag 
=z noch an dem Beispiel der Kettenlinie aufgezeigt werden, an der die Tangente rollgleitet. 
Hier ist x 
Or Sing, y=y1+8, K = —— 
und ergibt sich 


= U Sins— l—Asy1+2 
:. DER Ym = (2 5, 4) +3» 


und 
Ems, Nm — (il KEN 4) (1 ni 8°) 
s. Bild 2a und b. 

Endlich zeigt Bild 3a, b noch die beiden Scharen der Polkurven, die zum Rollgleiten eines 
Kreises & auf einem festen Kreis ( gehören. In diesem Fall entarten allerdings die Evolute von € 
und der Abschnitt der Normalen von € zu einem koinzidierenden Punkt. 

'Es liegt nahe zu fragen, ob sich die raumfeste und die körperfeste Kurve des Beschleuni- 
gungspoles im Falle des Rollgleitens auch in so einfacher Weise durch die Rollgleitkurve € und 
Schiene (© darstellen lassen, wie die beiden Kurven des Geschwindigkeitspoles. Das ist aber nicht 
der Fall, wie schon das einfache Beispiel des auf einer geraden Schiene rollenden Kreises lehrt. 
Es zeigt sich nämlich, daß die beiden Kurven des Beschleunigungspoles durch die Kurven € 

_ und Callein noch nicht bestimmt sind, vielmehr bedarf es dazu noch der Kenntnis des Bewegungs- 
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ablaufes. Die bekann- 
ten Gleichungen der Be- 
schleunigungspolkurven 
enthalten im Zähler und 
Nenner quadratische For- 
men der Veränderlichen 
© und ®, und die Ko- 
effizienten dieser beiden 
Formen sind in der Regel 
nicht zueinander propor- 
tional, so daß zur Bestim- 
mungder Beschleunigungs- 
polkurven die Kenntnis 
von &(t) unerläßlich ist, 
also der zeitliche Ablauf 
des Rollgleitens. Darum 
hat die Berechnung der 
Beschleunigungspolkurve 
kaum Interesse. 


5. Zusammenfassung 


Die bewegliche ebene 
Kurve € mag sich längs 
der raumfesten Kurve oder 
Schiene (so bewegen, daß 
die beiden Kurven stets 
gemeinsame Tangenten be- 
sitzen. Ferner sei das Ver- 
hältnis der beiden Bogen- 
längen, die längs der einen 
und der anderen Kurve 
zwischen ihren zur An- 
fangszeit =0 und zu 
einer beliebigen Zeit £ ge- 
meinsamen Berührungs- 
punkten gemessen werden, 
ein von der Zeit unab- 
hängiger Wert A. Die 
durch diese Voraussetzun- 
gen gekennzeichnete Be- 
wegung wird als Rollglei- 
ten bezeichnet, der Zah- 
lenwert A ist das Maß für 
das im Rollgleiten ent- 
haltene Rollen und Glei- 
ten: A= 0 bedeutet reines 
Gleiten, A= 1 reines Rol- 
len. Die raumfeste und 
die körperfeste Polkurve 
.diesses _Bewegungstypus 
werden allgemein mit 
Hilfe der Gleichungen 
von C und €& dargestellt 
und in besonderen Fällen 
diskutiert. 


Eingegangen am 28. Mai 1950. 
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Bild 3. Raumfeste (a) und körperfeste (b) Polkurve beim Rollgleiten eines Kreises & 
am Kreis C 


os 
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Bemerkungen zur Entwicklung gegebener Funktionen 
nach Legendreschen Polynomen 
mit Anwendung auf die numerische Integration 
gewöhnlicher linearer Differentialgleichungen 
Von Erwin Fehlberg in Frankfurt (Main) 


Teil I: Zur Berechnung der Koeffizienten, die bei. der Entwicklung gegebener Funktionen nach Le- 
yendreschen Polynomen auftreten, wird die zu entwickelnde Funktion abschnittsweise durch Parabeln 
dter Ordnung ersetzt und das Produkt dieser Näherungsparabeln mit den entsprechenden Legendreschen 
Polynomen exakt integriert. 

Teil II: Für die gewöhnliche lineare Differentialgleichung 2ter Ordnung wird die erste Randwertaufgabe 
numerisch gelöst, indem die Koeffizienten der Differentialgleichung nach Legendreschen Polynomen ent- 
wickelt werden und für die Lösung eine nach Legendreschen Polynomen fortschreitende Reihe mit un- 
bestimmten Koeffizienten angesetzt wird. Die Berechnung der allgemeinen Lösung wird so auf die Auf- 
lösung eines linearen Gleichungensystems zurückgeführt. Die spezielle Lösung der Randwertaufgabe kann 
numerisch so angegeben werden, daß darin die rechte Seite der Differentialgleichung sowie die Randwerte 
noch völlig offen bleiben. 


Part I: If a given function shall be ewpanded in Legendre's polynomials, the author evaluates the 
coefficients by replacing, in parts, the function by parabolas of the fourth order, and ewactly integrating 
ihe products of Ihe parabolas of approwimation and the corresponding Legendre’s polynomials. 

Part II: The method is applied to Ihe solution of Ihe first boundary problem of the ordinary linear 
differential equation of Ihe second order. The coefficients of the differential equation are expanded in 
Legendre's polynomials, and the solution is written as a series with indeterminate coefficents proceeding 
in Legendre's polynomials. The general solution is obtained by solving a set of linear equations. The 
special solution can be given numerically without fizing the right hand side of the differential equation and 
the boundary values. 


Premiere partie: Pour &valuer les coefficients, qui se prösentent en developpant une fonction donnee 
par polynömes de Legendre, cette fonction esi remplacee, par sections, par des paraboles du quatri&me ordre. 
Le produit de cettes paraboles d’approximation et des polynömes correspondants de Legendre est integre 
exactement. 

Deuxiöme partie: Le premier probleme des valeurs bordees est resolu numeriquement pour une Equation 
differentielle ordinaire du second ordre; les coefficients de l’equation differentielle sont developpes par 
polynömes de Legendre et la solution est Ecrite sous la forme d’une serie en polynömes de Legendre. 
L’&valuation de la solution generale est tellement reduite sur le probl&öme d’un syst&me d’equations lintaires. 
La solution sp£ciale du probleme des valeurs bordees peut ätre donnde numeriquement de manidre, que le 
cötE droit de l’öquation differentielle et les valeurs bordees sont encore libres. 


Yacrs I: lust Bsruncnenna KoahunmeHToB, BXONAIUNX B Pasıokenne naumoli Pyukunn, Ha 
NOIHHOM»L Jleskanıpa, nannası PyHRUmA BAMEHSIETCH B OTMENBHEIX YYAcTKaX TApaboNamu 
YETBEPTOTO TOPAAKA, A 3ATEM IPOHSBENECHHE ITAX IPHONMSKEHHEIX TTAPabo.ı Ha COOTBETCTBEHHBE 
HOJIUHOMBL Jleikanapa TOYHO UHTerpnpyerca. 

Yacrp Il: Iror meroy npumenserca npu pemeHun mepBoi kpaeBoi sanaum, OTHoCAIMeHcH Re: 
K IPOCTBIM JIHHCHHBIM AubepeHmmanbEsIM YPaBHeHHAM BTOPoro nopsra. ls 3Toro Koahn- « 
IMeHTBI Aubepennmansuoro ypaBHeHnns pastaramrca Ha. moımHoMsIr Jleskanımpa a perenme 

" IPUHUMAETCH B BUNG PIAa TOAMHOMOB Jleskanıpa c HeHSBecTHEIMN Koshuıumenramu. Tarum 
o6pasom oÖlNee pemmenme MAHHON IPOÖMeMEI CBOAUTCHA K PenıeHnnm JIMHeiHoli CMCTEMEI 
ypasHenni. Macruoe pemenne KpaeBoH Banayı MO3KeT ÖBITb UMCHEHHO NAHO B TAakOK_hopMe, 
B KOTOPoH mpaBası CTOPoHa AndepenmmassHoro ypaBHeHnA u Bra Rue BHAYCHHA OCTAITCH Elle 
cOÖepIIeHHO OTKPEITEIMH. 


I. Zur Entwicklung gegebener Funktionen nach Legendreschen Polynomen 


Entwickelt man eine Funktion f(r), die auf dem Intervall Be £ z.B. stetig 
sei, in eine nach Le gendreschen Polynomen fortschreitende Reihe: 


Fe) a> 7 D, (7) 


so sind die Koeffizienten a, der Entwicklung bekanntlich: 


2n+1 Be | 
el, " SO Pu@)ar,  n=01 2 302 
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In einer früheren Arbeit!) wurde ein Verfahren zur numerischen Berechnung der «a, an- 
gegeben, das darauf beruhte, daß man das Intervall —1<r<-+1 in Teilintervalle der 
Breite Ar—=0,1 einteilt und für je zwei aufeinander folgende Teilintervalle die Funktion f(r) 
durch eine Parabel 2. Ordnung ersetzt, die in den drei Endpunkten z;, 75:1, Tx.ı2 der beiden 
Teilintervalle mit den dortigen Funktionswerten fr, fr-+1, fs+2 übereinstimmt. 

Für solche Parabeln an Stelle von f(r) läßt sich das Integral in (2) exakt ausführen. Dieses 
Verfahren zur Berechnung der Koeffizienten erreicht eine relativ hohe Genauigkeit dadurch, 
daß nur f(r) angenähert wird, während die vor allem für großes n stark schwankenden Legendre- 
schen Polynome P,(r) in ihrem genauen Funktionsverlauf erfaßt werden. 

Wenn nun die Funktion f (r) ihrerseits auf dem Intervall —1<r<s-+1 stark schwankt, 
kann es erforderlich werden, den Funktionsverlauf durch Parabeln von höherer als 2. Ordnung 
anzunähern. 

Das in der oben zitierten Arbeit für Parabeln 2. Ordnung angegebene Berechnungsverfahren 
soll daher im folgenden auf Parabeln 4. Ordnung erweitert werden. 

Die an den Stellen Tas Te+1/2s Tet1> Tar2/25 Te+2 (Tr+ı —_>= Are 0,1) mit den 
dortigen Funktionswerten von /(r) übereinstimmende Parabel 4. Ordnung f,(r) hat die Gleichung: 


hAaa=ht+ 7 —_ (Kgyndy 3 Reue rue! (T— u) (T— Term) 


(de) 
2 


= 
Fade (m) (tern) ter BR 
2 
Tr ee - 7) (T— Ha) TH) (Tea) 
ie (4) | 


wenn: fu) = fr, MrTr+ 2) — It) = Afr,; Afksım — Afı = A®f, usw. gesetzt wird. 
Einsetzen von (3) in (2) ergibt: 


Tk+2 76 +2 


In = nn 2 fe: i' SET er (T— 7) ) Pulr)dr 
. ? T]. 
; er? 
+ nn : f (T u (7 — Tr41je) 2, (7) dr 
v 215) e 
A "%+2 \ 
4 a . (T— %)(T— Tr) (T— Te) Pu(r) dr 
(5) RR 


7642 


2 RER ih F— %) Tre) (TH) (T — Tara) Pulddrz . 


I 
N J 
Setzt man abkürzend: 


7.+ +2 . T_L 


fe 1 EEE (kr Re jr. f P.(z)dr, 


5 SE, 76-42) 


m 


‚sr @- ” B „(ddr= are) ee 17 (re — 7%) Pı(r)dr, 


LEE 0 F(k+2) 


er 2): . angew. "Math. Mech.: 24 (1944), S S. 71-76: „Eine Bemerkung zur numerischen Differentiation 


15 u reg ausgeführt am Beispiel der Kugelfunktionen als Approximationsfunktionen.““ 


dr 
+ 


% 
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7.19 EE 


k 
[ @—H) C— Hp) Pa d)dT=L, (% %R42)=(—1)*: [ GT) (T— TH) P(r)dr, 
T, T 


— (k+2) 
Tk+2 
3) N: Te+ıy2) (T— Te+1) P,@)dr=M, (T,, Te+2) 
% 
= (IH. [ FT) (T— Te419) (Tr — Tan) Pu(r)dr 
T_(£+2) 
"£+2 


J (7 Tr Tı) (7 E Te+1j2) (7 x Tk+1) (7 m TR+ 3/2) I, (d)dr= er N,.(T 748) 


f.. (ET) (T—- T_a41 9) (7 — T_%+1,) (T— T_(#4312)) Pn(T)AT, 
—(k-+2) 


hr. 


so findet man wegen: P,(r) = 


On+1): 
(2n--1)- 


Pam (®)] 2 
Pu (7.)] ’ } 


1 [Par (rd) — 


Ju (%. Te+2) = [Pa+ı (TRr2) — Pa-ı (Tr+e)] — [Pr+ı CR — 
K,(%, Te+2) —=2-Ar- [Pr+ı(Ta+2) = Pa (T«+2)] 
= Yn+ılTz, Te+2) — In (Ts T«+2)] P 
Daltps Tp+2) = (Ar)? [Pn+1(Ta+2) — Pn-ı (Tx+2)] 
1 
+ 2) "Ar: Yn+ı (7 s Te +2) — Ja-ı (7, Tr+2)] —2* [Ka+ı (T.. Te+2) 
— Ki (7. » Te+e)] 
M,„(t; T.+2) =3° (A r)® : [Pa+ı(Te+2) EEE Po (Tr+2)] 
1 3 
9, (Ar) - Insı (7, Tar2) —In-ı (Ts %a+2)] + 2, ? Ar -[Knrı (7.72) 
— Knı (Ts 2 —)3. [Ia+ı (7; ’ Te+2) — Lu (7. D Te-+2)] , 
(2 n-+ 133 N, En ’ 1)=, 2, (A zZ) [Pr+ı(Ta+2) — Pa (Tr+2)] 
“4 (Ar) - an (Tas Ta+2) — In-ı (Ts Tar2)] — 2° (AT)? [Kası (%r, Te+2) 
n—1 (Tg, T42)] + 3° AT [Imzı (%, Ta+2) — Dn-ı (% Tere)] 
— 4 [Marı (7 » Te-+ 2) — M,-ı (7. » Tkra)] 


Die Zahlenwerte dieser Funktionen sind in den Tabellen 1 bis 5 zusammengestellt. Die 
Tabellen 1 und 2 stimmen mit den entsprechenden Tabellen der oben zitierten Arbeit überein 
und sind hier nur der Vollständigkeit halber nochmals aufgeführt. [Die Funktion L,(T,, Tx+2) 
hat jetzt eine etwas andere Bedeutung als die Funktion Ly,r+1,n (T4, Tx+2) der oben zitierten 
Arbeit; sie würde der Funktion Zy,21172,n(T% Tk+2) entsprechen]. 


(2n-+-1)- 


(2n-+1)- (4). 


Tabellel 
(2n-+]) 3 (T.; 0 (—1)r- (2n +1): J,(T_@+2)> T_g) 
> (7. > T, +2) 
(0; 0,2) (0,2; 0,4) (0,45 0,6) (0,6; 0,8) (0,8; 1) 
R N 
0 -+-0,200 000 000 —-0,200 000 000 40,200 000 000 +-0,200 000 000 -+-0,200 000 000 
1 0,060 000 000 —+0,180 000 000 —+-0,300 000 000 +0,420 000 000 | 0,540 000 000 
2 —0,480 000 000 —0,360 000 000 —0,120 000 000 +0,240 000 000 +0,720 000 000 
3 —0,203 000 000 — 0,525 000 000 —0,595 000 000 —0,245 000 000 +0,693 000 000 
4 -+0,587 520 000 +0,123 120 000 —0,503 280 000 . —0,686 880 000 +0,479 520 000 
5 0,374. 924.000 0,718 212 000 0,174 460 000 — 0,738 892 000 0,158 796 000 
6 —0,601 036 800 -+0,315 806 400 +-0,760 468 800 —0,315 369 600 —0,159 868 800 
7 —0,544 926 300 —0,600 646 500 +0,599 878 500 +0,334 897 500 —0,375 140 700 
8 +0,539 473 920 —0,713 647 080 —0,194 528.280 +0,796 207 920 —0,427.506 480 
9 +0,688 168 076 +0,195 201 538 — 0,819 840 310 +0,773 187 042 —0,317 185.096 
10 —0,420 272 456 +0,854 589 083 _ —0,625 406 332 +0,291 447 859 —0,100 358 154 
11 —0,785 888 617 +.0,316 858 074 0,191 599 461 —0,328 336. 783 | +0,134 088 177 
12 +0,264 015 299 —0,659 703 344 —+0,794 811624 pp 697 316 154 


+0,298 192 575 
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Tabelle 2 
Zn HYDR, (4) = (Urt 2n+ 1) K,(T_a40) Te) 
(Ti %x+2) 
(0; 0,2) (0,2; 0,4) (0,4; 0,6) (0,6; 0,8) (0,8; 1) 
0 +0,020 000 000 -+0,020 000 000 —+.0,020 000 000 —+-0,020 000 000 0,020 000 000 
1 0,008 000 000 +0,020 000 000 —+-0,032 000 000 —+-0,044 000 000 —+-0,056 000 000 
2 —0,047 000 000 —0,033 000 000 —0,007 000 000 —+-0,031 000 000 —+-0,081 000 000 
s, —0,026 880 000 —0,056 280 000 —0,057 680 000 —0,014 280 000 —+0,090 720 000 
E: +0,054 420 000 +-.0,001 836 000 —0,059 388 000 —0,063 732 000 -+0,084 564 000 
5 +0,048 998 400 +0,070 514 400 —+-0,001 601 600 —0,083 820 000 —+-0,065 577 600 
6 —0,049 794 940 +0,048 289 020 —+-0,070 915 780 — 0,057 524 740 0,039 445 380 
7 —0,069 765 120 —0,045 655 020 —+-0,077 989 080 —+-0,003 933 180 —+-0,012 849 840 
8 +0,035 346 992 —0,085 151 497 0,009 401 102 —+.0,067 133 741 —0,008 315 428 
9 +.0,085 474 099 —0,009 906 285 —0,072 861 975 +0,096 394 310 —0,020 368 480 
10 —0,013 666 049 +.0,083 360 652 — 0,089 697 842 —+0,075 041 198 —0,022 523 873 
11 —0,093 415 372 -+0,067 612 681 —0,022 723 998 —+0,014 953 480 —0,016 706 663 
12 —0,011 984 784 —0,039 459 807 —+-0,066 289 781 —0,050 985 629 —0,006 629 003 


Für die Koeffizienten «a, ergibt sich dann: 
1 <0,1> 


DAR Sen) In le 10) 


k 


4, — 
45 ne -[Af. + (—1)"- Af-ı] -2n+1):K,(%, Terre) 
2 
a A DIOR ER DENT ERS GE 


+ gas [A HD Aa en HD) Mu lm Turn) 
Be Are) AR Fl) N. lt: Tre)! 


Hierbei bedeutet f_, den Funktionswert an der negativen Stelle 74, Af + =f_ «-+1) Se 
Ar Ale SEITEN USW. R 
Bei Benutzung der Tabellen 1 bis 5 kann die Berechnung der Koeffizienten a, bequem 


mit der Rechenmaschine durchgeführt werden. 
Das folgende Zahlenbeispiel zeigt die durch die NE von Parabeln 4. Ordnung 


erzielte Erhöhung der Genauigkeit der Koellizienten An» 
Beispiel: f(r)= Ps(r) für -I<ts+1 und Ar=Nl. 


un) 4, a4 Ag Ag do da 


Parabeln 2. Ordnung. . -+0,00069 -+0,00789 -+0,02219 -+-1,02402 +-0,01309 —0,01251 —0,04044 
4. eB : 0,00000 -+-0,00001 -+-0,00002 -+-1,00003 --0,00002 +-0,00001 0,00000 
Exakte Wertes? °: 0 0 0 1 0 0 : 0 
Tabelle 3 
(22 +1)2 (7, %ı)=(- 1». (2n +1): L, (T_+2)> T_,) 
(Ts T.4e) | 
in (0; 0,2) (0,2; 0,4) (0,4; 0,6) (0,6; 0,8) (0,8; 1) 
: | 

0 0,001 666 667 0,001 666 667 -+0,001 666 667 +0,001 666 667 0,001 666 667 
12 +.0,000 800 000 +.0,001 800 000 —+0,002 800 000 +.0,003 800 000 +-0,004 800 000 
Fr 2 —0,003 836 667 —0,002 536 667 —0,000 236 667 0,003 063 333 +-0,007 363 333 
\ 3 —0,002 669 333 —0,004 914 000 —0,004 638 667 —0,000 443 333 +-0,009 072 000 
4 —+0,004 191 000 —0,000 567 000 —0,005 507 400 —0,004 834 200 —+-0,009 768 600 
« 5 +0,004 801 133 +.0,005 678 640 —0,000 983 693 —0,007 473 547 —+0,009 444 600 
6 —0,003 340 909 +0,005 084 469 —+-0,005 356 403 —0,006 493 435 —+-0,008 237 619 
7 —0,006 690 526 —0,002 653 758 +0,007 514 382 —0,001 939 078 —+-0,006 404 314 
3.8 -+0,001 512 318 —0,007 753 921 —+-0,002 850 113 0,004 121 227 -+0,004 273 134 
9 +0,007 931 766 —0,002 908 084 —0,004 948 316 0,008 618 995 0,002 185 269 
+0,000 992 754 +.0,006 282 503 —0,008 852 378 +-0,009 108 734 +-0,000 434 896 
—0,008 240 335 +.0,007 821 226 —0,004 939 690 +0,005 114 430 — 0,000 779 477 
—0,003 786,691 —0,000 892 939 -+-0,003 693 006 — 0,001 542 040 —0,001 382 712 


“ 


Z. angew. Math. Mech, 


108 Fehlberg, Bemerkungen zur Entwicklung Bakebener Fonlktionee Bd. 31 Nr. 4/5 April/Mai 1951 
Tabelle 4. 
(Zu +1) M n (7, 7.12) Er; (her? : M, (T_«+9 3 T_,) 
(Fr » %.19) 
BEN (0; 0,2) (0,2; 0,4) (0,4; 0,6) (0,6; 0,8) (0,8; 1) 
Pe 
0 +.0,000 100 000 0,000 100 000 0,000 100 000 —+0,000 100 000 +-0,000 100 000 
1 +0,000 052 000 +.0,000 112 000 0,000 172 000 —+-0,000 232 000 —+-0,000 292 000 
2 —0,000 227 000 —0,000 145 000 — 0,000 003 000 —+.0,000 199.000 +-0,000 461 000 
3 —0,000 172 400 —0,000 299 800 — 0,000 270 400 — 0,000 000 200 —+-0,000 594 800 
4+ +-0,000 237 900 —0,000 056 700 -—0,000 344 100 — 0,000 266 460 —+-0,000 685 260 
5 —+.0,000 306 262 —+0,000 328 086 — 0,000 095 634 —0,000 451 154 0,000 728 838 
6 —0,000 169 354 —+0,000 332 746 —+-0,000 291 950 — 0,000 435 979 —.0,000 726 422 
7 — 0,000 418 220 —0,000 114 386 0,000 468 256 —0,000 199 709 0,000 682 645 
8 --0,000 037 640 —0,000 467 741 0,000 238 675 —+-0,000 163 648 0,000 604 941 
) 0,000 480 202 —0,000 240 924 —0,000 232 377 —+-0,000 486 021 0,000 502 599 
10 0,000 134 309 -+0,000 325 265 —0,000 540 648 -+0,000 606 232 —+-0,000 385 941 
11 —0,000 473 564 +-0,000 515 829 —0,000 394 471 —+-0,000 450 961 —+-0,000 265 631 
12 —0,000 317 274 +-0,000 051 088 +0,000 110 773 —+0,000 075 514 


Tabelle 5. 


—+0,000 151 998 


(Zn +1) N, N I 5 Tor 5) a al 162 N, 638 (k+ 2) > T_,) 

(0; 0,2) (0,2; 0,4) (0,4; 0,6) (0,6; 0,8) (0,8; 1) 
0,000 002 333 0,000 002 333 0,000 002 333 —+-0,000 002 333 —+-0,000 002 333 
-+0,000 001 400 0,000 002 800 0,000 004 200 0,000 005 600 —+-0,000 007 000 
—0,000 005 169 —0,000 003 069 0,000 000 431 —+-0,000 005 331 —+-0,000 011 631 
—0,000 004 607 —0,000 007 270 —0,000 006 013 0,000 001 123 —+-0,000 016 100 
-+-0,000 005 011 —0,000 002 339 —0,000 008 729 —0,000 005 339 —-0,000 020 179 
+0,000 008 063 -++-0,000 007 245 —0,000 003 822 —0,000 010 889 —+-0,000 023 576 
—0,000 002 714 0,000 009 115 0,000 005 782 —0,000 012 298 0,000 026 014 
—0,000 010 744, —0,000 000 908 -+0,000 011 981 —0,000 007 965 —+-0,000 027 315. 
—0,000 001 232 —0,000 011 418 -+0,000 008 392 —+-0,000 000 973 —+-0,000 027 457 
+0,000 011 853 —0,000 008 416 —(0,000 003 287 —+-0,000 010 803 —+-0,000 026 576 
+0,000 006 104 +-0,000 006 056 —0,000 013 608 —+-0,000 016 804 0,000 024 922 
—0,000 010 904 +0,000 014 359 —0,000 013 066 —+0,000 015 576 —+-0,000 022 769 
—0,000 010 992 --0,000 004 842 — 0,000 000 885 0,000 006 948 +-0,000 020 327 


Die in der oben zitierten Arbeit angestellte Fehlerbetrachtung läßt sich in ähnlicher Weise 
auf die Approximation von f(r) durch Parabeln 4. Ordnung übertragen. 
Es sei noch darauf hingewiesen, daß das hier entwickelte Verfahren zur Berechnung der. 
Entwicklungskoeffizienten «a, nicht auf eine Entwicklung nach Legendreschen Polynomen 


beschränkt ist, sondern in gleicher Weise für die Entwicklung nach Funktionen eines beliebigen 
Orthogonalsystems anwendbar ist. 


PR 


Su die Differentialgleichung (1) sei die erste Randwertaufgabe: 


gestellt?). 
Gesucht sei die Lösung von (1) und (2) auf ı dem Intervall —1< :Ss+1- 


®) Auf die Beseitigung des Gliedes mit y’ sei hier verzichtet, da diese zur Folge hätke, daß die sich dann Kr 
b(2)-y=%(a) in b(x) den Differentialquotienten: a’(x) enthalten 


würde, d.h. 5(z) bei evtl. nur in Kurvenform vorgegebenen Funktionen a(®) und b(x) nur auf dem Wege 
über eine Differentiation gebildet werden könnte. fi 


ergebende Differentialgleichung y’’ + 


Be ®) Es sei darauf hingewiesen, daß die folgenden Ausführungen sich in völlig dere Band Ders auf gewöhnliche ‘ 
An lineare Differentialgleichungen von höherer als zweiter Ordnung sowie auf andere u EEE Aula 
Er aufgaben übertragen. } 

be ; 


= —1, 


y=y(—I1), 


=-+J1, 


y=y(+1) 


II. Anwendung auf die numerische Integration gewöhnlicher linearer Differentialgleichungen 
Die folgenden Ausführungen beziehen sich auf die lineare Differentialgleichung 2. Ordnung: 
y" + ala) yY’+bla)y= ecke) - 


in der a(x), b(x) und c(x) stetige Funktionen von & seien?) ?). 


A), 


.@&) 


üben 385 Ban ra sr I ManklEe 27 ERS re Er 
„far 1 BB Se ER RR 


n 
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Die Funktionen «a(x), b(x) und c(x) entwickele man auf diesem Intervall nach Legendre- 
| schen Polynomen: 
$ a a (Bee bir) =D Baer .te), cl) Ir De 2.) 
N ' n m 
Für die gesuchte Lösung mache man ebenfalls einen nach Legendreschen Polynomen 
fortschreitenden Ansatz: ; 


HKS Pl) nee) 


Mit (3) und (4) gehe man in (1)ein. In den Entwicklungen (3) und (4) möge n von O bis 7 
laufen, da dies im allgemeinen zu einer genügend genauen Lösung von (1) führt. Außer Gliedern 
mit P,(&) treten dann in (1) Glieder mit P,„(x)- P,(&), Pn(x) Pn(x) und P}(x) auf, die man 
sämtlich wieder nach P,(x) entwickeln kann: 

a) P„(®)- P„(x): Aus der Rekursionsformel 

u 2m—]1 m—1 
ale) P,(«) ep (&) ee 


mM 


 Pmze (x) 


der Legendreschen Polynome gewinnt man: 


Pula): Par. TI. pn lo): Parıla) 


m 2m 1 
> 2m—1 N m— 1 () 
a FR y 9 nt 1 2 Be (x) Pr (%) =; Mn ? Pr 38): PR (x) 
als eine für unsere Zwecke geeignete Rekursionsformel. 
N n+1 x 
| MierR,Le= Pound’ 2,2, 221 ent 2n+i - P„;ı lassen sich für m>2 alle 


Produkte P,,: P, nach (5) rekursiv durch die Legendresche Polynome ausdrücken. 
b) P.(&) Pa(x): Au P,(e)=(2n— 1) Pu_ı(2) + (2% —5) Pus(&) + (2 n— 9) Pu (%) 
+(2n — 13) Pa-1(x) *) (#7). "folgt: 
P(&) Pa(2)= (2% — 1) Pn(&) Pa-ı(8) + (2n—5) Pu(®) Pn-3(%) 
+(2%n— 9) Pn(&) Pa-5(2) + (2 rn — 13) Pn(&) Pa-(&) 9) 


als eine für n<7 in Verbindung mit (5) geeignete Rekursionsformel. 


o) Pal): Pila)= (an 1) (an —3) Para) +2: 2m —3)(2n— 7) Paula) | (7) 
0 +3@n—5)An—11) Ps) Men] 


Nach Eingang mit den Ansätzen (3) und (4) in die Differentialgleichung (1) und Einführung 
der Entwicklungen (5), (6) und (7) kann man zu einem linearen Gleichungensystem für. die 
gesuchten Koeffizienten Y, der Lösung gelangen, indem man einen Koeffizientenvergleich für 
die einzelnen Legendreschen Polynome durchführt. 

Der Koeffizientenvergleich für P,(x), Pı(x), Ps(x), P;(x), P,(x) und P,(x) ergibt ein 
Gleichungensystem, das sich in Matrizenform?) folgendermaßen schreibt: 


(6) 


r ET. ER we (8) 
mit: 
C, Y, do Aoı °  * dor 
C, Y do di dir 
= , y == | ’ und D m 

C, Y, \a,. dene dar 
Die Koeffizientenmatrix D läßt sich additiv zusammensetzen aus drei Matrizen: 
ß ee. 22222...) 
E von denen 3 = (75° &un  Ynur-aus: festen Zahlwerten besteht, während U = (** "%yn' ) und 
y B=(°--Pyan*:) von den Entwicklungskoeffizienten der Funktion a(x) bzw. b(x) abhängen. 
; 4) Glieder mit P,(z), in denen v negativ ist, sind wegzulassen. 


: 5) Die Verwendung der Matrizenschreibweise zur textlichen Kürzung geht auf einen mir von Herrn R. Zur- 
 mühl in Darmstadt gemachten Vorschlag zurück, dem ich für seine Anregung Dank schulde. 


\ 


a 


ef Eie?} BR r EN f 
j g ne 7 + at. 107 EN) 
Falts n * f ’ FRA % er 
BEN uw r a 
N ni » ER A ER 
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Ba 


Die Elemente von & sind Null bis auf die folgenden: 


(mm = (2m +3)2m-+-1); mm =2 (2m +5)(2m+1); | ER (10). 
| Im m = 3 (2 m+7)(2m+1) 


Zur Berechnung der beiden anderen Matrizen führt man die ol EN Abkürzungen ein: 
I7(A,)= Sr An ’ 


m—+]1 Mm j 
El) = ya At Ans 
Ss m+2 m+Il, m+ 1 m 3 m m—l 
22m+5 2m-+3 An tom ENDETE Im—3' Ama» 


I | t 3: m +20 ml MN 
LlAn)= 5 2m+7 2m+5 2mt3 "Anıst 9 mt5 2mF+3 2m Pe 

3 m-+]1 m m—1 B) m Mm—L m 2 . 
OWOTTTEL: im oma Dad ee 


Land. mtt.m+3  m42.mHl, 
Ne um LTR AamI Se 
mr mF2 mal) 37m 
2m+7 2m+5 2m+3 2m—1 
RS a | m m—1 


L 2 de Fe 


am ame nee 
Bi Sam Em ge A ee 
2 Dm3 Im Lam 3 Omen 


3 m ' m—|1 Mr } MB 24 
2 -—- 12m 332 ma m Te 


L,(4 ee m+5  m+4 m+3 m-H+2 A N 
m 8 2Amt1i 2m+9 2mt7 2m+5 2mt3 
ee en m+3 m+2 m-+I1 m Ey 
85 2m +92m+72m+52m+32m— 1 
15 m +3 m+2 m+I1 23 m—1 "u 
4 2m +72m+52m +3 2m — 12m —3 mtr 
15 m +2 m+1 .m m—1 mi— 2, 2. RER 
2m 52m FI In ZI 2m an a ae 
3 m+l1 m m—1l1 m—2 ma: ; ; VA 
8 2m+32m—12m—32m— 53m" INT 
7 3 63 m m—1 m—2 m—3 mA ’ ERS 
m 13m 33m 5m 7m gen: | ER 2,6 


amt mt ai, 
mu 2m tB2mtI2mtI2m 7 2mL52mL3 ni a 
63 m+5 m+4 m+3 m+2 m+I m EEE | - 


B: 2m+112m+92mt72mt52mt32m—1 1 = 


105, m+4 m+3 m+2 m4+1 m m—1. EN e 
BAR E 2m+92m+72m+52m+3 2m —12m— a 

25 m+3 m+2 mtl m m-1 m— LP: FE 

4 2m+72m+52m+32m—12m— men DE 

105 , m+2 m+1 m, m—1l.m—2 m—3 An FREEEN 

16 2m+5 2m +3 2m—12m—32m—5 2m— —?” a 
Bach m+1 m m—1i m—2 m—3 mA, T 


“ am Frag mr 
Re Mm Em zl m—2ı m—3 m — —4 
Am-ıom —32m—52m— 72 m — 32 


’ a 3 Wr u aan a ee 
in n ne vr . 
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L,(4 429 m+7 m+6 m+5 m+4 m+3:m+2 m+1 4 
k De 2m+15 2m +13 2m +11 2m +9 2m +7 2m +5 2m +3, "tr? 
Rh. m+6 m+5 m+4 m+3 m+2 m+l m 
16 2m +132 m +11 2m +9 2m +72 m +52 m +32 m—1 
RE m 5 m 4 m+3 m+2 m+Il m m—l 
16 2m +112m +92 m +72 m +52 m +32 m—12m—3 
> are m+3 m+2 m-+I Mm m—1 en 
16 Am 92m + 72m +5 2m +32 m— 12m — 32m 5" 
2 
a 
= 


a 


An +5 


An +3 


175 m+3 m+2 m-+1 m m—l m—2 m—3 


16 2m+72m+52m+3 2m —12m—32m— 5 am, 
189 m +2 m+1 mM m—1l m—2 m—3 m—4 #4 
16 2m 52m 32m — 12m — 3 2m —-52m— 72m —9 ""” 
231 m+1 m m—1l m—2 m—3 m—4 m—m—) A 
o 16 am am 1 Am 32m 5 2m - 72m 92 m— 11.7 
429 m m—l m—2 m—3 m—4 m—I m—6 2 
a nen mom Ar 
Dann wird, wenn man die entsprechenden Abkürzungen in den. B„-Werten mit L,(B,) be- 
zeichnet: 
0 Bmo = Lu(Bn) 
&mı = Lo(An) Pmı = Lı(B,) 
Am? = 3 L,(A„) Bm 2 Tr L;(B,.) 
&m 3 = Lo(An)+ 5 L2(An) Pms = L;(B,.) (11) 
Imst — 3 L,(4,) re 7 L;(A„) j ; Ps = Lı(B) 
ms = Lo(Am) + 5 L,(A4„) +9 Lı(An) Pms = Ls(Bun) 
Am6 — 3 L,(4.) = 7 L;(An) + 11 L;(An) Bm 6 TE Le(Bn) 


Am? ='L,(An) =+ h) L;(An) Ei I L,(A) =: 13 L;(A,) Bm U = L;(B,) 
Durch (9), (10) und (11) ist dann die Matrix ® in (8) bestimmt. 
Zerlegt man die Matrix D in die beiden Teilmatrizen: 


Io Agı ; (gm dog °°° Ir 
io di 12 ds ir dir 
Dı=| le 
IB 14) \dag dys da; 
und den Vektor y in: 
Y, 
Y; 
(7) _ und a SEN 
' Y, 
so wird aus (8): 
De N (12). 
Auflösung nach 9, ergibt (|D,|+0 vorausgesetzt): 
2, dmde— 27'019 
oder wenn man setzt: D,I=6© und. - 19, =7T: 
} ee an: (13). 
2 ’ Die Matrizen S und T mögen die folgenden Elemente haben: 
m ud vo — sn dyı too  Foı 
v0 v0, 


\ 
DIE ? Bi Ne ? & 
—2 = dv dyo — 3 5 dyı r tzı 


an v=() 


j 


"Außer der Entwicklung der Funktionen a(%), b(%) und c(x) nach Legendreschen Polynomen 


+ \ tp nd ’ Re - 2 Be won; 
« \ e ne j VRR ach ARE u 
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Führt man weiter die Matrizen: en I 
0 RE RRTH 0 [ ra \ L, 3 
Dit Are 0 era BER Be; 
BD EEE 305 too kai E 
© = . . und ey = . u Z 
zog er 55 iso n He 


yes Eiyjrane.. ae ee a 
Po(#) . > 
ı Pı(«) j 
Führt man den Vektor: p(x) = “| ein, so kann man die gesuchte Lösung (4) auch in 
_ der folgenden Form schreiben: | 
x YIEp'(E):y ee ee RAR En 


wobei p/(x) den zu p(x) transponierten Vektor bedeuten soll. 

Einsetzen von (14) in (4’) ergibt: PR er 
le mer ar N TR (15) 
oder: 2’(2) 2Y N EI) EPITIS ET Tr Tr . (16). 


Von der Gl. (16) gehen wir aus, um jetzt noch dis Randwertaufgabe (2) zu erfüllen. Setzt man 
s=+1l,y=y(+ 1) und = —1, y—y(—1) in (16) ein, so erhält man zwei Gleichungen für Y 
und Y,, die man in Matrixform folgendermaßen schreiben kann: 
MEER FE re Er 
wobei: ; ; 
(1, tra io 1 Feen "): 
VE a a TH ee sr 


5 5 Ay ... N ‚5 r 
N ER = Ki ; = PR \ ve; t (4 ı) ist 
5 5 5 . : ey & 
2’ (—1)”s, 0 SI (-1N)s,1 ne 3 (—1)s,5 A ) 
v0 v=0 v0 x 
Aus (17) folgt für 9 ee Sn a 
Venus Be 
Einsetzen von (18) in (15) ergibt schließlich für die gesuchte Lösung von iD) und (2): r 
y()= pa) (S-T-MEN) c HT. Mid. 2... Rau) 


ist als wesentliche Arbeit die Berechnung der Matrix D sowie der zu D; 1g2iproken Matrix Dir vr 
durchzuführen. ae. 
Ein Vorteil gegenüber anderen numerischen Lösungsverfahren besteht bei der hier ent- 
wickelten Methode darin, daß die allgemeine Lösung gewonnen wird. Man braucht also, wenn 325 - 
neue Randwerte y(--1) und y(—1) vorgegeben werden, nicht die ganze Andıng von neuem wo 
zu beginnen, sondern hat nur in (19) für r diese neuen Werte einzusetze = 
Desgleichen hat man, wenn sich die rechte Seite c(@) ‘ vor (1) ändert, d.h. der Vektor 
nur diese neuen Werte von c in (19) einzusetzen. h “ TUR, SEN 
Selbst wenn sich die linke Seite von (1) ändert, spart man, falls sich satz Ri 
oder nur b(x) allein ändert, einen Teil der Arbeit bei der Neuberechnung ve Er 


die Matrix VA bzw. B neu zu berechnen braucht. EN r 


e- Teig. hr 


Beispiel: "+8: y+lla+D)y= ze, 
Wegen = Pı(@). und 4(@®+1)=$ Pula) +4 P;( e 
4=0, A=ı L=09I, L=0, A=0, 4 

B=}, B=0, B=} B=0, See 


Math. Mech. 
Ir. 4/5 April [Mai 1951 


Für die Matrizen 3, X, ® und D ergibt sich 
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nach (10), (11) und (9): 


: 00ER. 21. 0 VRR ER N 
| PORN) OHR. 0 E81 105080: 
2 Ban ne „100.20 505 0 
w oa oe 7=10.0:0 307.0 7 
E; 0 05.209.080 06,204.:99,...@ 0.0.0%.0,40 9.0 
= 0.202 000 000,620. 05.183 DDR 0 
LO ROTE 0 EB OA LOL IL 2.00 92200 
0.3. 3,.02105008.00,10 U ERONIEN.0 45.001 84 
Ber 0 00.00 a_|} 9 2 0 4% 0 8 0 
| RR ER SE RT O2 10205170. .70 10.5.7161 
| 00500 0 00% 0 a2 01080 
UBER EL BEE ze ER 0 00, 4 09 Sr 0 14h 


I. 


I 22 Y 
Ja 2402198 R 


Gl. (12) zerfällt in diesem Beispiel in die beiden folgenden Gleichungensysteme: 


2 IR 40 
LET [e EU. 4 


Br de 1 NT, 

, AN, bzw. ö 

a " \ 1825 45 84 Y, 0% 0 12 } 
Bern: se 70 1 LG, 08 25 di 
2 \#5 rer 1ödrin/ \Y, C, 0.0 


mit den Lösungen: 


RL E, 0,296 556 7 —0,082 814 6 
7.) = [0,018 8880 0,032 899 2 
\ 2 0,000. 529 8 —0,001 266 9 
ee 
8 = 0,063 098 0 —0,041 197 6 
E © = (-000 2266 0,017 626 9 
0.777 X 0,000. 080 1 —0,000 593 9 


Die beiden in (19) auftretenden Matrizen x. 
Br 3 Werte: gi E 


0,012 430 4\ » /C, — 0,085 049 8 0 
—0,025 079 2) \C, En 0,000 8129 0 
Ü 


0,010 260 7/ \C, 0,000 0345 0 
0,008 650 3\ /C} 0,084 217 6 
0,016 6625] -|C,) + 0,002 754 5) 9, 
0,007 068 3/ \C, —0,000 052 7 


MI und S-TMIN haben! die folgenden 


N nn | 0,545 972 2 0,545 972 2-\ 
Br EREE ‚01 0,54 3752 ° —0,544 375 2 | 
BN. RR 0.1 0,046 4348 0,046 434 8 
REED &.m-ı —| —0,045 8460  -+0,045 846.0 |. 
PR sr 5 V.0,000.443 8 0,000 443 8 | 
a en DEN 0,001 4995  —0,001 499 5 
N 0.000 0188 0,000 018 8 
BN: - 0,000 028 7 


N ne wre 0,000 028,7 


0,303 7773 0.) 
"0... —0,065 2722 
0.322 3929 0 
0068 595 1 
00191319. DER 
70° 0,003 406 4 
A 00005193 BI 


—, 001 265 0 0 


0.055882 0 .0.0026077 0 


0 0,026 309 2 0 0,001 027 7 
—0,087 567 9 0 0,012 208 6 0 
0 —0,043 4133 0 0,008 563 8 


‚0,032 944 6 0  —0,025 0771 0 
‘0 0,017 699 4 0 —0,016 659 7 
0,010 260 8 0 


000000885 0 —0,0059535 0 0,007.0683 
angew. Math. Mech. 24 (1944), 8. 71/76: 
2 0 ‚222 565 \ 
0° 
a 0,149 977 
8 


Eingegangen: 12.4. 1950 


a en 
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Stellt man jetzt z. B. die Randwertaufgabe: r = Hank ( ). so wird aus (19): 


0,257 353 

0,544 375 | 
0.264 205 
y(x) = p’(®) | _0.045 846 
| 0,022 311 
0.001 500 
| 0.000 753 
0.000 029 


Die Differentialgleichung des Beispiels ist andererseils exakt integrierbar eh liefert als allge- 
meines Integral: 


ya), (& -&oj > 30, 04° Sin. er 4): 
Für die vorgelegte Randwertaufgabe wird aus diesem allgemeinen Integral die folgende spezielle 
Lösung: 


yla) = email (4 116 624 Co| - 5 + 1,232 045 Sin —:) 


Die folgende Tabelle liefert für einige «-Werte einen Vergleich dieser exakten Lösung mit der 
nach dem obigen Verfahren gewonnenen Näherungslösung: 


| 
2; | Yexakt Y Näherung 
—0,6 —0,06572 —0,06574 
—0,2 0,01366 0, ‚01368 
0,2 0,25803 0,25805 
0,6 0,62006 0,62004®) 


*) Kommt es nicht auf so hohe Genauigkeit an, so kann man die Reihe (4) früher abbrechen, Läßt man 
z. B. in (4) n nur bis 4 laufen, so wird jetzt ®, eine Matrix von nur dritter Ordnung. Hierdurch vereinfacht 
sich die numerische Rechnung ganz wesentlich, In dem obigen Beispiel würde man statt der vorstehenden 
Tabelle die folgende erhalten: 


C Yerakt | YNäherung 
0,8 | 0,066 | 0,064 
Ip:g | 0,014 4.0) 0,016 
0.2 | 0.258 | 0.257 
0,6 | 0,620 | 0,616 z 


\ 


Die Berechnung vollständiger elliptischer Ta erale dritter 


Gattung durch Reihen 
Von M. Kolscher in Darmstadt ’ 


Die Berechnung dieser Integrale ist an sich möglich durch Rückführung auf die söhelliertan Werte der 
unvollständigen elliptischen Integrale erster und zweiter Gattung. Jedoch ist hierbei der Genauigkeit eine 
Grenze gesetzt durch Dichte und Stellenzahl.der vorhandenen Tabellen. Außerdem wird das Verfahren un- 
brauchbar, wenn der Modul nahe bei 1 liegt. Für diesen Fall ist eine Reihe zur direkten Berechnung von 
G. Hamel angegeben, welche aber versagt, wenn auch der Parameter noch oe bei I liegt. Hier Ye 23 
wir Reihen, welche in jedem Fall eine Berechnung erlauben. 


The compulation of integrals of the kind mentioned in the title is formally Rack by BR : 
to the already tabulated values of the incomplete elliptic integrals of the first and the second kind. The accuracy 
is limited, however, by the intervals and the number of decimal places to be found in ‚the existing tables. B 
the method fails, if the modulus is near to I. In this case, G. H am el has stated @ series for the. im 
computation, which again fails, if Ihe parameter is near to 1, too. A tie are ask pern 


compulation in either case. > 


L’evaluation d’integrales elliptiques incompletes du troisiöme espdce se Doc, 
des intögrales du premier et du second espöce Mais l’ewactitude y est limilde par 
la „‚densite‘“ des tables ewistantes. De plus la methode ne fonetionne pas, a le mo 
ce cas une serie de Hamel permet une Evaluation directe, sinon le para st 
donnons des series, qui permettent l’övaluation en tout gas 5 
Bsrymenenue uHTerpauoB YKa3aHHOrO B ‚aaraanım muaR 
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mePBOrO H BTOPOro Pona. TOYHOCTB ITOrO METONA OAHAKO OTPAHMYCHA NONHOTON IM TOUHOCTBO 

HMEIIMHUXCH TAaOImm. KpoMe Toro HTOT CHOCOR HENPHTOACH B CRYYAE, KOTIAa MOAYAIB O1UBOR 

x eneuune. us aroro cayyası T. Tamerıp gar ps, NO3BOAMOLMa HENOCPEACTBEHHOE BEIUHC- 

 ACHHe HTHX UHTErPJIOB, HO KOTOPHIÄ TOME HENPHTONEH, ecim m apamerp ÖANSOR K ENeHuNE. 
= BTech AAMTCH PAABL, MOSBONNOMME BEIUHCHEHNE BO BCEX CIYYasıX. 


I. Anwendungsbereiche der Reihen 


R Als Normalform des vollständigen Integrals dritter Gattung werde benutzt 

® NET: 

8 | M(R, = rain oo RR 2 2 AR 
3 0<skt<1, 2>0. 

F Ist A? >1, so soll der Cauchysche Hauptwert verstanden sein. Wir geben eine Reihe, 


welche nach Potenzen von %? fortschreitet und für k? <1 konvergiert sowie eine nach Potenzen 
von K’?=1-—.%? fortschreitende Reihe, welche für X? <1 konvergiert. Praktisch sollen sie 
aber nur für kleines k? bzw. k’ 2 Anwendung finden, so daß wenige Glieder bereits hohe Genauig- 
keit ergeben. 

_ Liegt %® nicht in diesen günstigen Bereichen, so läßt sich durch Anwendung der Landen- 
schen Transformation ein kleinerer Modul oder durch Umkehrung der Transformation ein 
größerer Modul erzielen. Das Ergebnis der Transformation ist angegeben. Es wird vorgeschlagen, 


für k <0,8 (IId) und dann mit k; <ıe ! ; die Reihe (IIb) zu verwenden. Für % > 0,8 wird mit (IIe) 
kr < 31 und hierfür ist Reihe (ide) zu verwenden. Die angegebenen Reihenglieder liefern bei 


einmaliger Transformation mindestens 8 Dezimalen. 


II. Formelzusammenstellung 


Re a) Hilfsgrößen: 

E | = für  o>0 
RSEET ve 

L . en ERrER BE ah) 

Ale ER r=(0 für oes0. 

4 b) Reihe für kleines A®: 


5 5 DE: A,(nAN). 
Nn=0 


_ Darin ist für A2> k®: 


16 7 8 
BR 51 1 Re: a Is 
a Sa KL 
Ray? 1, h Rd ) 
alt 1 gg N 0 jag* 


Si +3 Ze Has? Fame, 


er +3 Ei + At net +0” Ach 


nl 


gr+ 


.—— 


PA AP EA GE 5 


y 
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1 
A=-— 8\s 
Are (5,8 RE it} A, 
H 16.116.198 256 


5 U D AR ) 

# pr: kun Se 

Me 65 

| | Alt): | 

e BA. c) Reihe für kleines u: "x 

: i a ee Ce : (a a N 

Ri M=M 4m $- Ze a e)— Di Eon (k,A>0). a 

En Nn=( = j 2 wa: 

? 2 Darin ist N: 

h Y ® ! h Deren 5 ; N i 

Be. ae EN 1 In ezFarz an für osask Te 
ER | MAN... ee la21 13 

Br. Kar 1—k $ en.) Z 

B | | En Ey — 1-arc BER Ras A für kzA< 1 a 

En N, und A, wie vorher mit Ersetzung von r, A® durch 1/A, e. Die erste Form der A, ist für |o| ZK'®, 


av die zweite Form für |o| <k'* anzuwenden. Die Berechnung der w, erfolgt für el Kk? aus den 
RN Hilfsgrößen. 


B=43:@-4=, 
FR 1 1 ae 3 
Bı= gar 7) Kent at 62), 
1 1 a, BETT 4 
Be— og (—8-+ 112 + 2742 — 10 4° — 20 4%) 
mittels \ \ RN 
1 2 
VOR: NT 
1 er | 
ey E DEE FRE u I 
Für jo|< A? ist 4 BL | EL eR ir ne = er 
SHE Up: n + tagt ren n.\ 
mi Ka, Te Aa Ss N = —_ RO he 
ee REN MT . | 


d) Transformation zur N 


1—K 


ns h 


range 


. 1—k- 
a. @TItE i 
Br Im Bereich K<A<1 ist nur +A verwendbar. Pr 


a iR 
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N f) Sonderfall: 


Ist A so klein, daß hierdurch in (c, d, e) ein zu großer Dezimalenverlust eintritt, so aus (1) 
durch Entwicklung nach A?: 


1m 12%: (244%) (IR) , 2-43 Kg M 
1 BE ERS TI kanal. 2000 GE ee | | Lu 
HA 3a Er 15 Ki rn: 


UI. Die Gewinnung der Reihe für kleines k? 
Es werde gesetzt 


1 1 
An ae) BEE Sa 
= (— 1): i — Ir - = a N ut 
RS. N El) n 1—2n (2). 
Aus (1) et durch Reihenentwicklung der Wurzel für (IIb): 
ER, . er. 2 hi sin?” @ 
ade r" Ru; su Ne 7 ELIA 


0 
Zur bequemen Berechnung der R, berechnen wir mit |ö/<1 
r/2 
_4 dp 
(1 — ösin?@) (1 — A?sin? 30° 


Es wird für A?<1 


RR DE OR a a MH N 
2 ö PLN 2, _9 Y1—ö 
IR u 
und für A4?>1 der Cauchysche Hauptwert 
% RE 1 
a R=— -— — 
i | e ae SI Feet Ze 
” A? 


Hieraus sind durch Boraerehehkuhwiekkung nach ö die R, und damit die bei (IIb) angegebenen 
Werte 4, zu finden. 


Es ist leicht zu sehen, daß diese Reihe für alle A? im Bereich A? <1 konvergiert. Ist aller- 


en dings 42 klein, so tritt bei Berechnung der A, ein mit n zunehmender Dezimalenverlust ein. 
E Dies stört nicht, wenn A? >k? ist. Ist aber A? <k?, so können wir auf die Darstellung der A, 


der kein Dezimalenverlust eintritt. 
Übrigens gilt für M die ee 


(A? 2). +: NE. SEN RR (3) 


La 


6 mit: K als vollständigen pluptischen Integral erster Gattung. Setzen wir hierfür die bekannte 
Reihe ein, so folgt 


E: | &: Ki 
| = n | Ih) (a 1. a2 Er) a 1= Mr Daran Reha (4) 


| wir stellen fest, daß unser M die einzige Lösung dieser bei A=A singulären Differential- 
King: te Ye auch dort noch regulär Babr: 


N A 2 JE IV. Die Gewinnung en: Reihe für kleines k' 


k=1 1 bleibt zwar M endlich, aber —— aBr wird en unendlich. Zur ne 


“a dk 


x 4, 
ze: 


in geschlossener Form verzichten und erhalten durch Reihenentwicklung nach A? die Form, in 


re 
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1 Eee 
ar 1 1—kx 
3. oe 6: 


Wir berechnen das Hilfsintegral 


‚| ’ 
Fe S n — 1+® Krk 
Dr 3 ra nren VE 


Dt 


Dies ist elementar möglich und liefert 
h N Ss Re Ss’ + Ss” +58” 


mit j 
VEN TE 1 1%, MAL -A ar Oak und a1 2 
| m We Te 
"—I+1 ; “ 
oder 2 
a! 1 k A—k 
S ei as DR YA ; aretg]/ ;y fü k<zA<1 
1+4 B 
sowie BELA R N: 
BU a [Er NUR - 
RR EDER DIR UN TIERE am SE 
er 14829  PRd4] ee Een 


5 29 RE RN 
Entsprechend wird a 
S_= S, En 8% + Dunn 


wobei diese Werte durch Beihenentwicklinn nach ö aus 5, 8”, 8” zu ı ggrinnen sind. Bine 


chend wird ferner 
M=M'+M'" Se m RE ES TREE (6). 
So folgt aus der ersten Form von $’ 
ee le „GESA + =) 
"U (1-+ A) 24 ur a7: 
und daraus nach © Ss | =, 
- 7223 BE n\1-+ u 24 ee 


2 Ba 
Br a ihn Eee 


von 8’ verfahren. e. j 
Dieselbe Rechnung ist auch mit 8” möglich, aber nicht mit 8”, das si in sehr u: plizier 
Form enthält. Aus beiden erhalten übrigens &,, Sn Nenner mit ne von 1—4. 


” Darin wird 


wird M’’ nur eine Reihe geben, welche für ® Be <ı konvergiert. AL 


Dies kann durch Zusammenfassung mit M’ BEnnEEn, werden, 
Zunächst setzen wir 4 re 


In —= 


an I 


E 
E 


Gemäß o (7) setzen wir ein. 


TE 
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0) i ee) 
j a 


ey 


a a, | BE EN SEAE 
re ee RR A VEREIN nen 


Die Ausdrücke in den geschweiften Klammern verschwinden für ö6&= 1-4. Daher sind P, Q 
auch hier regulär. Durch a: nach e wird 


Pe a Re “ Der, ma Oo 


worin p, 0 auch bei ö=1 Kr ne Funktionen von 6, A sind. Insbesondere wird 


1—71—Ö 


Re A IErt fer a 1+4 () sin ( ö\] 
la ar tra" 1-3)}: 
Werden nun die P,, Q, nach Potenzen von ö entwickelt, so ist daraus zu finden 
„ 11 N 4 < < e 
5 Sn =hr 22 e’ a! "on 
Insbesondere wird x 
> ;® 1 2 
Po=y: 90=5, WITT: 


Damit erhält auch M’’-+ M’” aus (5) eine entsprechende Form. Führen wir noche = 1 — Y1—h2 


als PO nrNe nach %’? ein, so wird endlich 


m" + M"'— = v,. ken —_ > TE En SER ERS RE (7) 
- B x n=) 
und a x 
\ MD. 6; ed RN EA SEHALAT ORRTR (8). 
Die Berechnung weiterer Glieder auf diesem Wege ist aber unbequem. Daher führen wir in (3) & 
als Variable ein. 'Es folgt." , 
- LH ‚x: 
f ee N ER ET (9). 


Mm wind AV ws 


we v, Ww die äh ma rn ae sind, nn M’ nach (IIc). 


‚ Mit ‚der für beide Formen von M ‚BUBEN, Relation 


I : aM : Ba .1+k, 1 1/k+R 


a 5 ER RE N ARE UNE N 5 ı. 
R Ra ae Da ( ) 


Ever dank. (an. 
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Die Gl. (10) geht aus (4) hervor, wenn darin k®, A durch K'2, o ersetzt werden einschließlich ? 
des entsprechenden Anfangswertes. Daher ist die Lösung bekannt. 


In (11) läßt sich die rechte Seite in eine Potenzreihe der Form k’- 3 B,- %'?* bringen 
n=0 
und liefert danach bei bekanntem Anfangswert w, die weiteren Koeffizienten w, von W. Ebenso 
wie bei den anderen Reihen ist auch hier X’? = |o| nur scheinbar Konvergenzgrenze. Für jo <K? 
ist zur Vermeidung des Dezimalenschwundes Reihenentwicklung nach o durchgeführt. _ 


V. Zusätze 
Die Legendre Normalform ist 


rl2 
7 dp I6R 42 
ITS.) | % = 7 El) + a MA, I). 
‚2 Re (1 — 42 sin?) YL— k2sin? 4 A 
Das Integral 
a N RN 
yl— k2sin?o ä 
TOD ; 0 
Me en B 
wird durch (IId) mit 4?=— u? in ein Integral der Form (1) mit <A} < 1 übergeführt. 3 
VI. Zahlenbeispiel 
ne k 1 K 
Zur Kontrolle wurde der Fall A—=0,9; k=0,8 durchgerechnet. Mit (Ild) wird k, = wi, 4 
= a Hiernach wurden von (IIb) die angegebenen sechs Glieder in der ersten Form PR 
1 LA Es 
berücksichtigt. Zum Vergleich wurde auch nach (Ile)gerechnet, wofür k, = Y = E wird. A 
1 ; 
B Danach wurde mit allen Gliedern von (IIc) in der zweiten Form, wegen Orr 1° k’r, ge 
E rechnet. Beide Rechnungen gaben übereinstimmend 
# M (0,81; 0,8) = 2,6400 284. | 
B; Bei anderen Werten von % ist die gleiche Genauigkeit mit weniger Gliedern erreichbar, so daß 3 
der ungünstigste Wert gewählt wurde. Der Wert für A wurde so gewählt, daß die Hamelsche 7 
Reihe nicht mehr konvergiert. | di. 
“a Eingegangen: 24. 4. 1950. s E 
N. s x g ‘ E 
> Analytische Untersuchungen und konstruktive Erweiterungen % 
= - 
hr zu den graphischen Integrationsverfahren von M ei ß ner R 
E Von Rudolf RR in Braunschweig) Br 2 
N ” Für die graphischen Integrationsverfahren von Grammelund Meißner existierten bisher weder Er 
; eine analytische Behandlung der Konstruktionen noch eine daraus sich ergebende Fehlerbetrachtung. De 


analytischen Rechnungen führen für die@Grammelschen Polaren zu den Verfahren mit den kleinsten 
Fehlern und für die Meißnerschen Linienbilder zu neuen Konstruktionen, die in einer eigenarligen 
Analogie zu den Gr am melschen Methoden stehen und die nur mit einem Fehler in der EP TETE 
des Fehlers bei der Simpsonschen Regel behaftet sind. ET, 


For the graphical methods of integration due toGrammeland Meissn er, ale an "analykidal 
treatment nor a computalion of the error based on the former, were existing till now. The computation leads { 
to conclusions that keep, in the case of ram m el’s polars, the errors as small as possible, and give, in Na 2 k 
case of Meissner's line diagrams, new constructions. T’hese are, in a particular manner, a „2 K” 
Gra m mel’s methods, and are affected only with an error, the order of zoNEh is the aaa as that of 2 i MS Rs N 
sons’ rule. x 


Pour la methode graphique d’integration de Grammel et de Mei is 8 ner il 
ni un traitement analytique des constructions ni une methode pour l’ Evaluation 
ee traitement. Par les operations analytiques on obtient le proced& au moindre er 
deGrammel et, dans le cas du diagramme de Meissner, des construct 
remarquablement analogues aux methodes de M eissne 1 Perreur est 28 meı 
de la rögle de Simpson. { ! 


1) "Diese Arbeit und zwei demnächst erachbinende kibalen u 
Arbeiten am Tin b Ibn t FRI Bene an 3 1 SM Dt LS 
derTechnischenHochsc eeTBr Var Ken 


SAY 
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Aus rpabmueckux MeToXoB HHTerpmpoBanns T'pammena u Meüchepa no cux op eme 
HeT HH AHAJIHTHYeCKOrO MHCCIeNOBAHUA THX KOHCTPYKUNÜ Hm OMeHOK NOrpemmocrei, 
BEITEKAIIINX 3 TAKOTO NCCHEeNOBaHuA. AHANNTHYecKHe PacyeTs IPUBONAT B CAyyae Nonsp 
Tpammena K METONaM C HAHMEHBIIeH BOSMOSKHOH NOTPeEMHOCTBR, a B ca yyae JINHeÜHBIX 
auarpamm Meiicnepa K HOBBIM KOHCTPYKIUUSM, KOTOPBIe AHANOTHYHEI METONaM T'pamneıs 
H KOTOPBIM COOTBETCTBYET UOTPeNIHOCTB TOTO-3Ke NOPAAKaA, KaK U B c/Iyyae TITPHMeHeHnHsI 
IpaBna ÜuMmıcona. 
Einleitung 
In den folgenden Ausführungen wird von den Arbeiten. von Meißner [1] und Grammel 
[2], [3] über graphische Integrationsmethoden ausgegangen. Durch analytische Behandlung 
der Grammelschen Konstruktionen ist es erst möglich, die günstigste Art dieser Konstruktionen 
herzuleiten, bei denen man den Fehler abschätzen und vor allem möglichst klein machen kann. 
Die Vorteile dieser Konstruktionen werden dann in geeigneter Weise auf die Meißnerschen 
Linienbilder einer Funktion übertragen. Aus der analytischen Formulierung der Konstruktions- 
verfahren gehen auch praktisch verwendbare, numerische Integrationsformeln hervor. Während 
hier also Methoden für bestimmte Integrale behandelt werden, soll es einer späteren Arbeit 
vorbehalten bleiben, die Verfahren auf Differentialgleichungen auszudehnen. 


Das analytische Verfahren 


_ Zum Zwecke der Fehleruntersuchung wird den geometrischen Konstruktionen die analy- 
tische Rechnung gegenübergestellt. Dabei wird ein Vektor P der Länge y und der Amplitude & 


"in der Gaußschen Zahlenebene durch P=yei* dargestellt. Je nach Zweckmäßigkeit wird im 


folgenden sowohl von dem Vek- 
tor P als auch von dem Punkt P 
als Endpunkt des Vektors ge- 
sprochen. 

Für zwei Vektoren pP m yelıı 
und P,=y,ei% wird das skalare 
Produkt eingeführt: 

Pi‘ P3= yıy3 cos (&ı — 9). 
Für den zum Vektor P, gehörenden 
Einheitsvektor ep, gilt also z.B: 

Li * ep, — Yı C0S (2 —— %). Bild 1. Grammelsche Polaren 


. 1. Die Grammelschen Polaren 
1,1. Die Darstellungen einer Funktion und ihrer Ableitungen (Bild 1) 


Als Bild der Funktion y=f(x) wird die Grammelsche Polare (C) definiert durch die 
Kurve, die der Endpunkt des von & je Vektors 


P(«) eit Kr ER re. 


je 
beschreibt. In ähnlicher Weise werden die Bilder der Ableitungen von f(x), die Orthopolaren 
(Cr), definiert: 

el) zu er De TE ER TE RR 9) 


= 2 2 


Die Haupteigenschaft der Polaren und Orthopolaren läßt sich nach ranımel folgender- 
maßen aussprechen: 


Satz: 


Die Tangente in einem Punkt der Polaren (C)) [bzw. einer Orthopolaren (C7;)] geht 
durch den zugehörigen Punkt der Orthopolaren (C,) [bzw. der darauf folgenden 
Orthopolaren (O;+1)]- 


„Zugehörige Punkte‘ sollen stets solche genannt werden, die verschiedenen Polaren bei 


‚gleichem Argument z angehören. 


Beweis: Einen Tangentenvektor {(&) in einem Punkt P(x) der Polaren (C) erhält man durch 
..Differentiation des Vektors P(x) nach «: 


LEN EL ir 
ER ne Fa Ha) FO 
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Daher ist 

P(x)-+ At(x)= P’(x). 
mit A= f(z)/f’(z), womit der Satz bewiesen ist, 


1,2. Die Integration nach dem Sehnenverfahren (Bild 2) 
Bei der Grammelschen Darstellungsweise der Funktion und ihrer Ableitungen besteht 
das Integrationsproblem darin, aus der bekannten Orthopolaren die Polare zu konstruieren. Die 
gegebene Orthopolare ist also für die Funktion f(x) zu zeichnen: 
{ 3 i(®% te) 

ie 


P (a mit den Punkten Pı= = , 


Fi 


R“%=Mtke, 
f=f(o+ te), t<0<1. 


eist dienoch näher zu erläuternde Schrittweite 
% +8 


des Verfahrens. Um das Integral Ri f(a)da 
1 


wobei 


“m; 
rau: zu bestimmen, soll aus demVek- 
e'(®k + ‘) 


tor Pr ER "der neue Vektor PL .ı = Par 

k k+1 
dadurch bestimmt werden, daß man auf der 
Sehne des Kurvenzuges fortschreitet. Aus 
Bild 2 folgt sofort 


PR=(l —/4,) Pt Dis 


Um A, zu bestimmen, wird diese Vektorglei- 
chung skalar mit ep, multipliziert, also ist: 


5 % 
Re} BR BTUENE h 
ER, "ep, fr sin e i 
Daraus folgt nach kurzer Rechnung: 
1— “; 
Bild 3. Grammelsches Sehnenverfahren Ir 1 = Ik = [ te ” af a 
t F 
Man sieht, daß nur für {—=1/2 der Näherungswert des Integrals vom Anfangswert auf dem L 
Vektor Do Gep,, d.h. von o unabhängig ist. Grammel spricht z. B. nur vom „Bogenmitte- 
punkt‘ eines Shah aber nicht exakt von der Winkelhalbierenden. Für t=1/2 gilt y 
4 = 96 + 2fesin , re: REN" ( 
” Reiht man n Schritte in gleicher Weise aneinander, so ergibt sich R 


5 4 N, 
n=9+2[lfo tft" tA-ılsing 
Es ist also f rer, 
.. © c > 
N Be #E a: ß 
‚rayde > n—-9%=2[lfho tft :'' +A-ı]sin 2 IR ie PN 
Man kann das Anelakeen noch so modifizieren, daß man die Orthopolare ı 3 kai we man ar % 


aufträgt Ya e 
ver (®o 27 b Äh ie II n ar ; Ä [ER Si A 
um. nr Be 5 


So erhält man als Integrationsformel die bekannte Maclaurinsche 


[1ada- BFH ft + Fe Be id, a 


13. Die Integration nach dem Tangente %, 


Wie in (1,2) sind die Punkte P; und P; ı ı der ‚Orthiope 
bekannt. Man findet 5 +1, indem man auf der Tangent © 


NY 
3%, 


Pı= 


a U un a A En 


Nr. rs Aprit/taı 1951 Ludwig, Zu den Integrationsverfahren von Meißner und Grammel 123 


er |) — 


rungspunkt“ von PP; ,, wie Grammel sagt, d.h. exakter bis zur Winkelhalbierenden [wie 
eine Rechnung ähnlich der von (1,2) zeigt] und dann in der Tangentenrichtung des gesuchten 
Punktes diesen erreicht. 


1. Schritt: 
Pi=AyPr+(1— A) Pr, 


” i (2. + 5) 


a ERBEN ZEN 
cos —- 5 + Ik sin. 
2, Schritt: 

Pr+ı =yPrrı +(l —uı) Pı, 


e' (ap-t+ &) 
P, 4= 


| RE, EREFEE ä 


Also liefert die Integration über ein Inter- 
vall: s 


Kr — m (fr se 


und über » Intervalle: 


{ Bild 8, Grammelsches Tangenten vorfahren 


7 %y £ 
a io) da ne ger (Jo “r afı ar A Ep 2 1 +/n) ig Tr A (7). 
Durch eine ähnliche Modifizierung der Polaren wie in (1,2), indem man für die Schrittweite 


tg 5= a setzt, erhält man die bekannte Trapezregel: 


2 
3% h Br 
jfode= 5 (kraft +2 — 1 ++ 2 El a RR (8). 
17} e 

* ) u r 
\ 2. Die Meißnerschen Linienbilder 
Er Die a AR einer Funktion und ihrer Ableitungen durch Polaren und 
Linienbilder (Bild 4) f 


Die Me Onsrachen Polaren unterscheiden sich von den Grammelschen dadurch, daß 

die Länge des Vektors gleich dem Funktionswert selbst und nicht gleich dem reziproken Wert 

desselben ist. Dies vereinfacht das Verfahren. Die BeSneTuehen Polaren werden also durch I: 
folgende Vektoren definiert: 


Bi TE P(x)—=f(a)ei® een a BERN URN (9) 
Ei; Poa)=ifü@er (| 
- Darüber Hinans werden aber von Meißner noch die Linienbilder eingeführt, deren Punkte 


J a man durch vektorielle Addition zugehöriger Vektoren (d.h, Vektoren für dasselbe Argument ®) 
aufeinander folgender Polaren erhält, also: 


gi 5 ® Ola) = Pa) + P’(a) (2), \ ned 
gl la)= Pirla) + Patdle) (nd) 


peienscaten der Polaren und Linienbilder kann man in den folgenden drei 
nenfassen; i 


nale in einem Punkt des Polarbildes (©) [bzw. 01] geht durch den zugehörigen 
‚Polarbildes (C,) [bzw. Or, 1] der ersten [bzw. (%-H-1)-ten] Ableitung, 


x: ee) Hiraler = — P(x)+ P'(e). 
ei: 2 über ı die ‚ Polarsubnormale ı) 


Fan ae l va Bi 
-. we” A) 
3er‘ re TR Fi 
BAT a PL Bir, 
£ ar Bl. 8 5 
MEN SET 
ı ar 
ER 0 RE N + 


ee" BA Pla, BEE EEE I SEE so 
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Satz3a:| Die Tangente in einem Punkt Q(«) des Linienbildes Z (bzw. L;) ist das Lot im End- 
punkt des Vektors P(x) der Polaren € (bzw. O7;) auf diesem Vektor. 


Satz3b:| Das Linienbild Z ist die Hüllkurve aller Lote in den Endpunkten der Vektoren Pr 
der Polaren C (entsprechend für .L; und (7) auf diesem Vektor. 


Bun Tangente des Linienbildes: 7 = a | Fr | | 


Normale des Linienbildes: N=:<T=— [ f(x) +f” (x)] = ber: 


Bild4. Meißnersche Polaren und Linienbilder 


Satz 4:| Der Krümmungsmittelpunkt für einen Punkt des Linienbildes Z(bzw. Z;) ist der zuge- 
hörigePunkt des darauf folgenden ersten ‚rsten[bzw. (k-+1)-ten] Linienbildes Z, (bzw. L; + ı). 


Beweis: Aus Q(x)=[f(&) + if (JE dr in folgt nach der bekannten Formel für den 
Krümmungsradius ebener Kurven in Parameterdarstellung: 


BE 
in. 


Einige Hinweise sollen noch gegeben werden über Besonderheiten im ı Kurvenverlauf der 
Polaren und Linienbilder, deren Kenntnis für Konstruktionen vorteilhaft ist. 


Über die Schnittpunkte von Polaren und Linienbildern in zugehörigen Punkten güt: 


TER ERRRR 


Satz 5:| Polare und Linienbild haben dort gemeinsame Punkte und Tangenten, wo die erste A 
Ableitung der Funktion verschwindet. 


Beweis: Schnittpunkt von Polarbild C und Linienbild ep 


P(«)=Q(2)= P(x) + P‘(e), A, Es u | 
also PXr2)=0 oder F(@= | 


Geometrisch sieht man leicht ein, daß beide ASHEVeN an ‚der Beirdeteten Stelle 
eine gemeinsame Tangente haben. 


Satz 6:| Das Polarbild der ersten Ableitung’ und das Linienbild einer Funktion San 3 
sich in Punkten, für die die Funktion verschwindet. Fe AFTER 


Beweis: Schnittpunkt der 1. Polaren €, mit dem Linienbild L: 
P’(«)=Q(2)= P(2)-+ P'(«), 5 
also P(x)=0 oder Sal) EN 


Aus einer Figur folgt sofort für den Schnittwinkel der Da Kurven 
Tr 
t jo I 
A 


A a a ee Ba A TE 
u er U Mr N i a 


- 
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Eine Singularität des Linienbildes behandelt der 
Satz 7:| Wenn f(x) + f’’(x)= 0 ist, so hat das Linienbild der Funktion im allgemeinen an 
| dieser Stelle eine Spitze. 
Beweis: Der Tangentenvektor 7(x) des Linienbildes (x) 


I zit Hr 


kann verschwinden, d.h. kinematisch betrachtet, die Bahngeschwindigkeit wird Null. 


NE 


2,2. Die Integration nach dem Sehnenverfahren für Linienbilder (Bild 5) 


Bei der Darstellung einer Funktion durch die Meißnerschen Polaren und Linienbilder 
besteht das Integrationsproblem darin, aus der 1. Polaren C, unter Verwendung des Linien- 
bildes Z-die Polare C zu konstruieren. Hierbei wird im Gegensatz zu Meißner von dem Satz 3 
über die Krümmung der Linienbilder kein Gebrauch gemacht. Solche Konstruktionen lassen 
sich weniger leicht analytisch bis zu praktisch brauchbaren Fehlerabschätzungen überblicken 
und gestatten nicht ohne weiteres die später beabsichtigten Anwendungen auf Differential- 
gleichungen 1. Ordnung. Die hier angewandten Konstruktionen weisen eine eigentümliche 
Analogie zu denen von Grammel auf. Sie sind ebenso einfach, erlauben aber die direkte Ab- 
lesung der Funktionswerte ohne Berechnung reziproker Werte. 


Mit den Punkten der Polaren: 


Pi = ifre' **, Prrı=ifrie Rt, P= ge” 
ist gleichzeitig der Punkt Q&,— Pı-+-P; des Linienbildes gegeben. Den Punkt Q; . ı erreicht 
A' 


man auf der Sehne, die parallel zum mittleren 
Vektor P; ist. Also 


Qrrı = Pı+ı E Q1epr 41 = a7 -r 1, > 
wobei o, und o, zu bestimmende Skalare sind 


Da E7A a “Lr 
und e_,=te | Si >) ist. Durch skalare Multi- 


k 
plikation obiger Vektorgleichung mit N fin- 
det man: 


1—=g+ Bi ei) 5 


| und daraus 

Be 4 | tat) 

. RAmHm=latftfr)tsz|e 
Es ergibt sich also wieder die Gl. (7): 


9er — Ir (fr ERS) tg 2 


Also gilt der Satz 8: Bild 5. Sehnenverfahren für Linienbilder * 
5 
Die Integrationsmethoden nach dem Sehnenverfahren für Linienbilder und nach 
° dem Grammelschen Tangentenverfahren liefern die Bleiahen Benee, näm- 
SR RR lich die Aare eel)). 


2,3. Die Integration nach den Tangentenverfahren für Linienbilder (Bild 6) 


s. E „ Entsprechend wie in Abschn. 1,3 erfolgt die Integration über ein Intervall in zwei Schritten, 
“indem man in der Tangentenrichtung des Linienbildes bis zur ‚Mitte‘ geht und dann in der 
%  Tangentenrichtung des gesuchten Punktes diesen erreicht. Gegeben sind also wieder Pr: Pu 


= Er Di; Bi ı und Hapit auch @r. 
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1% Schritt: 
Q=Pı+ Q2ep: = P; % ee 


k 
liefert | 
si = 
92T cos m 2} ge sin 35 
2: Schritt: 
Pr+ 1 =: O3 @pz41 5 %% + O3 Pur? 
ergibt x 
ET 
= 2fsinz REITEN 4 


Man erhält also Gl. (4) aus Abschn. 1,2: 
he 
+1 — n—2fr SIn:n 


Bild 6. Tangentenverfahren für Linienbilder Also gilt der Satz 9: 


Die Integrationsmethoden nach dem Tangentenverfahren für Linienbilder und 
nach dem Grammelschen Sehnenverfahren liefern die gleichen Ergebnisse, 
nämlich die Maclaurinsche Regel?). 


2,4. Restabschätzungen für das Sehnen- und Tangentenverfahren 
Das kombinierte Sehnen-Tangentenverfahren für Linienbilder \ 
Für das Sehnen- und Tangentenverfahren mit der Festlegung der Schrittweite durch x 


sin 2 -5 bzw. 185 -, ergeben sich die bekannten Restglieder für die Trapez- bzw. 
Maclaurinsche Regel: so 
5 
rm = Sf" (+ ih) A | 


13 va 19% 
kan = Dr (+ Tah) : 3. 1,2 


So wie man die Simpsonsche Regel als Mittelwert aus der Trapez- und Maclaurinschen 
Regel herleiten kann, so sollen auch hier beide Verfahren gleichzeitig angewandt und ein Mitte- 
wert so bestimmt werden, daß der Fehler möglichst klein wird. e 1 


P2 - -. 


Für die Schrittweite 2A gilt: ‚ 2% ; \ 
ver = PO + lin +/)tsh, Bar re 


FO 2 — gi rw 
dabei gilt der obere Index ($) für Punkte nach dem Sehnenverfahren, (7) für Punkte nach dem 


- Tangentenverfahren und (M) für Punkte, die gemittelt worden. sind. Aue dem ar bei a 
der Simpsonschen Regel bildet man den Mittelwert: a Br w. 
1.7 | rm ER BEP 

Ir—z (gete + 2 064) Be Er i Kr Br 

oder: L ED N ee 


2; t2h R 4 i uhr, hr, 


ie Fa) da = ver u „20 Rs R= 3 +) tg "+ Ahnen... 


Zur Bieten des Restes R werden die einzelnen Ausdrücke in n (ar i 
nach Potenzen von h entwickelt. Man N ei e ER SANS 


B= Jan — I= wL” 


iRsiipre 
!) In der modifizierten Form nach 1. nr - N 
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Der Ausdruck in der Klammer ist dabei an einer Stelle & innerhalb des Integrationsintervalles Pr; & 


“ 


zu wählen. Das Ergebnis läßt sich zusammenfassen in den 
Satz 10: 


Bei der graphischen Integration über ein Intervall 2 A von x. bis &; + findet man 
für die Linienbilder durch das Sehnenverfahren den Punkt Pit, durch das Tan- 
gentenverfahren den Punkt Pi ». Teilt man die Strecke erh piD: innerlich 


im Verhältnis 2:1, erhält man den neuen Näherungswert Pi%,, der ein Analogon 
zur Integration nach der Simpsonschen Regel darstellt mit dem Restglied 


MS TOOL H DEN LEINE 


Bei Anwendung der Grammelschen Verfahren muß man allerdings wegen der reziproken 
Funktionswerte die Mittelwerte numerisch bilden. 


Beispiel (Bild 7): Das Sehnen-Tangenten-Verfahren soll für die Funktion 2 


A > % i 
angewandt werden und zwar für einen Bereich O<xz<z, also 
Bi 
Sie- et, Sim — Be: 
der sog. es Eine leichte nee ergibt: 
st | sin | ta. | 
0) v } 


Gewählt wurde als Schrittweite e=0,1. Die Wertetabelle kann ohne weiteres berechnet und 


3 <0,8 die Einheit e, — 50 mm 


die Polare gezeichnet werden. Im Original wurden a) für O< 
und b) für 0,3< 3= 1,2 e;= 200 mm gewählt. Für b) wurde die Polare um 0,8 gedreht und 


der Anfangspunkt der Integration beliebig gewählt (Verkleinerung im Bild 7 etwa 3/,). Man 
liest aus Bild 7 ab: 
817 2,4== 1,752, DET ID ODE 
Diese Werte stimmen mit den exakten Werten [4] 
$9?:2,4=1,7525 und Sir 1,8519 


gut überein (Sin durch quadratische Interpolation aus den benachbarten Werten). Man sicht KR 
an dem Beispiel auch den Satz 5 Den, 


2,5. Das Sehnen-Tangenten-Verfahren für Tipieahiuerahersin fach unterteiltes Sn t 
Intervall (Bild 8) | AuchE 
Ein Integral engl 3 
+kh { an L- “ re 


T Sean Wu 


soll näherungsweise graphisch ae werden in den folgenden zwei Be e 
1. Schritt: %-seitiges Sehnenpolygon: BR 
Nach Gl. (7) gilt 


®—9—(fo+2h+:+2f-ıt tg z ER 


2. Schritt: Polygonzug; bestehend aus dem Tangentenstück i im Ausgangs 
Sehnenpolygon und Tangentenstück im Andpunk ea 3 


Nach den Gl. (7) aa (12) gilt: B‘ 
= Irre, | 


ae ee ap 
gr = get cos h Fr sin BEN | 
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E) 


Daraus folgt: 
Er R [ h 
Rp lkhrh)sind (hi r2hr N Hafer )tgncoh, (17). 
Aus den beiden Näherungswerten für das Integral 
h r 
FL = (fo+ 2fı 4 or ar — fe ı+ fr) ig Do ee ne Nas 15 ee een (16*), 
3 / 
KBelfih na k(nteht. 2A tfı)tg 5 008 h. .(17*). 
soll nun durch geeignete Mittelwertbildung 
NEN a. (18). 


eine Annäherung möglichst hoher Ordnung in A erreicht werden. Die Taylorschen Reihen für 
den exakten Wert J und für J, und J, stimmen in den ersten beiden Gliedern, d.h. bis h?, überein 


Bild 8. Sehnen-Tangenten-Verfahren für %-fach unterteiltes Intervall (k = 4) 


Soll in (18) auch der Koeffizient von A? für J mit dem von J übereinstimmen, so muß man 
wählen 


Mit diesem Wert t hat man aber zugleich erreicht, daß auch der Koeffizient von h* aus J mit 
dem von J übereinstimmt, und das Restglied lautet: 


R% 2 ER oe a NW 3 MR (20). 


Satz 11: | Bildet man aus dem Näherungswert J, des Sehnenpolygons und dem Näherun 


polygon J, des Sehnen-Tangentenpolygons den Mittelwert 


J=th+(l—-0)5 


85" 


mit 


.,5k—6 

oR=i) 

| so ist das Restglied von 5. Ordnung in A [Gl. (20)]. 

Tabelle der Koeffizienten tund 1-18 
: 2. j Tabellel 


t 


130 Ludwig, Zu den Integrationsverfahren von Meißner und Grammel Ba. ne ; en. Mn 


2,6. Formeln für ein neristheä Integrationsverfahren 
Ersetzt man in (16*) und (17*) die EIBOUGME EL Funktionen durch das erste Glied 
ihrer Taylorreihen, also 


tg u sinh—h, cos h—1, | . 


so findet man folgende Näherungsformeln 
=; (ht2ht+t +2)» | 


IP hfte + (ht 24 + 2feat ho) 
Entsprechend Gl. (18) wird der Mittelwert 


Pi tAi—)g BR EEE N . (22) 


gebildet. Durch eine Rechnung, die der des Abschn. 2,5 erht ähnlich ist, findet man denselben 
Wert für t [Gl]. (19)]: 


. (21). 


5k—6 
6(k— 1) 


und es lautet die allgemeine Näherungsformel 


7 es OR EHNHEH] ,....n 


+12 (k-1)(k+ +9} 


BEE 


f N & j' 
aa a De a En er u Se an 


4 
mit dem Restglied 

\ k(ö5k— 6 » 

ee > nn Wr . (24). 


Die nachfolgende Tabelle 2 enthält die vollständigen Formeln für 1ISks6: 


h | 4 
akt fı) Sn 8 h 0 


h | ROLE 
alh+r2fitf) 2hfı 3(lhtifhtf). 


h 13% 7 ee BR | 
s(hr2ht2k+f) h(fı+ fe) as hr 3h+3hrh) n u Hl i 


h . h i 
jaWwt2htent2n Aapt2n Amon sn+ + 


h Erg am : ” | were 
lt 2 t+2h+2h+S) |FCh+2n+3H) ist AERO TEN f 
+2htß) +24 3) + 53419 [A N a 


h 
glhr2fitrhrkrf n BArzhpth 2 + If HI + + Ne 
+ fs) + fel En + Be Le LE ze | 


Newton benannte Formel 15]. 
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Berechnung der turbulenten Strömung im Rohr auf der 
3 Grundlage der Mischungsweghypothese 
Von W. Szablewski in Berlin 


Die Differentialgleichung der turbulenten Rohrströmung mit dem Prandtischen Schubspannungs- 
ansatz wird von der Wand her in zwei Schritten integriert. Für den übrigen Bereich der Rohrströmung wird 
die erhaltene Lösung als Approximation benutzt. Durch passende Wahl der empirischen Koeffizienten der 
Theorie erhält man gute Übereinstimmung mit dem Experiment. 


The differential equation of the turbulent pipe flow is integrated, by use of Prandtl’s relation for the 
shearing siresses, in Iwo steps starting from the wall of the pipe. The solution obtained is taken as an approxi- 
mation for the remaining part of the pipe flow. Good conformity with the experimental dates is reached by 
appropriate values of Ihe empirical coefficients occuring in the theory. 


A V’aide de la relation de Prandil pour l’effert de cisaillement l’&quation differentielle de l’Ecoulement 
turbulant dans un tuyau est integree en deus pas partant de la paroi. Pour le reste de l’&coulement dans le 
tuyau la, solution obtenue est employee comme approwimation. Si les coefficients de la th&orie sont con- 
venablement choisis, on oblient une concordance salisfaisante avec V’ewperience. 


Äußepennnansnoe ypaBHeHHe BUXPeBOrO IHOTOKA B TPyÖax HHTerpnpyerca OT CTEHKH 
Tpyö, Ipu yeMm AI BUXPeBEIX CABATAMIINX HANpsKeHnä TpuHnmaerca hopmyıa Ilpauırıa. 
llonyyenHoe pemeHue uCNONBSyeTcH B KAyecTBe HPHÖNN:KeHHOH POPMYJIEI 1A OONACTHU LOTORA, 
Ha KOTOPYyE He PACHPOCTPaHANIOCh UHTETPHPOoBaAHne ypasHeuus. Ilyrem nmonxonamero HoN0Op& 
OMUHPHYECKUX KOAPHIHEHTOB B Teopernyeckof hopMyıle MOCTHTAeTCH YAOBIETBOPHTENBHOE 
COBIAAeHHe C AAHHBIMH OIIEITA. 


I. Bezeichnungen 


u,v Geschwindigkeitskomponenten. t Schubspannung. 
P 0 Dichte. to Wandschubspannung. 
u Zähigkeit. ip. , 
vi n Kinematische Zähigkeit. N, 
p Druck. ZuE, Widerstandsziffern. 
x, r Zylinderkoordinaten. - dp 1 
fo Rohrhalbmesser. iss a FE 
d Rohrdurchmesser. iz 


) 
y=r,—r Wandabstand. | Y, 
en: V 70 Sog. Schubspannungsgeschwin- 


n=1 —; Dimensionsloser Wandabstand. 


| f digkeit. 

r - 1 Prandtlscher Mischungsweg. Toll: } Zitfer d 

h; x, u Beinen 5 De AT, u, nad = 

& nung Pac ro an d) o AR +9 DimensionsloserWandabstand. 
£ «u, Geschwindigkeit in Rohrachse. BE 

, w Mittlere Durchflußgeschwin- 6, Dicke der laminaren Unter- 
=) } digkeit. schicht. 

je R= U,d 9% 00 do : F > } 
E AT. » 1 Dimensionslose Dicke der lami- 
4 dj Reynoldssche Zahlen. j naren Unterschicht. 

3 IR Fi No 


II. Einleitung 
“ Im stationären Fall stellt die turbulente Rohrströmung (zusammen mit dem zweidimensio- 
nalen Analogon der parallelen Kanalströmung) im voll ausgebildeten Zustand wegen des Fort- 
“ falls von Trägheitswirkungen das einfachste Modell einer turbulenten Strömung dar. Die Diffe- 
_ rentialgleichung der Bewegung lautet in Zylinderkoordinaten, wenn wir den rotationssymmetrl- 
chen Charakter der Strömung beachten, 


e Ban: Eau; a 
; ET) a EHRT (1) 


vr 


9* 


ware Ba in EN =’ y0, 
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j r 
PER 


Zu der laminaren Schubspannung 
- du 
Tum =H- % 
(y=r,—r Abstand von der Wand) 
tritt noch eine scheinbare Schubspannung, die durch den turbulenten Mischungsvorgang erzeugt 
wird. Wir legen unserer Berechnung den bekannten Prandtlschen Ansatz [1] 


du 
ER = en(ze) a 


für die turbulente Schubspannung zugrunde. (Das Überstreichen deutet dabei an, daß es sich 
um die zeitlichen Mittelwerte einer stationären turbulenten Strömung handelt. Im folgenden 
werden wir die Querstriche fortlassen.) 

Hier bedeutet ! den sog. Prandtlschen Mischungsweg, der die Vorstellung fixiert, daß der 
Vorgang der turbulenten Vermischung durch sog. Turbulenzballen mit einemmittleren Mischungs- 
weg erzeugt wird. Die Verteilung des Mischungsweges Z über der Grenzschichtdicke, also die 
funktionale Abhängigkeit ! von %, bildet das zentrale Problem der Prandtlschen Theorie. 
Bisher existieren nur einige empirische Aussagen, die die Analyse experimentell ermittelter Ge- 
schwindigkeitsverteilungen an Hand des Schubspannungsansatzes ergibt. Hier sind vor allem 
die Messungen von Nikuradse [2] zu nennen. Es zeigt sich dabei erwartungsgemäß eine Ab- 
hängigkeit des Mischungsweges von der Reynoldsschen Zahl des Strömungsvorganges, die 
allerdings für größere Re-Zahlen sehr gering wird. 

Für die unmittelbare Wandnähe hat Prandtl [3] aus dimensionsanalytischen Gründen den 


Ansatz 
Imny NN WR 


gebildet und aus ihm für einen gewissen wandnahen Bereich, der charakterisiert ist durch die 
Annahme einer konstanten Schubspannung und die Vernachlässigung der laminaren Reibung, 
unter Zuhilfenahme dimensionsanalytischer Überlegungen das bekannte logarithmische Ver- 
teilungsgesetz der Geschwindigkeit für turbulente Rohrströmungen erschlossen!). Diese Formel 
gilt allerdings, ihrer Herleitung gemäß, nicht bis zur Wand, wo sie eine logarithmische Singulari- 
tät besitzt. Ferner treten nach Rohrmitte zu gewisse Abweichungen auf, die für die Darstellung 
eine Modifikation der im logarithmischen Verteilungsgesetz auftretenden empirischen Koeffi- 
zienten erfordern. 

In der folgenden Untersuchung werden wir eine Integration der Differentialgleichung (1) 
„von der Wand her‘ vornehmen. Genauer gesprochen, werden wir die Integration über die lami- 
nare Unterschicht, das sich anschließende Gebiet gleicher Größenordnung von laminarer und tur- 
bulenter Reibung bis in das Gebiet überwiegend turbulenter Reibung hinein erstrecken. Dabei 
werden wir uns für das letztere Gebiet auf einen Bereich beschränken, der definiert ist durch 


un tuld) EN RE EN 


0 9 
(=> wenig verschieden von), 

für En also die Mischungswegverteilung noch durch die beiden ersten Glieder der Taylor-Ent- 
wicklung beschrieben wird. Insofern werden wir uns also auf mehr oder minder größere Wand- 
nähe bei der Integration beschränken. Es dürfte jedoch keine Schwierigkeiten bereiten, unter j 
Inkaufnahme eines größeren Rechenaufwandes weitere Glieder der Taylor-Entwicklung für //r, zu ii 
berücksichtigen und so zu größeren Wandabständen vorzudringen. Die derart gewonnenen Funk- j 

tionen, die innerhalb des definierten Bereichs die strenge Lösung darstellen, ergeben nach den. 
bisherigen Untersuchungen in dem an den Gültigkeitsbereich anschließenden Bereich, eine gute 
Approximation der Geschwindigkeitsverteilung. 
Zur weiteren Berechnung der Strömungsvorgänge gehen wir dann so vor, daß wir die be- eh 
rechneten Funktionen, die natürlich die Randbedingungen an der Wand erfüllen, der äußeren 
Randbedingung u = u, (u, Geschwindigkeit in Rohrachse) für y= r, unterwerfen. - N 
Die empirischen Koeffizienten, die im analytischen Apparat auftreten, sind die Koeffi- 
zienten der Mischungswegverteilung entspr. unserem Ansatz: -(4). Des weiteren tritt als Smpige 
scher Koeffizient in den Rechnungen die dimensionslos gemachte Dicke ö, der laminaren Une 
schicht auf: e ’ e 


046g 


. (2a). 


;. wobei ug = 7 Ss 
RR 0.3.17. 


; DRS 
1) Vgl. L. Prandtl, Führer durch die Strömungslehre (1942), S. 112— 
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die sog. Schubspannungsgeschwindigkeit ist. Leider besitzen diese Koeffizienten keinen uni- 


Gründen annehmen, daß in Wandnähe die Länge Z außer von der Länge y (Wandabstand) noch 
von der Länge »/v,, die den Zähigkeitseinfluß zum Ausdruck bringt, abhängig ist. Immerhin 
- Jäßt sich für große Re-Zahlen nach den Messungen Nikuradses [2] die Aussage machen, daß 
der Zähigkeitseinfluß sehr gering wird. Insofern wir noch keine ausreichenden Kenntnisse über 


_ die erwähnten Abhängigkeiten besitzen, müssen wir %, u und — U 00 an Hand der Meßergebnisse 
bestimmen. v 


v5 III. Die turbulente Rohrströmung 
Die turbulente Rohrströmung ist das Objekt zahlreicher experimenteller Untersuchungen 


- gewesen. Wir werden uns auf die besonders ausgedehnten Messungen von Nikuradse [2] 
E stützen. , 
E Integration von (1) ergibt 

> BR dp 
Es. Ü 2da Ex 0) 
pr: die Schubspannung ist also linear über Sa Radius verteilt. 
« 4] e 
FL Transformieren wir auf n=1 _ (r, Rohrradius), so erhalten wir mit (2a) und (2b) 
n 0 

Bd 92 

» U U “ 

» ER de 
N, 1.dp r, BU EN a ER 

B 2 do waren DEN as tr’ \dn PR 

wo u, Geschwindigkeit in Rohrmitte, Re = = (d Rohrdurchmesser). 
_ Führen wir die Widerstandsziffer 
Breranl. 0; 
A Suse >44) 
22 fi 
ein, so können wir auch schreiben 
n U,\® u 
. Id— d— 
s) u, PR N SR 
ar Be TR ER ee 1) AT A ERTAER | 1): 
r dn TER: dn 8 ( ) | (6) 


N 


Ei; Auflösung nach Tan ar 


- 


De (Jan 
uf RO TEN ET —dn :.. (M. 


N 6 6) ergibt für di Wand 


# 


VSt N 17 TE% F 
. N u Ns, (8). 


versellen Charakter. Wie Prandtl [3] dargelegt hat, muß man aus dimensionsanalytischen. 
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ein, so folgt für diese aus (8) 


= g . 
5 : (du 
wie man sofort sieht, wenn man 79 = u 15 beachtet. 
. 0 


Wir gehen nun an die Integration der Differentialgleichung heran. i 
Beschränken wir uns auf kleine n, d.h. En 


‘ l | 
EEE al lul 
1—nzl | 2 
(letztere Einschränkung ist nach (5) äquivalent der Annahme einer konstanten Schubspannung], : G 
so haben wir nach (7) 2; 
" I — 
4 2a ’ 
BR Ans urL,; Re: 
| % Re) »:n? Re #an? K 


Hier läßt sich die Quadratur nach bekannten Methoden durchführen und ergibt 


4 PRO Fi 

Ef”: nu. ae Fr LE 2. 42m2 
Pr WR Re ae + Reim) 
Be , 11/A 1/R RENT RE 

Aa EA hy re je PATE 

# -F * Van zes ent + Re sen |+ const 


Setzt man die Integrationskonstante gleich Null, so erhält man m = dafür un=D2Diesone: 


. (11). 


I ; rt RER 
er stimmte Funktion gibt jedoch den experimentell ermittelten Geschwindigkeitsverlauf nur unbe- 
ki friedigend wieder. Die Ursache ist darin zu finden, daß die mathematische Formulierung die 
“u Existenz der laminaren Unterschicht nicht berücksichtigt. 
y Infolge der Haftbedingung besitzen auch die turbulenten Grenzschichtströmungen eine, 
SA wenn auch sehr dünne, der Wand anliegende laminare Strömungszone. Diese Jaminare Unter- 
KaRn: schicht der Strömung ist als frei von der turbulenten Vermischung anzusehen. Der Vorgang der 
Kr ‚turbulenten Vermischung setzt vielmehr erst am Rande derselben ein. Mathematisch heißt das, 
> daß unser Ansatz (10) e 
FR I BEN: 
N. — X , 10 
I; ro N _ - 
>; zu modifizieren ist in wer De 
y% r . _— *(n—n0)?) an a LE. an (2, Si . 
©, ro ee; zn 
ST wo n, die dimensionslose Dicke der laminaren Unterschicht ist. (no ist-allerdings, wie wir noch 

Br: sehen werden, eine gegenüber 1 sehr kleine Zahl, so daß man für nicht allzu kleine 7 unbereb LIE Wer 
Er. n— no = n setzen kann.) D * 
Er Entsprechend haben wir in (11) n durch n —n, zu ersetzen. Die Re rec, ie { 
a nunmehr so zu bestimmen, daß sie die Geschwindigkeit am Rande der. laminaren nern ER; 

| N Se 

: A a 10m) SIEH u Mei, 
ergibt! Nun hatten wir bereits in (8) gezeigt, daß an der Wand rule =, | Tit de 
kleinen Dicke n, der laminaren Unterschicht erhalten ns 


MIEMREA 


%. 16,9 


®) Anmerkung bei der Korrektur: Wie aus einer vor kuteam erschiene. 
Archiv XVIIL (1950), 8. 277 ersichtlich, ist der Ansatz (12) bereits in einer 4 
Mag. (7) 39 (1948), $.7 enthalten. Die in der vorliegenden Arbeit hen: 
im Jahre 1948 Herrn Prof. L. Prandtl schriftlie ES 


we Y s ER ge S X Fran } AN F J | BEN 
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ıser Integrationsergebnis, lautet mithin 


te JRR ee irre 3 ea—ne)) 


u x: Er 


ak v3 in [Re] “(n— +14 Re2 r 2 ea) + BE j . (14). 


“ für sn — m<1l,.. d.h. ern), 1l—n®el 

Er _ f 0 

3 Wir wollen diese Verteilung nach dem Vorgang von Prandtl [4] auf die sog. Schubspannungs- 
geschwindigkeit v, — 70 mschreiben. Diesen in der Theorie der turbulenten Strömungen viel 


gebrauchten Begriff hat Prandtl folgendermaßen eingeführt. Unter der Annahme einer im gan- 


zen Gebiet konstanten Schubspannung erhält man aus dem Schubspannungsansatz = ol8 En ): 


ER 


S Aus dimensionsanalytischen Überlegungen ergibt sich für wandnahe Punkte der Ansatz 1=xy. 
Das führt auf 


T u ei a 
Bu ’ ‚ u. Y d Yy hr 
Integration ergibt 
Bi au u_v, Emy+e). 
- % y 
_ Unter Heranziehung einer zweiten Länge _ , die dem Einfluß der Zähigkeit auf Z Rechnung 
Ü U 
f _ trägt, erhält man, indem man O durch 0, — d: mei ersetzt, die dimensionsanalytisch richtige 
Gleichung a 


3 Ka 7 U BER 
»-yı IE RL} 
Das ist das bekannte logarithmische Verteilungsgesetz der Geschwindigkeit für turbulente Rohr- 
strömungen. Dabei sind x und C', empirische Konstanten. Obwohl nur für wandnahe Punkte 
- abgeleitet hat es sich gezeigt, daß durch dieses Verteilungsgesetz die Geschwindigkeit über dem 
ganzen Rohrradius in guter Näherung beschrieben wird. 

Be: _ Wir rechnen nunmehr unsere Formel (14) auf die Schubspannungsgeschwindigkeit v. um. 
Für die dimensionslose Geschwindigkeit v/v, erhält man nach (9) 


ME, 27 1 
4 > 437 - Ux Ur A ( ) 
Ve E23 8 : 
ur Für den dimensionslosen Wandabstand = an | 
BT, ut, 2 a i am 
RR | Sera my! ESTER AREAL), 
So erhalten wir — mit der Dicke ö, der laminaren Unterschicht — 
Al + Et EEE NN 
ERNE Pr En. BE, 


dor: (8) 


4 


experimentell für wandnahe Punkte ermittelte ER uRT RR 


"heran. Man erhält in bekannter Weise 
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Um einen Begriff von den Größenordnungen zu geben, sei vermerkt, daß für ein mittleres © 


Re= En — 240000 Re = ODER z.B.L.Prandtl[5]). Bemerkung: Reyi— = Reli, 
» 
siehe (38). 
Für die Jaminare Unterschicht selbst haben wir, wie schon dargelegt, 
u: Rei e 
Se für | NSNo: | 
bzw. (19). 
I für Y BrA| 
Y% v VEREIN 


Die Gl. (14) bzw. vr sind ae nicht allzu wandnahe rer einer erheblichen Vereinfachung 


1 EA EEE 
ne Kr 7 > Sy können Far setzen 
e ur 


2 | KR 
Yı + Re g »? (7 — no)" #Re Vz #.(n — no) 


V+ An Eee at 


y v v 


fähig. Fürn— > 


% 
bzw. 


Damit erhalten wir 


EEE ER EN 1/A Rei 
\ re Varafanben Vet dn Er Re ; (20) 
I ee wei 
A { f 
x Re 8 a ! 


dabei haben wir 7 — n, gleich durch n ersetzt, was an dieser Stelle auf Grund der Einschrän- 


kung N—- mM2>— — gerechtfertigt ist, wie wir noch zeigen werden. Entsprechend 


#Re 4 : | | . 


we 1, wy RS et | 
“erde te ’ ie ne ih: 
ur N 
F Rei —. 
1 1 vY 7 8 -. 
REN SENT 


Damit sind wir aber wieder auf das logarithmische Verteilungsgesetz der Geschwindigkeit 
entsprechend (15) gestoßen! Wir haben gleichzeitig eine Analyse dieser bekannten Formel auf 
der Grundlage der Mischungsweghypothese erhalten. - 

“ Als empirische Koeffizienten treten hier auf der Koeffizient x der Mischungswegverteilung 


und die dimensionslose Dicke 20 ge der laminaren Unterschicht. Im allgemeinen werden diese | £ [ 


Größen (s. Einleitung) kennzahlabhängig sein. Für große Reynoldssche Zahlen wird allerdings ; 
der Einfluß der Zähigkeit sehr gering. Wir ziehen zum Vergleich das von Nikuradse m ö . 
eh 


— 5,02 hg Y4r5sı Et. 
Ko AR > \ a f 
RB $ 


nz 


At; h © 
Wir erhalten weiter +2 _ 


’ 


* IS 241 


ERUaAr = .. a I WE ui bi - 
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Wie wir schon angegeben haben, ist für ein mittleres Re— 240000 Re hier — 10000, was 


0 0,0014 ergibt. n, ist also gegenüber 1 eine sehr kleine Zahl, die laminare Unterschicht 
gegenüber der turbulenten Grenzschicht von sehr geringer Dicke. Aus dem angegebenen Wert 
von n, ergibt sich auch die Berechtigung der oben vollzogenen Vernachlässigung von 7, gegen- 
über n, wenn wir die Forderung stellen „— > — BEN 


x 
Re Vz 
Mit den so bestimmten Konstanten x = 0,417 und %4 00 — 7,01 haben wir unsere Vertei- 
v 


lungsfunktion (18) berechnet und in Bild 1 mit der von Reichardt [6] und Nikuradse [2] ge- 
messenen Geschwindigkeitsverteilung verglichen, Innerbalb der Jaminaren Unterschicht, die bis 


Be do _ — 7,01 reicht, haben wir nach (19) 
a a 
DRS 


(18) liefert dann 


— y POKEEEER 
y v Ux 
10 1 9,324 
15 1,176 11,156 \ Experiment 
20 1,301 12,180 0 Nıkuradse,VDJ-Forschg. 3561932) 
30 1, 477 13,447 2 « Reıchardt ZAMM 20 (1940) 

100 16,715 Theorie 

200 | B. 477 19,466 

Während die theoretische Kurve mit der 0, 25 en 7 er 3% > 


von Nikuradse gemessenen Verteilung 
sehr gut übereinstimmt, weicht dem- 
gegenüber die von Reichardt gemessene Verteilung, die in einem flach rechtwinkligen 
Kanal aufgenommen wurde, etwas ab. Wir bemerken zum Erscheinungsbild: Man kann 
hier drei Bereiche unterscheiden. Zunächst die laminare Unterschicht mit rein laminarer 
Reibung. Anschließend daran einen Übergangsbereich, der dadurch gekennzeichnet ist, 
daß laminare und turbulente Reibung von gleicher Größenordnung sind, und der in den Bereich 
mündet, in dem gegenüber der turbulenten Reibung die laminare Reibung keine Rolle mehr 
spielt. In diesem letzteren Bereich wird die Geschwindigkeitsverteilung bei der logarithmischen 
Auftragung durch eine Gerade dargestellt, deren Gleichung wir in (21) wiedergegeben haben. 
Gl. (21) ist identisch mit Gl. (20). Diese hat die Gestalt 


i EV Atmgteomst un ana... 20a) 


BIN EH 
eine Gleichung, die man bei Vernachlässigung der laminaren Reibung auch sofort durch Inte- 


Bild 1. Geschwindigkeitsverteilung in Wandnähe 


ı 
gration der Differentialgleichung (6) erhält, wenn man in dieser setzt- 8x, 1—n=*|1. Wir 
geben noch die Dieken der unterschiedenen Bereiche an: ”o 
a) Laminare Unterschicht 


0< #27  baw. len... 
\ N / Ei 
Re | R 


Für mittlere Reynoldsche Zahlen ist, wie oben bemerkt, 
} Re Vs = 10000, also 0<n<0,0014. 
j Bi Der Bereich gleicher Größenordnung von laminarer und turbulenter Reibung erstreckt 
ihr nach den Messungen (s. Bild 1) bis 2 = =» 50. Wir hatten diese Grenze umschrieben durch 


100 
5, >12. Wirhaben also en Intervall 7< "+? < 50 bzw. 0,014 <n< — 


Rey} 
8 
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Für mittlere Reynoldssche Zahlen ergibt das 0,0014<n<0,01. Also auch der Übergangs 
bereich ist noch relativ sehr dünn. ’ . 
c) Der Bereich der turbulenten Reibung erfüllt demnach das Intervall N ——<sn 3 1; 
Reh 
für mittlere Re-Zahlen N Er erfüllt also fast das ganze Intervall. 


2. Unser bisheriger Integrationsprozeß erstreckte 
sich auf ein Gebiet, das, wie wir gesehen haben, bis in. 
den Bereich der überwiegend turbulenten Reibung 
hineinreicht. Nach oben, also zum äußeren Rand der 
Grenzschicht hin, war das Gebiet beschränkt durch 


0,16 


% Iren, 1—n*=1 .. .(siehe (10)). 


0 


Wir gehen nun daran, die gefundene Lösung fortzu- 
setzen, indem wir den Integrationsprozeß in ein Gebiet 
‚ausdehnen, das gekennzeichnet ist durch 


o Re=105:10° 
e fe =3240:10? 
RE RE 


0 02 0,4 06 08 10 l 
Bild 2. Mischungswegverteilung im Rohr # v_“n +u n° ei are ten Valle (27). 
. 0, 


(Die oben genannte Einschränkung 1—n=-1 bzw. tr, wird fallen gelassen.) Wir ziehen 
also das zweite Glied der Taylor-Entwicklung für //r, in Betracht. Daß I/r, bei wachsender 
Entfernung von der Rohrwand nicht mehr durch x#n beschrieben wird, ist aus den Messungen 


y von Nikuradse [2] bekannt. Aus den gemessenen Schubspannungen und der gemessenen 
| Geschwindigkeitsverteilung ergibt sich nach dem Prandtlschen Schubspannungsansatz 

\ du 2 

| : T= Ruhe n) | x 


die im Bild 2 wiedergegebene Verteilung des Mischungsweges für große Reynoldssche Zahlen. 
(Bild 2 nach Nikuradse [2].) Für größere Re- Zahlen ist die Mischungswegverteilung kennzahl- _ 


unabhängig. Man kann aus dem Bild den Gültigkeitsbereich unserer Annahme er vun ab- 
schätzen. To 

Bei der Integration handelt es sich, da wir bereits in den Bereich rein turbulenter Reibung 
gelangt sind, um die Fortsetzung von Gl. (20). Nach (6) erhalten wir bei SL SREUNE der 
laminaren Reibung die Differentialgleichung 
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Außer durch = = xn-+- un? beschränken wir unseren Bereich nach oben noch 

0 a R a 2 
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Die Integration ergibt dann bei Vernachlässigung höherer Potenzen von n 
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Gegenüber (20) haben wir also in Fortsetzung ein in n lineares Glied erhalten. Di Integrations- 
konstante ist bekannt: nach (20) ist 


124\: Rei 
E ont — 1/44 — Va (2r2e VEl+ en. 
F Man kann Gl. (29) natürlich auch en, a man in die Differentialgleichung (28) mit dem 
t Ansatz 
| u 11/h l “ 
Ä lt {nn —Cn} und FAN HN 


eingeht. Man N dann 
4 
abetnl20w +22) te )+ 0 


was bei ‚Vernachlässigung höherer Potenzen vn n auf = rs S führt, 


In dieser Weise könnte man fortfahren. Wir een: uns at mit diesem Schritt, 
wobei wir hoffen, daß die gewonnene Funktion auch außerhalb ihres Geltungsbereiches die 
Geschwindigkeitsverteilung gut approximieren wird. 

Es handelt sich nun um die Bestimmung des Koeffizienten u/x. Während man den Koeffi- 
zienten x verhältnismäßig einfach und sicher aus der Neigung der sich in der logarithmischen 
Auftragung für die turbulente Geschwindigkeitsverteilung der wandnahen Punkte ergebenden 
Geraden (21) bestimmen konnte, stößt man hier auf größere Schwierigkeiten. Schreibt man (29) 
wieder auf die en v, um, so erhält man 


ee. 1 


vw 1, u 21. L "+ ine Fed 9 (BON 


a 3 2 Nie: Er 
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Das zusätzliche, in v; y/v lineare, Glied 
besitzt also keinen universellen Charak- 
ter mehr. Es gibt bei der logarithmischen 
Auftragung die Abweichung der Tan- o fe=105-103 
gente von der Geraden in erster Ord- Messung Nikuradse VD} 1932 
nung wieder. VREOLIE 
Prandtl [7] hat für die in Bld3 
dargestellte Mischungswegverteilung die 


Interpolationsformel 
: i 4 
2 =0,14 — 0,08 (1 — n)? — 0,06 (1 — n) 
r e (31) 
aufgestellt. Bzw. 
I . 
3 04 0Min HE. (92); 
< 0 
also 


= 0,40, u—= — 0,44 
u 


© fe=3240-10°? 


g | MEER 1,1 ergibt. Wir werden die Messung Nıkuradse VDJ 1932 
5 BE FEESETE, are Theorie 
weiteren Berechnungen mit diesem Wert 


| von u/x durchführen. 
Wir haben Funktion (30) für eine 
verhältnismäßig kleine und eine große 
Be en Zahl berechnet: 
 Re= 105-10° = 0,0180 
R= 3240-1? = 0,0096 | log my 
. (Widerstandszitiern nach Nikuradse) = 
Pe [ReYA= ReyA, siehe (38)] 
din Bi 3 mit den Messungen von Nikuradse e chen Der bekannten, bei größeren 
er AR, e Be - =, auftretenden, Abweichung der Geschwindigkeitsverteilung vom logarithmischen 
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Verteilungsgesetz (21) bzw. (22) für wandnahe Punkte wird dureh das zusätzliche, in » lineare, Bi 


Glied gut Rechnung getragen. 


3. Die Differentialgleichung (6) sowie die Lösungsfunktion enthalten als Parameter die 


Widerstandsziffer A. Diese wird festgelegt durch die Randbedingung _ run n=1, Wan 


unterwerfen auch die von uns berechnete Funktion (29) dieser Randbedingung. Es ist klar, daß 
unsere Funktion, da sie ja außerhalb ihres Geltungsbereichs nur eine Approximation darstellt, 
dadurch auch innerhalb ihres Geltungsbereiches mit einem gewissen Fehler behaftet wird. Inso- 
fern wir (29) der genannten Randbedingung unterwerfen, ergibt sich eine neue Bestimmung der 
Integrationskonstanten derart, daß 


4A Are) Eee. (33). 


Wäre (29) die exakte Lösung im ganzen Bereich, so würden die beiden Integrationskonstanten 
gleich zu setzen sein, was hier natürlich nur innerhalb einer gewissen Fehlergrenze möglich ist. 

Bevor wir diese Verteilungsfunktion mit Meßwerten vergleichen, wollen wir noch die Wider- 
standsziffer, wie es üblich ist, auf die mittlere Durchflußgeschwindigkeit w beziehen. Für w 
haben wir die Definition 


ee 2arr —dı 
Ri 1 
ER 
bzw. mit RSS ” 2 [End 
ro 2 mut, 
Mit unserer Funktion (33) ergibt das ' 
u Ni 7 ze el: 
2=142J3 MENTOR a Re ae .@4). 


Die auf die mittlere Durchflußgeschwindigkeit w bezogene Widerstandsziffer A steht mit der 


auf die Maximalgeschwindigkeit %, bezogenen Widerstandsziffer ja in folgendem Zusammenhang: 


Aw: — Au?; ee, 
es folgt 
me EN. AT (al 
} PRRErSı 
\ u 


und 
1 ; | 
»V< 
le ie, I 
Hr le 
% 8 Eger 
Wir merken noch an 
RER, RER 
v 3 v 


Rey) = Re VA 
Wir haben die Verteilungsfunktion (37) für eine verhältnismäßig ken ng 
sche Zahl berechnet N, a 
De 105 108° = 
 Re= 3240 - 10° 
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(Widerstandsziffern A nach Nikuradse [2].) Die Berechnung ergibt 


ulu, 
ne Re ton or | Re= 3240: 10° 
0,01 0,515 | 0,826 
0,1 0,735 | ..0,796 
0,3 0,849 0.884 
0,5 0,908 | 0,929 
0,7 0.950 0.961 
0,9 0,985 0,988 


0 fe +105:10° Experiment Nikuradse 
» fe -3240:103 VD/-Forschg. 356(1932) 
Theorie 


In Bild 4 werden diese N uteckinktiomet mit den 
von Nikuradse [2] gemessenen verglichen. Die theore- 
tischen Kurven liegen etwas zu tief; insbesondere für 
größere” Werte von n. Dieser Umstand entspricht, daß 


wir für die mittlere Durchflußgeschwindigkeit nach (34) 02 0% 06 08 10 
wu, RA 55 N RER (39) Bild 4. Gesehwindigkeitsverteilung im Rohr 

erhalten, während Nikuradse [2] durch Auswertung der gemessenen Verteilungen 
OU N a EN, (40) 


gefunden hat. 
4. Wir erwähnten schon, daß in der exakten Lösung die Widerstandsziffer A durch die 


Randbedingung 2 == Tsfür in -e:llfestgelegt wird. Ünterwerfen wir (29) dieser Randbedingung, 


1 
so wird diese Funktion, worauf wir schon hingewiesen haben, durch diesen Vorgang mit einem 
gewissen Fehler behaftet. Entsprechend wird die von uns so erhaltene Widerstandsziffer natürlich 
auch einen gewissen Fehler enthalten. 

Man erhält so für die Widerstandsziffer 


Re L 


11/2 A Vs 11/4 
:y5 141m (2»2e]/2 je — 90 + yAllıe) . (41). 


Lösen wir nach A auf, so 


RER rn 1 DIET Rel/s 
2 RN Dan pe 


Führen wir er wie üblich, statt der Maximalgeschwindigkeit %, die mittlere Durchfluß- 


geschwindigkeit w als Bezugsgröße ein, so erhalten wir nach (38) ReyA— Re Vi; ferner nach (35) 
und (36) 


Das gibt dann folgende Gleichung für die Widerstandsziffer 


| x Re/ a 
= 7 (Reya) tz N ls + In ” + — 2 no| - : (43). 
| ix 5 DIR) 

° . Hier ist nach (17) =. identisch mit der dimensionslosen Dicke en der laminaren 
"Unterschicht. Für die empirischen Koeffizienten x, ER ind; #00 do ‚hatten wir erhalten: nach (23) 

und (24) «= 0,417 und id _ — 7,01; nach (32) Fa in LA di 

Damit lautet die Widerstandsformel 

eo Ne (44), 
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 einstimmung mit Meßergebnissen. 
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während sich experimentell für große Reynoldssche Zahlen nach Nikuradse [2] 


ie (Beyi)— 0,55... Be a 
Va SRRERR 
ergibt. ] 
Ohne Berücksichtigung des Gliedes w/x, also bei Beschränkung auf den Ansatz — —=xn, 
hätte man statt —0,50 die Konstante —0,81 erhalten. ro 


IV. Zusammenfassung 


An Hand des Prandtlschen Schubspannungsansatzes wird die Geschwindigkeitsverteilung 
berechnet und eine Widerstandsformel aufgestellt. Als empirische Koeffizienten treten die Dicke 
der laminaren Unterschicht und die Koeffizienten der Mischungswegverteilung auf. Gute Über- 
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Eingegangen am 15. 7. 1950. 


Ein Beitrag zur Berechnung ee 
Von H. Jung in Stuttgart 


Durch die Anwendung der Fouriertransformation gelingt es, verwickelte Knickfälle einfach zu lösen. 
Außerdem konnte auf diesem Weg ein neues Iterationsverfahren angegeben werden, das, wie in einer 
Arbeit noch gezeigt werden soll, sich auch auf die Stabschwingungen bei elastischer Auflagerung aus- 
dehnen läßt. 


By application of ihe Fourier transformation, complicated cases of buckling are easily solved. The 
method leads to a new iteration process that can be generalized to oscillations of bars in a bearings, as 
will be shown in a fulure paper. 


En appliquant la iransformation de Fourier, des cas ments de ambars sont faciles & r&soudre. 
Un rösultat de la methode est un proc&dE nouveau d’iteration, qui peut ätre gen£ralise, comme sera demontr& 
dans une publication future, au cas d’oscillations de barres, dont les supports sont Elastiques. r 
Ilyrem npumenns mpeoöpasosanus Dypre aBTop OyeHb IPOCTO pemaeT CIO3KHBIE cayyan 
IPOrON5HOTO usruba. Kpome Toro emy YAaloch HAÖTN IIPH IIOMOIN 9TOrO IIPeo6OpacoBaHuA 
HOBEIH METON HOCHENOBANENBHBIX IPHÖNSKCHHÄ, ROTOPBIH, KAK HTO ÖyAET IIOKABAHO B AAlb- 
Heümeii pa6oTe, MOKET ÖBITB PHCHPOCTpaHeH Ha cıyyali KkonNebannäi CTEPKHA, JTEsKamero 
Ha YIIPyTOM OCHOBAHHH. } u 
1. Einleitung 
Wird ein Stab durch Kräfte, die in Richtung der &-Achse wirken, a so knickt 
er aus, wenn die Belastung einen kritischen Wert, die sog. Knicklast, überschreitet. Die Größe 
der Knicklast hängt wesentlich von den Einspannbedingungen an den Stabenden ab. Die beiden 
Grenzlagen für die Einspannung sind die gelenkige Lagerung und die ‚vollkommene Einspannung, 
die bei technischen Problemen nie realisiert werden können. 

Es soll aus diesem Grund gezeigt werden, wie durch Momente höherer Drene auf ein- R 
fache Weise die Knicklasten bestimmt werden können, wenn Stäbe elastisch EinRSCDAULE und 
elastisch aufgelagert sind. Kan 

Wie in Ziffer 2 gezeigt wird, lassen sich diese Momente höherer Ordnung in einfacher N 
Weise mit Hilfe der Fouriertransformation finden. Die Anwendung dieser Momente höherer 
Ordnung für Stäbe konstanten Querschnitts wird in Ziffer 3 gegeben. Mit den in Ziffer 2und3 Br 
gemachten Überlegungen läßt sich dann in Ziffer 4 ein neues Iterationsverfahren aufbauen, das 8 E 
gegenüber den bisher bekannten Methoden den Vorteil bietet, daß außer der Iterationsfunktion x 
auch noch die Randbedingungen mitberücksichtigt werden. Dadurch ist es möglich, solche 
Probleme zu behandeln, bei denen die Stabenden beliebig gelagert sind, ‘und der Stab ein \ Far} 
änderliches Trägheitsmoment hat. Anden in Ziffer 5 durchgeführten Bei spielen zeigt es si h 
daß der benötigte Rechenaufwand gegenüber den sonstigen Methoden kleiner is 3. 


!) Auf der Tagung für angewandte Mathematik und Mechanik in Darmstadt, 1950 wurde 
Bericht über dieses Verfahren gegeben. Eh By BURN 
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2. Die am Stab auftretenden Singularitäten und ihre Unterfunktionen 


- Ein elastischer homogener Stab werde durch Kräfte, die senkrecht zur Stabachse wirken 
' und in einer Ebene liegen, beansprucht. Die Eh diese Kräfte hervorgerufene Durchbiegung 


genügt der Differentialgleichung 
LEI) y" al = (a) RR PR a 


p(x) die Belastung, 
Ki! y(x) die Durchbiegung, 
B- SE J (x) das Trägheitsmoment, 
der Elastizitätsmodul bedeuten. 


wobei 


r Die Ausdehnung des Stabes sei — a<xz<b. In diesem Intervall ist auch p (x) als stückweise 
‚ stetig vorgegeben. Multipliziert man (1) mit e** durch und integriert zwischen den Grenzen — © 
und + © so wird mit p(@)=0 für -o<z<—a und b<z<+» 


& + Be 
7 BEI erde Inaetrden... 22... Ra A) 
Die Belastung p(x) des Stabes wird durch die Auflagerkräfte und Einspannmomente erweitert 


und mit 9 (x) bezeichnet. Durch diese Erweiterung erreicht man, daß die Querkraftlinie und die 
Momentenlinie sich durch Fourierintegrale darstellen lassen. Wird in (2) die linke Seite zweimal 


d partiell nach & integriert, so erhält man 

— 22 Se EJ (a) y" et* da— N pla)etedo— Aeita — Beitb iA Myerit®— iv AMzeitb 

’ | —o RG " (2a). 
u Mit der Fouriertransformation 

Me 

% g(A)= [ Gla)eit"de, 

3 +00 ET ED RRSERE EL 


a)= 5. j gA)e-irr dA 


—o©o 


‚erhält man aus (2a) % 
M(a)— ak fr Beeremdedn .. 2.2.2222. 0) 


5 Nur (4) ist ersichtlich, daß zur Bestimmung des Momentenverlaufes es nur darauf ankommt, 
“die Unterfunktion der Gesamtbelastung zu kennen. Die Unterfunktionen von Einzelkräften 
und Einzelmomenten lassen sich einfach bestimmen. 
ne Greift an der Stelle = £ eine Einzelkraft von der Größe P an, so ist die Unterfunktion 
0% 

N Bild 1) 


a PA) — lim, N a ae nn 326) 
3 ae % | : MR | | 
Dr Pi; Definiert m man ein an der Stelle x = a angreifendes Moment durch den Grenzwert (Bild 2) 
Re RR TUR.N, F | ei 3 . 
j Br DIR Fl: nz P | M“almPp 28, @ 
| P N 
E E-tE 
| x 2 
Bild 2 
urn ih (a— e) 
e° 4 te ln three, (6). 
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‚definiert. Damit lautet die Unterfunktion des Momentes zweiter Ordnung 
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% { Ba 


144 Jung, Ein Beitrag zur Berechnung der Knicklasten ir „angen en Ki 


Um zu zeigen, wie ein Moment höherer Ordnung einen Knick in.der elastischen Linie erzeugt, 
werde von der Differentialgleichung der elastischen Linie in der Form ausgegangen 


BEI(e) y' = = Mi. Re N ee (7). 
Wird auf (7) die Fouriertransformation angewendet, so erhält man die Durchbiegung 


oo 


b . 
Mi) , | 
ya)=5 j n ki r Ta EDEN 22.2... BMI: 


Die Funktion M(x) ist durch Momente zweiter und dritter Ordnung zu erweitern, die an den 
Auflagern angreifen, damit y’(x) und y(x) sich durch Fourierintegrale darstellen lassen. Um die 
Momente höherer Ordnung abzuleiten, kann im folgenden J(x) = konst. vorausgesetzt werden, 
da das Verhalten der elastischen Linie nur an der Stelle «—={£ untersucht wird. Damit geht (8) 
über in 


—o 


+» &+e 


ii 1 
u arena. - (8). 


Bekanntlich wird durch ein an der Stelle # = £ angreifendes Moment ein Knick in der Krümmung 
der elastischen Linie hervorgerufen. Läßt man an der Stelle — sein positives und an der Stelle 
C-+e ein negatives Moment angreifen, so erhält man eine Knickstelle in der elastischen Linie. 
Das Moment zweiter Ordnung wird durch den Grenzwert 

N=lmM :2e 


e>0) 


S+e 
N ih wie B 
n(A)= Im 7752: 2 et dE— ah, De NE 


Greifen an der Stelle &—e und ö+e jeweils ein Doppelmoment an, so ergibt sich mit demselben - 
Grenzübergang wie in (6) das Moment dritter Ordnung 


. R=lim N‘ Au, 


e>0 


Die Unterfunktion für einen Sprung in der elastischen Linie wird damit 


A) iA 


R 
Bee ET KALI 


Diese Singularitäten wurden schon von N&menyi?) über die Maxwellschen Einflußfunktionen 
abgeleitet. 

Ist ein Stab an der Stelle —£ elastisch aufgelagert, so wird zit c der Federkonstanten 7 
des Auflagers und A der Auflagerkraft | HEN 1a 


R 
rm m0A, 
EJ() X | 
Mit diesen Singularitäten (Einzelkraft, Einzelmoment, Moment zweiter und dritter Ordnung) : 


lassen sich somit die Auflagerbedingungen berücksichtigen. In Ziffer 3 soll nun an einigen 

Beispielen gezeigt werden, wie sich diese Singularitäten bei der Stabknickung anwenden lassen. 
3. Bestimmung der Knicklasten bei konstantem Trägheitsmoment f£ 

Wird ein an seinen Enden = —a und &—=--a beliebig gelagerter Stab durch ine 

Kraft, diein Richtung der Stabachse wirkt, beansprucht, so erhält man die Knicklasten als BR 
Eigenwert der Differentialgleichung a 


BI-y'=— Poy-Mlo)ı. So 
Die Belastungsfunktion M (x) wird durch die Auflagebedingungen hervorgerufen. Mit der Fourier- a 
transformation erhält man aus (11) Rn 
+ © 1 +a e / W F Kr 0% e y Br: 
&) = — —— M(E aumde, RE 
v5, 6) ‚MON 
EI Fe { A 


f 
—© > 


N NR N&meänyi, Über die Singularitäten der Blastiitätstheorie. 5 
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Aus (12) sieht man, daß es nur darauf ankommt, die Unterfunktion 


+a 
A WR Os a a 9 ur(19) 
zu kennen. Setzt man z 
+» 
1 Fa) 
y(2) = — ODER EN EUR 5 (13a), 
27 12 E 
EJ 
so soll y(x)im Intervall —o<2<--oo beschränkt und quadratisch integrabel sein. Dann 
ist auch die Funktion au: wegen (3) quadratisch integrabel und beschränkt und es ist 
EB 
R EJ 
für alle Az lim ua): N OR ET en . (14). 
Badge u 
EJ 


Die Az sind die Nullstellen des Nenners. Aus (14) lassen sich, wie an zwei Beispielen gezeigt 
wird, die Knicklasten bestimmen. " 


Bild 3 : Bild 4 


Durch Bild3 sei der Belastungsfall gegeben. Mit (6) wird die Unterfunktion, die die 
Endeinspannung berücksichtigt, 


iA (— M, ei4@ 1 M,ei*°) sinAa 


JA)= f j EM a ne (15). 
Geht man mit (15) in (13a) ein, so ist 
r +© 
ya) 4 a ee (16). 
a 
er) 


Mit der Bedingung (14) .erhält man aus (16) 


Een 
Baar de 


nd damit auch die bekannte Knicklast. 


Um die Verwendung der Momente höherer Ordnung zu zeigen, sei ein Belastungsfall mit 
elastischer Lagerung gegeben (vgl. Bild4). Da am linken Auflager noch eine Durchbiegung 
auftreten kann, sind hier Momente höherer Ordnung anzubringen. Die Unterfunktion für die 
Auflagerkräfte und Momente lautet mit (5) bis (10) 


fa) —, (Ae-ika— Aetta -iAM, e-ika_ 22N,e-ika ti R,e ta — iAMzeite) . (17). 
Für die Einspannung am linken Auflager macht man die Annahme 


N=aM,; 3 Rı=%4, 


wobei Cı und c, die Federkonstanten für die elastische Auflagerung und elastische Endein- 
spannung sind. Damit geht (17) über in 


JA) -. (Ae ia — Aeita Lg AM, eite —GAMzeit®— A2c,Mye'**+iAe, Aette) (18). 
Da das Moment beschränkt und quadratisch integrabel sein soll, um es als Fourierintegral dar- 


stellen zu können, muß die Funktion f(A) nach (3) ebenfalls beschränkt und quadratisch integrabel 
vorausgesetzt werden. Damit erhält man für f(A) die Bedingung 


lim f(A)= c, 
1-0 
® 10 
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aus der die Beziehung für die Einspannmomente _ 

M, = M, — 244 | 
folgt. Setzt man noch M,=x4, wobei x ein noch unbekannter Proportionalitätsfaktor ist, 
so erhält man mit (18) und (19) für die Durchbiegung 

+% 


1 N 4A 
2% 2 (2 — 


— 2sin Aa + Ado, erita + 2 Au eit“ 
+ x (Asin Aa— A2oje'?@)Jje-it® dA) 


y(@) = — 12 j“ (20). 


EJ 
Mit (14) folgen aus (20) die beiden Gleichungen 
— 2sinAa+ A?c,cosAa+2AacosAa-+x4?°c,sinAa—=(0, | (21) 
3c,sin Aa-—2AasinAa-+xAsinda — »A2a snda=0 | 


Setzt man in (21) Aa=z und eliminiert x, so bekommt man die transzendente Gleichung zur 
Bestimmung von % 


Zu te Dun de 


ERTL 302 2308 ; Eee TÄLER nr; 
2sin +. a t2?+227 -;, Sin 22 22% „05 22 el (22). 
Ist z, aus (22) bestimmt, so ‚wird IR 3 
P . 
Vr- ” € 
i AR 4. Iterationsverfahren 
“ { Ausgangspunkt sei die Differentialgleichung in der Form 
REN f x ’ 
® | een | 
Ei RB | | 
N wobei u. — E gesetzt wurde. 
AR Wendet man auf (23) wieder-die Fouriertransformation an, so erhält man 
'. pe 
> WER j j y(&) } 
2 legen 
Deh, E 
Setzt man in (24) u—=1, so wird mit %,(%) n ee, ; 
Be + © 5 
; | BERG | 
Pe; = > ih (&E— 2) A & 
en ( nl) 5, I» et. a 65) | 
\ Das aus (25) erhaltene 9,1 (x) genügt bei en Lagerung nur dann Be Randbedingungen, Me 
\ wenn %,(x) schon die Eigenfunktion selbst ist. Da aber y,(@) im allgemeinen nur eine Näherung 
r er der Eigenfunktion ist, wird in (25) eine Fehlerfunktion g,(x) hinzugefügt. Wie in Ziffer 3 gezeigt ‘ 
2 „ wird, läßt sich diese Fehlerfunktion, die nur durch die Randbedingungen bestimmt wird, in | 
“ einfacher Weise so ermitteln, daß Jr Bo 7 
3 UnlE) Aalen se oe 
; e Yn+1(%) = 3 En 72 | I(d) e! ( ©) de dA . . r DD See \ 2 p en y > 
Br! 1a 1 EEE 
j die Randbedingungen befriedigt. Für n— oo muß, damit y,(x) gegen die Er kon- 3 
vergiert g,(x) gegen 0 bzw. gegen einen festen Wert konvergieren. VIERHRE 


Aus der Bedingung (14) erhält man gerade die noch fehlenden Gleichungen, die zusammen Fi 
mit den Gleichgewichtsbedingungen die Auflagerkräfte und Einspannmomente ergeben. Dies rc 
Funktion %n+1(%) kann sehr einfach über die Residuen berechnet werden. 3 


Aus den %,(z) lassen sich die PR ech Konstanten bestimmen 


On, n—1 it Yn Yn—ı 4% Wi % Pe > ze . KAT, Re 


Bekanntlich bekommt man aus den &,,n—ı Näherungswerte für Ds > Eigenwerie, Es 


®») L, Collatz, Eigenwertprobleme. Lelpeie 1945. 


u 
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5. Beispiele zur Iterationsmethode 

An einem bekannten Beispiel soll nun noch gezeigt werden, wie mit Hilfe der Iterations- 
methode die Eigenwerte berechnet werden. „Es wird der einseitig starr eingespannte Stab mit 
einer Last Pan einem Stabende herausgegriffen (Bild5). 


3 Als Ausgangsfunktion wählt man 
z | Ya); | 
die die Randbedingungen noch nicht befriedigt. Mit 
& P=1 wird die statische Gleichgewichtsbedingung ; I ba 
Be MER Wert. (99)... ; ner 
| Die Unterfunktion der Fehlerfunktion ist damit Ba 
ke. 1 
R = [uerası niet] DR nRag): 
Bo NH ) 
Für die Durchbiegung erhält man mit (30) die Unterfunktion 
$ N 3 1 3 
1 2 Wr 
Jı@) = EFE - j Mei: dE+ [ze# dE+- N,e*—iAR, e. 
e & ö N) B Ö . (31). 
28 | Bet; MS, BUN EL Pe 
. - et, Ha) + NaetiaR, er 
; Die Bedingung (14) führt auf Ei Gleichungssystem 
eek N EAREN. 
: 5 } — M+-— + M=0, 5037 Ntk=0 ART ER En ae (32). 
B Als (29) und (32) lassen sich ie Momente höherer Ordnung und das Einspannmoment bestimmen. 
Die Durchbiegung wird rien 
R | #1 | 
= NW) | a ) etz dA 
ü ’ b 2 nd —o : 4 4 a ehe A ht a io (33). 
e 2 © a n 
| | 124.12 6 In 
Führt man mit (33) einen weiteren Iterationsschritt durch, so erhält man 
: Ys(2) = 67 ET Aa ne £ 
 Bildet man mit (33) und (34) die Schwarzschen Konstanten, so erhält man Näherungswerte 
für den ersten Eigenwert. RR: 
TREE f war 
) & 2 in = 2,475 J 
w- B: SI yıyada ; 
Apr Pe 1) IT”, 
\ 3 FR | 
4 
BEN S Yıya d% 
ia=" — -— = 2473.J 
J y,dx 


2® m durch Abb. 6 Benenen Beispiel wird nun noch das. Iterationsverfahren bei nicht 
stantem Trägheitsmoment durchgeführt. Als Ausgangsfunktion werde gewählt 


ER ng Bo nr nr „für —1<sız<s0, 
) 2 nur « eine sehr rohe Näherung, die =; ee tadingungen noch nicht befriedigt. Die 


ah eo a Ne ae ere, ‚mA or en (35). 
ktion Io} wird hierfür & 25 . 
ao (36). 


0<sa<si1 
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Aus (26) erhält man mit %,(%) und dem Verhältnis 9 — Ji: Jo Er BR. 


N A 


na = 7,9 [Tau +9 Murrae — Bao + more 


® 2 2.{87% 
| fa: | 
0 E = 
Wird in (37) die ae über & durchgeführt, so man 
1 1 eit ee —i4 eik eih Ms; | a 
ger, hr pr 2 ar ee n,08 (273) 


Bi Rn Bild 6 


Ist der Stab TE gelagert, so sind in (37a) noch Momente höherer Oral nach VB 
Er und (10) hinzuzufügen. Aus der BE (14) folgen die beiden Gleichungen 
EN AURTB 40: B FE Be 
Mit (35) und (38) lassen ach die unbekannten Auflagekräfte und Einspannmomente berechnen 

Die Durchbiegung ist durch die Unterfunktionen EN 


0) A 
AN 7 ln a En el 4 Tea) 
F "M eit a 2 a ! 
re “2 72 34 TER, En B ar 
Denn Wertet man = intkekäe . 
a ae 
“mit Hilfe der Batanen aus, so erhält Son ERS, = rt EEE, e 
Bit: A (+ 2)’ — zur 018250 2 
: Yı(2) = B ; de 2% x 5 
Bf ae), a N a 1] RL HE u a 
Die a Y6(2) und Y(@) berechneten re Konstanten sind i in Tabelle r zusammen. 
SH ey; N EI Fa Ba = 
Tabelle 1 RE EEE ET 


ur, 


rip . Y 44 h. 7 + { 
1 = a En y er 11 
er pe pe 
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KLEINE MITTEILUNGEN 


Die Summe der reziproken Potenzen der Null- 
stellen von M,, „2 (2) hinsichtlich =. 


Im folgenden soll eine Gesetzmäßigkeit hergeleitet 
werden, die für die Summe einer beliebigen negativ 
ganzzahligen Potenz der unendlich vielen Nullstellen 
von Mx, „2 (2) hinsichtlich z bei festen Werten von x 
und besteht. Es dürfte zum besseren Verständnis 
erwünscht sein, den Beweisen dieser Formel zunächst 
einige allgemeine Tatsachen über die Verteilung dieser 
Nullstellen vorauszuschicken. 


1. Allgemeine Aussagen über die Nullstellen 
von Mx»,uj2 (2) hinsichtlich z. Es wird hier unter 
dem Symbol Mx,„/2 (2) die durch die Gleichung 


ul u; :)|ru+m 


‚2 .gral ut +2; 14m.) 74m 

& a) 
definierte Funktion verstanden, die abgesehen von 
dem Faktor z(1 + „i2 für alle h = E 6 #0,—], 
— 2... 


eine in z ganze transzendente Funktion und 
bekanntlich die um die Stelle z=(0 herum existente 
Lösung der Wh.D.Gl. 


a RT! 
| el 


7 (da) 


d?y 
FE 


ist. Für R(p) > — 1 zählt daher zu den Nullstellen 
von M»,uj2(2) stets auch die Stelle z=0. Diese 
triviale Nullstelle bleibt im folgenden außer Betracht. 
Wir erwähnen noch als eine weitere Folge von (1), 
daß wegen der reinen multiplikativen Verzweigung 
der Funktion an der Stelle z=0 jede Nullstelle des 
Hauptzweiges aus dem Bereich -n< are?) <S+r 
auch eine Nullstelle irgendeines beliebigen anderen 
Zweiges dieser Funktion ist. 


Aus der D.Gl. (1a) folgt in bekannter Schlußweise, 
daß alle evtl. vorhandenen Nullstellen nur einfach 
sein können. Außerdem müssen die evtl. Nullstellen 
der Funktion in unendlicher Zahl auftreten, denn sonst 
wäre die Funktion Ar ef 22: M, „2a (2)), die für 
|2]>®im ganzen Winkelbereich |arc (2)| Sr ver- 
schwindet, eine rationale Funktion, was sie aber auf 
Grund ihrer asymptotischen Entwicklung gewiß nicht 
ist. Die unendlich vielen Nullstellen können dann 
aber nach eine ekannten funktionentheoretischen 
Satz nur den Punkt ® zur Häufungsstelle haben. Aus 
den Rekursionsformeln, die die verwandten Funkti onen 
M „21/2, uj2+1 mit M,, „j verbinden, und aus der Ein- 


fachheit der Nullstellen folgt ferner, daß die Nullstellen 
dieser Funktionen von denen der Funktion Mx, „2 (2) 
selbst stets verschieden sind. Aus der bekannten 
Beziehung | 
ni 
el Pi) 
M_ un? Bere 


BL) 


"M, ua 


‚geht schließlich noch hervor, daß die Nullstellen von 
M x, uj2(2) entgegengesetzt gleich denen der Funktion 


-M—. x, uj2(2) sind. 


Da die Nullstellen von My, uj2 (2) bis ins Unendliche 


_ ausstrahlen, so vermag die asymptotische Entwicklung 
_ dieser ion für große Werte von z sehr wohl eine 


wr .. 


Information über die Lage der großen Nullstellen . 


von Mx,uj2(2) zu geben. Aus ihr folgt mit den Ab- 
kürzungen 


ee 


x (3%) 

mit B=0 für I#1u>0 
mr tir (a0)... (38), 
EEE (37), 


daß sich die weit entfernten Nullstellen z=x,„ der 
Funktion M%,uj2(2) für n — oo bei reellen Werten 
von u — 1 nach Betrag |z„| und Phase 9, aus den 

beiden Formeln s 


—1 
[znl 2%, In |m| m x: (20+°5>)4 2 4a) 


mn rl z— 2 inll-2a)/lel) (4b) 


: 37 m 
ob. Vorzeichen: — 7 <m<H ee 
unt. Vorzeichen.: — 2 <m<Ht BR 


berechnen lassen. Die Nullstellen aus dem Bereich 
— nn <oare(2) <S-+n können also im Falle eines 
beliebigen komplexen x in zwei Gruppen zerlegt wer- 
den. Die Nullstellen der ersten Gruppe nähern sich 
asymptotisch der negativ imaginären Achse der z- 
Ebene, und zwar von rechts oder von links her, je 
nachdem x, 20 ist, die zweite Gruppe von Null- 
stellen unter dem gleichen Verhalten der positiv imagi- 
nären Achse. Zwei demselben » entsprechende Null- 
stellen beider Gruppen liegen aber für ein x,==0 nicht 
symmetrisch zur reellen Achse. Erst für ein reell- 
wertiges x bilden zwei solche Nullstellen ein konjugiert 
komplexes Paar. Der Abstand zweier aufeinander 
folgender Nullstellen nähert sich mit n— © dem 
Wert 2x. Ist x,—=0 und also x=i-r mit 77.0, so ist 
nach (3%) auch «=(0, und es liegen dann zum min- 
desten die weit entfernten Nullstellen in diesem Falle 
genau auf der imaginären Achse der z-Ebene. Daß 
diese Feststellung unter der angegebenen Beschrän- 
kung sogar für alle Nullstellen zutrifft, kann man 
leicht nachweisen, indem man von der D.Gl. 


tn (Fett a © 


e 
(T=m+ it,) 


für die Funktion Mir, „J2 (— i£) /(— iZ)U2 ausgeht. 
Von den damit für ein reelles r selbst als reell erkannte 
Nullstellen der Funktion Mir, „2 (— i£) läßt sich 
dann noch weiter zeigen, daß jede dieser unendlich 
vielen Nullstellen mit wachsenden Werten von r an 
Größe abnimmt. 

Die allgemeinen Angaben mögen noch der Voll- 
ständigkeit halber durch die Ergebnisse einer Arbeit 
von A, Kienast [1] ergänzt werden. Es gelang ihm 
nämlich in der genannten Arbeit nachzuweisen, daß 
die Funktion. Mx, „p2(2)/2(1 +#)/2 bei reellen Werten 
von x und u zwar keine positiven, wohl aber n+1 
negative reelle Nullstellen besitzt, falls mit n = 0, 1, 
2,...—_n—1<(u+1lj]2+x<—n ist. Hin- 


IA TO EIN 7, 


‚gegen hat sie keine negativen, wohl aber n+1 posi- 


tive Nullstellen fürein —n — 1 (1+ u) —r»<—n. 
Dabei ist eine wesentliche Voraussetzung, daß u>—1 
ist. Auch über die Zahl der Nullstellen im Falle eines 
u<—1 vermochte er präzise Aussagen zu machen. 
Sein Beweisverfahren nutzt die Tatsache aus, daß die 
Funktionen K, K,),... Km mit Km = Mx -+m,u[2(2) 
-z—i+u)2 mit m=0,1,...,n eine Sturmsche 
Kette bilden. Diese Angaben dürften zum Ver- 
ständnis des Folgenden ausreichen, 
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2. Die Herleitung der Summenformel. Wir 


gehen zum Beweise der Summenformel von der 
Beziehung 

(0 +) 
wir Maus +1) dv 1 ‚Mu +1(2) 
2ri Mx,uja (vo) v(w—z) : M x, ul2 (2) 

Rn 

Mx,uj2 +1(Q,) 1 () 
_ M’x,ul2(a)) a,(a, —2) 


mit a, = ay/x,u) 


aus, die eine unmittelbare Folge des Residuensatzes 
von Cauchy ist. Das links stehende Integral, dessen 


‚ Integrand auf Grund von.l. (1) eine in v eindeutige 


Funktion ist, ist längs eines geschlossenen Kreises 
mit sehr großem Radius A, erstreckt zu denken, der 
sich zwischen den Nullstellen @,, und ay , | hindurch- 


N 


zieht. Auf der rechten Seite von (1) sind die Null- 
stellen a, von M, „2'®); die selbstredend von x und u 


abhängen, so zu zählen, daß für A — 1,2,3,....]a, | 
An Hand ur symp- 


sie ,,1Sie gl 


“ist. 


totischen Entwicklung von Mx, „j2 (v) für |v] > po läßt 
sich dann ohne Schwierigkeit zeigen, daß das in (1) 
unter dem Integralzeichen auftretende Verhältnis der 
beiden M-Funktionen für v— co und für alle Jare(v)| S 
z stets beschränkt bleibt. Für A„ > ® folgt dann 
aus (1) die es 


Mx,uj2 +1 (2) 


SM %,ul2 +1 


(a;) 


2 / 


M x, u]2 (2) 


A 


Mau (a) 


\ a,(2—a,) 
Nun läßt sich RER aus den Integraldarstel- 


(2). 


lungen, die für die Funktion M,, „2 2) bekannt sind, 
leicht die Rekursionsformel 


I+ u) ı 
(1+ u) Miu (2) = en = . "Ma, u]2\2) 
A In 
Ga Far ar Mx,u]2 +1 (2) 

y r u 
aufstellen. Sie kann nach Division durch z und 
Mx, uj2 (2) nach Zusammenfassung geeigneter Re 
auf die Form 
we { auisE ++ 2 „Mx,ul2(2 ve 
dz Pi u; a +2 

1\2 
(= )* M (2) u 
a BET EEE AT —1 
‚trau M x, ul2(#) BR 
gebracht werden. Aus (3a) folgt für jedes a 4: Es ist 
\ MN l+u 
Mu, ule+1l 0), ug (ar) = ( 


und daher in Rücksicht auf (3b), (2) und 1,(1) 


dz 


A 


NE Su 
2 


oo 


=1 


1 


2 ———— 
ar i+n 
En e Bi EU 


[4 
a’? — a,) 


= Ihe) 


(4). 


u>—]) 


Aus (4) entsteht durch Integration zwischen den 
Grenzen 0 und z nebenbei die zuerst von J. Horn [2] 
angegebene Produktdarstellung 


Mu,ufe (2) = HF @Rr(I+ u) 


.erzzi ji +, 


oo 


4=1 


2 


ZU-,, 


a). 


- 


2a} l;; 
)- 4 1 (4a), 


Kleine aa sea | 


Nondern Vorne Falle auch unmittelbar hätte : 


ausgegangen werden können. BEE 
Mit der ständigen Abkürzung 3° 
SH = E p=2,3,4 . (5) t 
Fe 
ist nun für alle |< ja lsIe|- 000 
© R oo k ar 
aD: n+1, ar 
Deren D: Sn t2 (6a). 
A=1 f n=0 N 
Andererseits besteht die in einem Kreise von nicht 
verschwindendem Radius absolut konvergente Ent- 
wicklung 
Le 


» “a 
In f(2) = I Zn Den (6b) 3 
leo p=ıP > 3 
mit en 3 ‘ 
= (Pt ; 
b, 4 
1 LScHE ee. d 
b, RE u 
Den = Te, 0 
Die, td, Ä j 
EUR (7) Fi 
_ndb-1 b—2 b, 2 
Mn) Er ee = 
dd dbv—ı &v—2 b; but 1 
p In b, bo. .o.o.. b, b, e . x 
Setzt man im vorliegenden Falle j ” 5 


al ;1+% >)» 

so ist mit | 
(= l(uat+ly=la+1)(u+2)-wrl+p—I) 
RN EL UNE R 


m 


b NER N RN, S # _ ) A use ER 
BE 
+14 t = ee 
(#- P) ra e=h28. Kay 5; 5, , 


4 ER 


und es folgt (4), (5) Eon, (6b) Er 
PR ya #=23,4.). ST: 


Löst man die Deinste (7) von Hebrecheng lern ar 
Ecke her auf und berücksichtigt ie die ‚Gl 
ra man zu der a LA, 


zu 


Er 


nr "% Ba: Di "R ”; 
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. (10a) 


2 ar Dia 
| Mer) e-® 
S, = 
Wr AH ar 
Mac 
S 


S(u+2)?(u+: 


a. . — 
“ x + 
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ers am: A: 
°(u + 3)(u + 4) 
) Ie- ey) Tu 4 19) + + Diw+ 2) 


3) (u + 4) (u + 5) 


EN erene- 
S- 


= Ar dc ehe 22; me Fee FE dr 
En 
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anf ar) 


= 


Ur) (100) 


(u + 1) (u+ 2)?(u + 3)?(u+ 4)(u + 5)(u + 6) 


und darin ist - 
A=2-[21u° + 181u2 + 5134 + 473] 
— (u + D%(u + 2)(14u2+ 850 + 131) 
3. Allgemeine Gesetzmäßigkeitenim Aufbau 
der 8, und Sonderfälle. Wir besprechen noch in 
Kürze einige allgemeine Eigenschaften der reziproken 


Potenzsummen Sy), wie sie in den angegebenen 
speziellen Bormeln bereits deutlich hervortreten. 


a) Alle Potenzsummen $, sind gerade oder ungerade 
Funktionen von x, je nachdem p selbst gerade oder 
ungerade ist. 


Diese Eigenschaft folgt unmittelbar aus GL. (2) 


und 2. (5). 
2 

ß) Alle S, müssen den Faktor x? — £ = ) auf- 
EEE nach GL. 1.(1) 
(2) = + me. tar 
‚ ul2 { 


weisen Da nämlich für «= + 


ist, so hatin diesem Falle die Funktion M. x, „j2 (2) keine 
im 'Endlichen gelegenen Nullstellen. Daßin den Gl. 2. 


2 
(10&...e) der mit dem Faktor x? — N & 


plizierte Term im Zähler der Brüche ebene stets 
so geschrieben worden ist, daß darin Potenzen von 


Kam Ba) 


multi- 


auftreten und nicht die Potenzen von x? 


F allein, ist freilich eine willkürliche Maßnahme, die aber 
in bezug auf die Gestalt des letzten Summanden zweck- 

F mäßig zu sein scheint. 

- y) Der Nenner Ny(u) von Sp hat die allgemeine 

2 


Form: 
Not) = (u + um. (ur 7) 
De (u + p)lolm 


Darin bedeuten die [p/n] diejenigen größten ganzen 

positiven Zahlen, die eben noch kleiner oder gerade 

‚ en p/nsind. Es sind also in (1) stets [p/1]= p und 

[p[p] = 1. Dieses Gesetz wird durch den Aufbau der 

‘Nenner in den obigen Gl. 2. (10a...e) nahegelegt und 

läßt sich dann mit Hilfe von 2. (9a) durch den Schluß 

von n auf n+1 auch leicht allgemein beweisen. 

5 Auf Grund dieser Ergebnisse kann demnach allge- 
Sin, gesetzt; werden: 


| ar Ui: 2% er) er (2a) 
' „4 T ) .) 


. das Laguerre-Polynom L) (2). 


6) In diesem Ansatz kann über den Summanden 
mit dem Koeffizienten a(@ von vornherein eine ganz 
präzise Aussage gemacht werden. Dividiert man näm- 


lich die Gl. 2. (9a) durch AT — x und läßt dann 


Du Lee 


1 

Sr gehen, so erhält man die Relationen 
lim Sp Kr 1 
a1 + L gie 

en | Be uTı ers = (u + 1)» 


S; ; 1 
3 x ans =(u+1)(u+ 


277 .: 


. (3a) 


1)» (3b). 


“Ho 


Führt man aber denselben Grenzübergang an den 
Gl.(2a,b) durch, so ergibt sich schließlich im Hinblick 
auf (3) das Formelpaar 


a(@p) — Na (w) ! (43) 
(ur 1)(u+ 1), 
Na 
ep+n _____Nep+ı(M) 
Bi (“+1)2-(u+lep+1' (eh). 


Auf die nämliche Weise würde man auch für die 
höheren Koeffizienten ar "zu allgemeinen Darstellungs- 
formeln gelangen. 


Bemerkenswert sind die Sonderfälle x= 0 und 


mEL mit m=1,2,3,... Für«=0 wer- 


“= m-+ 


den nero die Nullstellen von M,, a/2 (2) identisch 
mit den Nullstellen der modifizierten Besselschen 
Funktion /ay2(2/2). Es ist also in der üblichen Be- 
zeichnungsweise in diesem Falle 


us H2i ig A=h2 .-.)] 
iin AL.) 
In der Tat liegen jetzt für «= 0 alle a, nicht nur auf 


der imaginären Achse der z-Ebene, sondern sie lie- 
gen dort noch obendrein symmetrisch zum Ursprung. 
s leuchtet daher ein, daß im vorliegenden Falle die 


EICHE 


‚ungeraden Potenzsummen verschwinden. Die geraden 


Potenzsummen S;,, S, S, der Gl. 2. (10a, c und e) aber 

gehen dann in die zuerst von Rayleigh [3] auf- 

gestellten Formeln für die reziproken Potenzsummen 

der j,,, über. 

I+u 
5, 50 entartet be- 


Wird andererseits «= m + 
kanntlich die Funktion Mx, „/2 (2) im wesentlichen in 
In diesem Falle 
bestehen also die Sy jaweilet nur aus p Summanden, 


(2b) “ ‚and man rechnet leicht aus, daß z. B. für 


in 1 


Ba mezi EN NA 
e 5 (1-+ u)? En 
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(6b) 
(9==72, 3, e) 


ist. Auf diese Weise ergibt sich eine einfache Kontroll- 
möglichkeit für die allgemeinen Summenformeln. 
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Ein übersichtliches Schema für die Matrizen- 
multiplikation. 


Für die numerische Rechnung ist es häufig be- 
quemer, zwei zu multiplizierende Matrizen X und ® 
nicht nebeneinander, sondern nach Art des folgenden 
Schemas anzuordnen: 


[8] Önt 
[A] [AB]; [ajı @j2 a7; * ajn] Er | 
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mit den gleichen charakteristischen Zahlen wie X und 
zeigt, daß die Transformationsmatrix Tim allgemeinen 
die Form 


haben muß. 
Schreibt man an Stelle von T-AT=® nun 
AT—= TP und ordnet dies in der folgenden Weise an: 


[®B] 
[ T I[AT] 
[A] AT] 

so erhält man zur Bestimmung der Elemente t;; uns 
P»u (soweit sie verschieden von O0 oder 1 sind) dad 
folgende einfache Schema. s 

Die Matrizen AT, ® und T werden in der durch 
die Pfeile gekennzeichneten Reihenfolge spaltenweise 
ausgefüllt. Zunächst bildet man die erste Spalte 
von X Tin ® (sie ist gleich der ersten Spalte von W) 
und schreibt sie noch einmal in ® ein. Dann füllt 
man die erste Spalte von ® aus, danach die zweite 
Spalte von T. Jetzt kann man die zweite Spalte 
von AT bilden, darauf die zweite Spalte von ®, dann 
die dritte Spalte von T. usf. Zur Bestimmung der 
noch unbekannten Elemente in den Matrizen ® und T 
sind immer nur einfache Additionen, Multiplikationen 
und Divisionen erforderlich. Ist ® gefunden, so läßt 


Pıı Pıa Pıs *'" Pın (3) 


1 Pe2 Pas ''" Pan A 
Re ) 1 Pas °"" Pan | 
(4) RS ones. 


wenn 

1 0 0 ) Ay terre Gınkan 

) RER N) 0 Green Ganlan ' 

0) lyo bag 0) Ası oma nich ben Agnlnn 

Y 0 Ina Ins Inn U EEE ER Anntnn ! (2) 

4yı Aja dia dın Aıı .er Tree Amnlnn =, 
Ası (ga (gg den Ayı SU@ late ale,e Dualen® aanlnn 
Qzı Aygg As; Azn Da Re at Agntnn 
Anı Ang (nz Ann Ar EB Anninn (1) 


Jedes Element der Produktmatrix AB = € wird 
gefunden durch Bildung des inneren Produktes der 


„ links.stehenden Zeile von X mit der darüberstehenden 


Spalte von 8. Man sieht bei dieser Anordnung un- 
mittelbar, daß A ebensoviel Spalten wie B Zeilen 
haben muß, damit die Multiplikation überhaupt aus- 
führbar ist. 


Auch bei der Multiplikation von mehr als zwei 
Matrizen ist die Schreibweise zweckmäßig. Sie läßt 
erkennen, daß jede Matrizenmultiplikation sich in ihre 
einzelnen Zeilen (Spalten) auflösen läßt, sofern man 
mit der Matrix ganz links (rechts) zu multiplizieren 
beginnt. 


Besonders angebracht ist das Schema bei dem Ver- 


fahren von K.Hessenberg zur Aufstellung der 


charakteristischen Gleichung p(k) = |XU—k€|] einer _ 


gegeben quadratischen Matrix X. Hessenbergtrans- 
formiert die Matrix X durch Multiplikation mit einer 
geeigneten Matrix T von rechts und mit ihrer Rezi- 
proken T-! von links in eine Matrix der Gestalt 


Pıı Pıa Pıis Pu ''' Pun-ı Pın] 
DPaa Pas Par '*'" 
?=|0 r be Pa '"" Pas n-ı Pan 


KT ee 


. und außerdem das Element %4o,»+o durch eine N F 
Fall 2.)- Sind alle „= 0(u=1, 2, ...n, so 


sich p(k) = | U — kE| = |® — as rekursiv 
auflösen, indem man der Reihe nach die Hauptunter- 
determinanten D,, D,, ..., Dn ermittelt, wo Do, die 
durch Fortlassen der letzten n — o Zeilen und Spalten 
entstehende o-reihige Unterdeterminante bedeutet. 
Sind die %, berechnet, so werden ebenso einfach die 
Eigenvektoren t» zur Matrix ® bestimmt und mit 
Hilfe von | | 
$ b=ry, 

zurücktransformiert. 

Es kann nun vorkommen, daß bei der Berechnung 


von X ein Element der Hauptdiagonale gleich Null 


wird. Män hat dann zwei Fälle zu unterscheiden: 
1. Mindestens ein Element der v-ten Spalte von T. 
ist #0. 
2. Alle Elemente der v-ten Spalte von Tsind = 0- 
Im ersten Falle wird die Matrix T, im zweiten Falle 


die Matrix ® abgeändert, und zwar: 


Fall 1.) Ist 1» = 0 und fy+o,, das erste von Null 
. verschiedene Element in der v-ten Spalte von T, so 
- ersetze man die Null an der Stelle v, v» + o (also ober- R 
halb der Hauptdiagonale) durch ein vorläufig noch » 


unbekanntes, aber von Null verschiedenes Elem 


eine 
un 
u‘ 


ersetze man die Null an der Stelle 4 durch 


Legt 
- 


ie ae 


q 


f 
. 
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und in der Matrix ® die 1 an der Stelle p»,»—ı durch 
eine Null. Die Determinante |® — kE| zerfällt dann 
in das Produkt zweier einzelner Determinanten und 
ist besonders bequem auswertbar. 

Ein einfaches Beispiel diene zur Erläuterung: 


130102 21702 
1221,34] 
0 1-2-5 
0.0.0.3 
80:04 0 2. 042 
0-100 -1-2-3-1 
0192.71.0 2.5.4-3 
0:=1::0” 1 17,2-.374:|NUX 
DJ 272170.%8 22.12:0% 2% Kontrolle: . 
-1-5-3 -1 -1-2-3 -1 Spur X = Spur B 
20 4-3 25 4-3 u 
178 374 1.234] AT 


Bei der Berechnung der vierten Spalte von T ergibt 
sich mit den ursprünglichen Matrizen T und ® zu- 
nächst t,,= 0, so daß T-! nicht existiert. Man hat 
daher entsprechend der obigen Fallunterscheidung T 
und ® durch Setzung von t„—=1 und 9,=0 zu 
modifizieren. Jetzt bildet man 


IB -RE|=(k—3)(B® — 2 —8k+12)= 0 


und errechnet zum Beispiel für k=3 zur Matrix ® 
den Eigenvektor 


-2 2 y 
SR) N, 0 
y=| j |;nunistg= Tu, folglich: | ] 
-1 -1 

Sales 03.040 2 |t 
0-1.0:0 0 
07.277150 1 
U... -1 


Hier ist zufällig ieh 


Hannover, Sigurd Falk. 


Ist eine Revision der Gewichtskoeffizienten- 
Berechnung in der Methode der kleinsten 
Quadrate erforderlich ? 


(Eine Stellungnahme zum Einwand von 
S.L. Piotrowski) 


Kürzlich ist in einer amerikanischen Publikation!) 
der Nachweis geführt worden, daß die herkömmliche 
Berechnung des Gewichtes der Unbekannten in einer 
Ausgleichung nach vermittelnden Beobachtungen mit 
Linearmachung der Fehlergleichungen nach dem Satz 
von Taylor ein nur unvollkommenes Ergebnis liefert. 
An einem Beispiel wird dort gezeigt, daß die übliche 
Berechnung zu den Gewichten p, = 1:0,20 und 


py —= 1: 0,060 für die Unbekannten & und y führt, 


während die genauen Werte p, = 1:0,09 u. py=1:0,028 
betragen. Diese Abweichungen sind so beträchtlich, 
daß man die Frage stellen muß, ob z.B. bei den 
grundlegenden Aufgaben der Ingenieur-Geodäsie wie 
auch der staatlichen Vermessungstätigkeit, nämlich 
bei den trigonometrischen Punkteinschaltungen, die 
eines der Hauptanwendungsgebiete der Ausgleichung 
nach vermittelnden Beobachtungen darstellen, eine 
Revision der seitherigen Rechnungsmethoden not- 
wendig wird, und wie diese dann-wohl durchzuführen 
wäre, 
? r 


Zum besseren Verständnis soll zunächst der Grund- 
gedanke von Piotrowski in etwas abgewandelter 
Form kurz umrissen werden: Zwischen den um v; 


2) $S,L. Piotrowski: Some remarks on the weights of un- 


_ knowns as determined by the method of differential corrections. 
j Pte Acad, Sciences, Vol. 34, Nr.2. Harvard Reprint, 304 
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verbesserten Beobachtungswerten L; und den (zwei) 
Unbekannten & und y bestehen die Beziehungen 


nra=laay, HEL23.NFR. 
Mit der Zerlegung 

ı=%-+4%;\ (2) 

Yy Be Yot Ay j re ‚Iamley a Teheran =); 


worin x, und y, beliebig wählbare Näherungswerte 
sind, ergeben sich nach dem Satz von Taylor die 
folgenden Fehlergleichungen: 


worin 
aM) PR } 
7 0x ’ ı— oy 3 Ka Un, Yo) . (4), 


Unter Benutzung der Gaußschen Summen-Symbole 
ZeB, 


(8), 


a; 


n 
5 IB Fr) 
[vv] für z v; 
= 


folgen aus 
[vv] = Minimum 


die Normalgleichungen 
lau. 057 [850] 0 
oder 
[aa] Ax + [ad] Ay+ [ao -—D)]=0 | 
[ad] Ar + [88] 4y+ lo —-DJ=0J 


Die Gewichte pz; und p, der Unbekannten x und y 
sind bekanntlich 


pa = 1:[00]; 


a 


(6). 


?y=1:[ßPl . . . (7a)» 
worin 

VL ST 
024; Tag 
Sieht man in (2) die &,,y, als konstante Teile an, 
so erhält man aus (6): 


04x Age HN 
[aa] Fi .“ [eb] 57, =4j 


RT DZ N (7b). 


oder 
[aa]a; + [ab] P;=a;, 
P) i 84 Be 
Ax AU = 
Lab]. z7, + Ne _ b; 
ocer 


[adJa/ + [68] = bj 


Nun denke man sich, um die Einflüsse der in (3) ver- 
nachlässigten Glieder zweiter und höherer Ordnung 
zu eliminieren, die Ausgleichung mehrfach wiederholt, 
indem die sich ergebenden Werte &, y immer wieder 
als neue Näherungswerte z,, y, angesehen werden. 
Im Endzustand wird dann sein: 


uf! u=-/M-L.... 
so daß aus (5): 
[a(f— L)]]=0 und [b(f-ZL)]=0.. (10). 


Die Differentiation nach Z; liefert jetzt [im Hinblick 
auf (#)]: 


M0ol+ alle SL + 


0 
+ lab) + Yu DW 5L, = 9; 
er Zi 
0% 
b ltr Dir —— 
tab] + DM z7, + 
++ en DSL; =D; 
worin 4 ” 
1 = EN 
fee = dx: 17) ne öy: Szu — dx 0y (22) 


® 
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(11) geht in i. ) über, wenn 
Be ea | 
02: Nagy alas). 


Falls (13) nicht erfüllt ist, dürfe zur Weiterrechnung 
nicht (8), sondern es müsse (11) verwendet werden. 
Soweit Piotrowski. 
II. ; x 
Für die Weiterrechnung lösen wir (8) nach &5 und ßj 
auf; es ergibt Br 


= Zatl bled; | 


D . (14) 
P;= z'dile a] — ajlab]) 
mit : | 
—=[aa] [bb] — [ab ..... (15), 
woraus a 
LI ER! ed _ = 
[Ha == und [#’ß’] = PR (16). 
Ersetzen wir in (11): 
[4A] =[aa] Hofes], [ABI=Tabl + ofen (1m) 
[BB] = [6b] + [ofyy] 
so hat man entsprechend (14): 
=, (a [BB] — v[ÄB)) | 
und . (18) 
= 5 0LAA- aldm) | 
mit 
Di= [42182] 1A BR mer (19). 
Dann habe 
[& 0] = Tr its a) [BB] + [66] [A BY 
— 2[ab] [AB]BB], | 
(20). 
BBI= 5; =, (lbb][A AP + [aa] BB? | 


— 2[a5] AB][AA)) 
Setzt man jetzt in (19) und (20) die Beziehungen (17) 
ein, und vernachlässigt man die Glieder von der 


Ordnung v? als kleine Größen gegenüber den v-Glie- 
dern, so erhält man 


D=D+AD,AD= [aa] [ofyy) + [55] [ofer] 
Le) 


] a ® zu] 
: 6) 
„[&%6] = (Dr aD% (1db]+ 2 [vfy)), 
D AH NBe), 


[? Bl = (Dr AD) (laa] + 2[vfzx]) 
Mit der Bedingung (13) geht (22) wieder in (16) über. 


Mit einem hinreichenden Grade von Annäherung findet 
man für (22) auch: i 


Bal=z{t s 27) [05] + 2{ofy)): 


SCHE: 5 . (23) 
BA=7 ( —=2 7 0) + 2[vfzs]) 
bzw. 
[aa] = 5: &]+ ag [Yo] AD); | 
(24) 


Ba=Lp Pr (Kofıs) -[BP]: | 


Die in (24) an zweiter Stelle rechter Hand stehenden 
Ausdrücke geben somit die ‚Korrektionsbeträge [xa] 


— [x'o’], bzw. [ßß] — [3’#’] für die Gewichtgrezi-. 


Bıökan an, die durch dıe Glieder höherer Ordnung 
verursacht werden. 
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PasaE 
Im Falle des Vorwärtseinschneidens — um einen 
markanten Fall Rn: — hat man 


I SENSE 


fÜ = are tg 


—.% 


worin &, y die Koordinaten des Neupunktes, 2 Yi 
die der Festpunkte sind. L; ist der auf (x;, y;) beob- 
achtete Richtungswinkel 9; nach dem Neupunkt 
(®, y). Die Ricktungskoeffizienten sind dann ?) 


of ,„,Sino; 
= Eu —o0 we | : 
Ar "x . (26), 
au 
Kr Ole , 
worin vr; A 
s=N aan, are 
„180% - 3600 “AT: 
= — (in Sek.) | | 
Durch Differentiation von (26) ergibt sich ®) ® 
„sin S 
fa = + DICEeE DE 
en sin = F 
Im = BI Zu Dual er mar: 
(28). 


q [23 
In + wu? («- ©) 
„, C08 2 - 
& go: 
Damit erhält man aus (24) und (27) die folgenden 
Schlußformeln: 


[x] =[%’ DE ne A-—[a’0’]AD) . . (29), 


[BBI= [BR] + ZH A-I[B'PJAD) . . (80), 


worin: 


— 


ee 
ß RaSESun 
AD Auol— Lac) +21ah] |»: er “2 @]| 


. VE HR 

Wir wollen nun. (29), (30), (31) auf das in O. Kg - 
gerts Handbuch der Vermessungskunde angegebene 
Beispiel ?) des Vorwärtseinschneidens aus drei Fest- 
punkten anwenden und entnehmen dort: Y 

[a a/—= 74, [bb] = 86, [ab] = — 2, so daß D= 6360, 
all: 14 = 0, 01351; BRI ==”, 01163, 

2 — 610 14°; 9, = 16° 42°; 9, = 316° 407, 
=+ 10”; ne 13; werTr 

Nach (28) ergeben sich der Reihe nach 

für dv feu: + 0,00212; — 0,00149; — 0,0026, 

für 9 fey: — 0,00135; — 0,00225; + 0,00013. Rz Be Ks ww 
Damit wird: Z=— 0,0153; AD= — 0,0048. vi 

Hiermit erhält man für die Geichtefiinten % ET 


Seitherige 
Berechnung: [%%]= 0,01 851 

Korrektions- ER 
glied:  [aaJ—1Xr/a7 

verbesserte EEE TOT Jr 
Werte: [aa] 001351 | BB 0, $, 


Die praktische Erprobung zeigt eindeutig, ‚der‘ 
Einwand von Piotrowski bei dem 
handelten Fall des Vorwärtseinschneid 
los ist, Die Er € 
einer Größenordnung, ae 
wichtsreziproken keineswegs 

Vgl. W. Tora und 
ma a ‚1, 8,8 
verfahren ermittelten Ent erı 

1936), 8.59 


AU 


er en, 
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F 
Bei regelmäßigen, vollständig symmetrischen Be-1 

stimmungsfiguren fällt die Bedingung (13) wegen (28) 
mit (10), bzw. mit (5), die automatisch nach einer 
Ausgleichung erfüllt sind, zusammen. Die vollkom- 
mene Symmetrie hat also ein Verschwinden jener Kor- 
rektionsglieder für die Gewichtsreziproken zur Folge. 
Sicherlich darf man annehmen, daß bei nur drei Be- 
stimmungsstrahlen in der Praxis die Symmetrie am 
‚ehesten gestört wird, so daß bei mehr Bestimmungs- 
strahlen wahrscheinlich auch keine größeren Korrek- 
tionsbeträge als die oben berechneten zu erwarten sein 
dürften. Andererseits wachsen jene Korrekturbeträge 
aber auch mit der Messungsungenauigkeit, d.h. mit 
- den v-Werten, an. Diese sind nun im obigen Beispiel 
durchaus nicht klein, so daß auch aus diesem Grunde 
keine wesentlich höheren Korrekturwerte zubefürchten 
sind. Aber selbst wenn diese bei abnorm asymme- 
trischen Figuren wider Erwarten einmal auf den 10- 
bis 100fachen Wert anwachsen sollten, sind sie bei 
trigonometrischen Einschaltungen so geringfügig, daß 
sie außer acht zu lassen sind, und es daher bei der 
seither geübten Berechnungsweise der Gewichtskoeffi- 
zienten verbleiben kann. 


Was übrigens die theoretische Strenge der Ableitung 
von Piotrowski betrifft, so darf nicht übersehen 
. werden, daß (7) im Grunde einer Entwicklung nach 
dem Satz von Taylor entspricht, die nur auf die 
linearen Glieder beschränkt ist. 


Frankfurt/M. H. Wolf. 


Beweis der Vektorformel 
axbxo)=la:-db—l(a-b)ec 
Zum Beweis genügt es, anzunehmen, daß b, c linear 
unabhängig sind. Klar ist, daß eine Beziehung der Form 
gilt 
ax(x a=Ab-+uc SSH); 
Wir betrachten jetzt die folgenden Vektoren: 
Ba. DAX. CH .C. XUD7Xc)E 
cr 
UI SE TexX (6, X c)] 
ss ee 
& 
von denen je drei aufeinanderfolgende ein Rechtsdrei- 
bein bilden; das zweite beginnend mit c ist überdies 
orthogonal, und daher auch das dritte, woraus folgt 
&= ocmit oe >0, und genauer &= (b X c)?:c. Fer- 
ner ist 
DieRoE URL TER (DE X IC) (BEXIC)2. 
Multiplizieren wir also (1) links und rechts skalar mit 
c*, so ergibt sich R 
e*H-fax(bxoj=ick-b 
m —— ——— u 


U HDIRS ERDE ART BR C)? 

(ac) (BE x ec)? 

d.h.}= (a: c). Vertauschung von b und c ergibt ana- 

log u= — (ab). 
Berlin-Charlottenburg. 


E. Mohr. 


BUCHBESPRECHUNGEN 


B.L. van der Waerden (Prof. a. d. Univ. Amster- 
dam), Moderne Algebra. Unter Benutzung 
der Vorlesungen von E.ArtinundE.Noether. 
Erster Teil. (Die Grundlehren der Mathematischen 
Wissenschaften in Einzeldarstellungen mit besonderer 
Berücksichtigung der Anwendungsgebiete, Band 
XXXILL.) 3. Aufl. VIII + 292 S. Berlin- Göttingen- 
Heidelberg 1950. Springer-Verlag. Preis brosch. 24,— 
DM, geb. 27,— DM. 

Das Buch bietet eine glückliche Synthese der klassi- 

. schen Algebra mit den modernen Begriffen der Grup- 
pen-, Körper-, Bewertungs- und Idealtheorie sowie 
der Theorie der hyperkomplexen Zahlen. Die Dar- 
stellung ist klar und für die Einführung mit großem 

- Geschick abgefaßt. Die dreiersten Kapitel dienen der 
; - Vorbereitung und enthalten die Grundbegriffe über 
3 Mengen, Gruppen, Ringe, Körper und Ideale. Dem 
schließt sich wesentlich die Theorie der kommutativen 
Körper an; sie wird bis zum Begriff der Bewertung 
eines Körpers geführt. Zahlreich eingestreute Auf- 
gaben unterstützen die Probe auf das Verständnis des 
Textes. £ . 


Dresden. 


Dr.H, Hasse (o. Prof. a. d. Humbold-Univ. Berlin), 
_ VorlesungenüberZahlentheorie. (Die 
Grundlehren der Mathematischen Wissenschaften in 
 Einzeldarstellungen mit besonderer Berücksichtigung 
der Anwendungsgebiete, Band LIX). XII + 474 S. 
mit 28 Abb. Berlin-Göttingen-Heidelberg 1950. 
Springer-Verlag. Preis 42,— DM, geb. 45,— DM. 
...% Das. Buch folgt der im Vorjahre erschienenen 
„Zahlentheorie‘‘ des Verfassers. Die im Vorwort er- 
örterte Frage nach der Existenzberechtigung dieses 
- neuen Werkes möchte ich nach einer anderen Richtung 
ergänzen: Die ‚Vorlesungen‘ füllen eine wesentliche 
Lücke in der vorhandenen Lehrbuchliteratur aus, 
insofern sie in die „elementare“ Zahlentheorie in 
icht faßlicher Form die modernen Begriffe der 
ruppe und des Ringes einführen (Abschn. I u. II). 
den vielfach erhebliche Vereinfachungen 


pr m 


Draeger. 


a EHE 


der Beweise erzielt. Der III. Abschnitt bringt eine 
geradezu glänzende Analyse des berühmten Dirichlet- 
schen Primzahlsatzes, und an ihn schließt sich in 
engem innerem Zusammenhang eine Theorie der 
quadratischen Zahlkörper an, deren Fassung wieder- 
um eine moderne Zusammenschau der Probleme er- 
möglicht. An Vorkenntnissen in mathematisch- 
algebraischer Hinsicht werden wesentlich nur die 
Grundlagen der Algebra vorausgesetzt. Auf die vielen 
Einzelvorzüge des Werkes einzugehen, ist in diesem 
Rahmen unmöglich. Dem Interessenten an dieser 
reinsten Mathematik kann nur wärmstens empfohlen 
werden: Nimm und lies! 


Dresden. Draeger. 


Dr.-Ing. Lothar Cremer (Dozent a. d. Univ. Mün- 
chen), Die wissenschaftlichen Grund- 
lagen der Raumakustik, Band III. Wel- 
lentheoretische Raumakustik. XIV 
+ 3558. mit 87 Abb. Leipzig 1950. S. Hirzel-Ver- 
lag. Preis geb. 21,50 DM. 

Es handelt sich um den III. Teil eines zusammen- 
fassenden Werkes des Verfassers über ‚‚die wissen- 
schaftlichen Grundlagen der Raumakustik“. Wäh- 
rend der I. Teil von den anschaulichen Darstellungen 
der Experimentalphysik über die Wellenausbreitung 
ausgeht und die Einflüsse der geometrischen Raum- 
form aufzeigt (Echo, Schallkonzentrationen und ihre 
psychologische Wirkung) und der II. (noch nicht er- 
schienene) Teil aus statistischen Betrachtungen die 
Nachhallerscheinungen und Materialeinflüsse ab- 
leiten wird, geht der III. Teil von einer strengen 
Formulierung der Schallwellentheorie aus. Er dürfte 
damit den Mathematiker am meisten interessieren. 
Zunächst werden die Schallfeldgleichungen abgeleitet. 
Von dieser Grundlage ausgehend werden dann die 
raumakustischen Probleme behandelt, die ohne höhere 
Mathematik nicht dargestellt werden können: mit- 
schwingende Wände und Raumresonanzen mit ihrer 
Anwendung in einfacher zu übersehenden Räumen, 
im Rohr, im Trichter, im Rechteckraum, in der 
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Flachkammer zur Messung der Schallschluckung, 
auch bei schrägem Einfall und durch parallel zur 
Ausbreitung laufende Schluckflächen. Zum Schluß 
wird die Schallabsorption durch innere Reibung und 
Wärmeleitfähigkeit der Luft sowie durch Anregung 
von Molekülschwingungen untersucht. Durch seine 
allgemeine Grundlage wird sich das Buch weit über 
den Rahmen der speziell interessierten Fachwelt, 
besonders auch bei den technisch interessierten Phy- 
sikern Freunde erwerben. 


Dresden. 


Reichardt. 


C- Truesdell, An essaytowarda unified 
theory of special functions -based 


upon the functional equation Ira) 


—= F(z,& +1). (Annals of Mathematics Studies 
Number 18.) IV + 182 S. Princeton 1948. Princeton 
University Press. Preis 3,— $. 


Die vertrauten speziellen Funktionen f(y, &) der an- 
gewandten Mathematik besitzen eine Fülle funktions- 
eigener und verknüpfender Beziehungen, die durch 
ihre Gestalt auf einen tieferen Zusammenhang hin- 
weisen. Seine Quelle erkennt der Verfasser in der 
gemeinsamen Differenzen-Differentialgleichung 


0 
RA BEFUND. 


Ist weiter die Koeffizientenbedingung 


0) 0 
20 Alu, a+ 1A a)4+ 5, InB (ya) |- 0 


erfüllt, dann lassen sich, wie Trues dell zeigt, aus 
den Koeffizienten die Funktionen p(y,&) und 
q(y) = z so bestimmen, daß p (y, a) f(y, a) = F (2, &) 
der im Buchtitel angegebenen, schon vielfach unter- 
suchten koeffizientenfreien Funktionalgleichung ge- 
nügt. Von ihren Lösungen ausgehend gelangt man 
je nach Wahl des Anfangswertes F(z,a) und der 
Koeffizienten durch Rücktransformation zu den 
Funktionen f(y,&) samt ihren Eigenschaften und 
gegenseitigen Beziehungen. Hierdurch gewinnt die 
Theorie der speziellen Funktionen eine bewunderns- 
werte Einheitlichkeit und Durchsichtigkeit; daß sie 
auch bisher unbekannte Beziehungen aufdeckt, be- 
legt der Verfasser mit eindrucksvollen Beispielen. 


Dresden. P,E. Böhmer. 

Dr. R. Jaeckel (apl. Professor a. d. Universität 
Bonn) unter Mitarbeit von Dr.H. Schwarz und 

„Dr. E Sehüller, Kleinste Drucke, ihre 

Messung und Erzeugung.- X+ 302 S. m. 
301 Abb. Berlin, München 1950. Springer-Verlag und 
J. F. Bergmann. Preis geh. 39,60 DM. 

Der Verfasser gliedert seine Darstellung in sechs 
Hauptteile: Kinetische Gastheorie; Vakuum-Meß- 
instrumente; Pumpen; Erzeugung, Erhöhung und 
Aufrechterhaltung des Vakuums ohne Pumpen; 
Vakuum-Verbindungen und -Leitungen; Vakuum- 
zubehör. 

Der .erste Teil gibt eine Zusammenstellung von 
Formeln der kinetischen Theorie der Gase, die in den 
folgenden Abschnitten gebraucht werden. Die näch- 
sten drei Abschnitte über Vakuum-Messung- und 
-Erzeugung bilden den Kern des Buches und haben 
auch im Umfang das größte Gewicht. Hier sind die 
Geräte und Verfahren beschrieben und ihre Wirkungs- 
weise durch zum Teil ausführliche Rechnungen im 
einzelnen theoretisch untersucht; Die beiden letzten 
Abschnitte geben praktische Hinweise für den Auf- 
bau von Vakuumapparaturen. Sie bilden eine Er- 
gänzung zu den vorherigen Kapiteln und geben dem 
Werk den Charakter eines allgemeinen Hilfsbuches 
der Hochvakuumtechnik, da der im Untertitel des 


Buches erläuterte Umfang der Darstellung in diesen 
Abschnitten überschritten wird. Besonders hervor- 
zuheben ist, daß zahlreiche, auch ausländische Neu- 
entwicklungen dargestellt worden sind. Ein viel- 
seitiges Literaturverzeichnis ist beigegeben. Das 
Erscheinen dieses Buches ist sehr zu begrüßen, um so 
mehr, als zusammenfassende deutsche Darstellungen 
über die Hochvakuumtechnik längere Zeit vergriffen 
waren, ; 


Dresden. Recknagel. 


! 


Mathematisch-physikalische Semesterberichte zur 
Pflege des Zusammenhanges von Schule und Uni- 
versität, herausgeg. von H. Behnke, Münster und 
W. Lietzmann, Göttingen in Verbindung mit 
dem math. Forschungsinstitut in Oberwolfach. Bd.1, 
318 S. Göttingen 1949/50. Vandenhoeck und Rup- 
recht. Preis 19,60 DM (Subskriptionspreis 17,— DM). 

Zum Geleit für das Neuerscheinen dieser Semester- 
berehte sagt H. Behnke: ,‚‚,Den Lehrern der 
Mathematik und der verwandten Wissenschaften an . 
den höheren Schulen soll unser Bemühen um diese 


Blätter in erster Linie gewidmet sein“. Der „Gym- 
nasiallehrer als Mittler zwischen der Forschung und ; 
dem Leben der geistig und fachlich interessierten € 
Schichten unseres Volkes‘ soll durch die Berichte j 
über den Fortgang der Forschung auf dem laufenden i 
gehalten werden. 


Die Wiederaufnahme dieser hier durch Berichte 
und dort durch Ferienkurse geübten Tradition kann 
nicht warm genug begrüßt werden. Was immer schon | 
galt, ist heute angesichts der Lage und auch der i 
Forderungen der Gegenwart womöglich in erhöhtem 
Maße wahr. Der Lehrer, der es mit seinem Berufe L 
ernst nimmt, kann um der letzten und höchsten 
Unterrichtserfolge willen in seiner Fachwissenschaft 
gar nicht weit und tief genug verankert sein. Daraus 7 
folgt die Notwendigkeit der entsprechenden Studien- 
pläne ebenso wie diejenige seines wissenschaftlichen 
Weiterarbeitens auch nach dem Studium, Ihm dafür 
fortlaufend geeignete Handreichungen und Anregungen 


zu bieten, ist darum ebenfalls ein Erfordeınis. h 

Entspricht der Inhalt der vorliegenden Semester- 
berichte dem? Diese Frage ist unbedingt zu bejahen. 
und den Berichten daher weiteste Verbreitung zu 
wünschen. Der Inhalt ist denkbar vielseitig, und den 5 
meisten Aufsätzen sind weiterführende Literatur- 1 


angaben hinzugefügt. Der Inhalt erstreckt sich von 
der Darstellung mengentheoretischer Probleme bis R 
zu einem Aufsatz über mathematische Geräte und ° 1 
umfaßt zugleich mehrere Artikel über die moderne | 
Physik. Beiträge über F. Klein, Toeplitz, Bavink und | 
„Goethe und die Mathematik“ ergänzen den sach- j 
lichen Inhalt nach der persönlichen Seite — für den 
Lehrer von besonderem Werte, zumal wenn dabei, | 
wie in dem Aufsatze über*Goethe und die Mathematik, 
zugleich Einsichten in das Wesen der Mathematik 
und in mathematische Strukturprobleme vermittelt. 
werden. 

Die Darstellung ist zumeist konzentriert knapp. und 
Mitunter — ich nenne als Beispiel den Beitrag von i 
Hermes über die Funktionenanalysis — wäre eine 


Fingerzeige. 
Dresden. 
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C. Zwikker (Technical Director Light Division Phi- 
lips Eindhoven), Advanced Plane Geometry. 
XI -+299 S.m. 273 Abb. Amsterdam 1950. North- 
Holland Publishing Company. 


Verf, bietet einen Gesamtüberblick über alle Ge- 
biete der Geometrie, die irgendwie für die Anwendung 
in der Technik wichtig werden können, also analytische, 
projektive, algebraische, Differentialgeometrie, Fra- 
gen der Potentialtheorie, der Kinematik. Im Ein- 
zelnen bstrachtet er aus der projektiven Geometrie die 
Eigenschaften von Punktreihe und Strahlenbüschel, 
Dualität, die Schnittpunktsätze; aus der Dreieckslehre 
den Fragenkreis des Feuerbachschen Kreises; vom 
Kreise die Kreisverwandtschaften und die projektive 
Erzeugung des Kreises, jedoch keine Büscheltheorie. 
Aus der algebraischen Geometrie werden Polareigen- 
schaften und -verwandtschaft, sowie eine Theorie der 
Kurven 3, Ordnung mit Hilfe der elliptischen Funk- 
tionen behandelt. In der Theorie der Kegelschnitte 
beschränkt sich Verf. auf die metrischen und affinen 
Eigenschaften des einzelnen Kegelschnittes und geht 
dabei ausführlich auf kinematische Eigenschaften und 
die Keplerbewegung ein. Aus der Differentialgeometrie 
behandelt er Evolventen und Evoluten, Enveloppen, 
Fußpunktskurven, Kaustiken. Orthogonale Trajek- 
torien einer Kurvenschar behandelt er hauptsächlich 
im Hinblick und in Anwandung auf Potentialtheorie. 
Ein Kapitel bringt unstetige und geknickte Kurven, 
in erster Linie als Grenzkurven in einer Kurvenschar, 


‚in diesem Zusammenhang die Schwarzsche Poly- 


gonabbildung. Im einzelnen geht er dann noch auf die 
Spiralen, Lemniskaten, Radlinien, Konchoiden ein. 
Ein Kapitel über Zahnprofile von Zahnrädern bildet 
den Schluß. In jedem Kapitel wird die Kraft der er- 
arbeiteten Methoden an schönen Anwendungen auf 
die Elektrotechnik, Kinematik u. ä. erprobt, die auch 
mathematisch äußerst reizvoll sind. 

erraschend ist die Darstellung. Verf. faßt die 
Ebene durchweg als Wertebereich einer komplexen 
Veränderlichen auf und stellt jede Kurve durch eine 
komplexe Funktion eines reellen Parameters dar. Die 
komplexen Zahlen bieten alle Vorzüge der Vektorrech- 
nung, sind aber schmiegsamer. Freilich muß Verf. 
dann in der Projektiven Geometrie auf die Benutzung 
homogener Koordinaten verzichten. 

In seinen Voraussetzungen und Anforderungen will 
das Buch eine mittlere Lage einnehmen. Die Elemente 
der Infintesimalrechnung und der. analytischen Geo- 
metrie, wie sie bei uns etwa in den beiden ersten Stu- 
diensemestern gelehrt werden, werden vorausgesetzt, 
um für eine Fülle feinerer Eigenschaften der Figuren 
Platz zu gewinnen. Andrerseits verbietet schon der 
weite Rahmen, in den einzelnen Gebieten zu weit den 
feinsten Verästelungen zu folgen; es wird dann immer 
auf Spezialwerke verwiesen. Das Buch wendet sich 


. wohl in erster Linie an den angehenden Ingenieur, dem 


die mathematische Denk- und Ausdrucksweise bereits 
etwas vertraut ist. Dieser wird darin wohl alles fin- 


‚den, was ihm die Geometrie für seine Arbeit zu bieten 


vermag. Umgekehrt lernt der Mathematiker aus dem 
Buch eine Fülle von Beziehungen seiner Wissenschaft 
zur Technik. Eine verständliche Darstellungsweise 


_ und eine große Zahl von sehr anschaulichen Figuren 


machen das Studium des Buches zu einer Freude. 
Dresden, Ott-Heinrich Keller. 


A. Schmidt (Prof. a.d. Univ. Göttingen), Mathe- 
matische Grundlagenforscehung. En- 
zyklopädie der mathematischen Wissenschaften mit 
Einschluß ihrer Anwendung“. Band I,1, Heft1, Teil II, 


488. Leipzig 1950. B. G. Teubner Verlagsgesellschatt. 


Preis geh. 5,60 DM. 

Auf dem internationalen Mathematikerkongreß in 
Zürich 1932 wurde beschlossen, die Enzyklopädie der 
Mathematischen Wissenschaften neu herauszugeben 
aber ohne die Abschnitte aus der angewandten Mathe- 


 matik. Diese sollten später in einem besonderen Werk 
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zusammengefaßt werden. Zuerst wurde unter der Be- 
treuung der Professoren Hasse und Hecke’die 
Neubearbeitung des ersten am meisten veralteten 
Bandes ‚‚Algebra und Zahlentheorie‘ in Angriff ge- 
nommen. Leider waren die Zeitverhältnisse dem Un- 
ternehmen nicht günstig. Erst jetzt erscheint als 
erstes Heft dieses Bandes der Abschnitt über ‚‚Mathe- 
matische Grundlagenforschung“, der von A. Schmidt 
bereits im Jahre 1939 geschrieben wurde, der aber 
durch meist kenntlich gemachte Einschaltungen die 
Entwicklung in den letzten zehn Jahren, auch die im 
Auslande berücksichtigt. In diesem ersten Artikel 
scheint es dem Verfasser gelungen zu sein, die doppelte, 
nicht leichte Aufgabe der Enzyklopädie zu erfüllen, 
nämlich: dem Fachmann eine zuverlässige Übersicht 
über Inhalt und besonders über dieneuere Literatur des 
Gebietes zu geben und außerdem denen, die nur über 
eine gewisse mathematische Allgemeinbildung ver- 
fügen eine überschlägige Kenntnis des Gebietes zu ver- 
mitteln. Dadurch, daß den Herausgebern auch dieses 
Mal ein internationaler Mitarbeiterstab für die ein- 
zelnen Sachgebiete zur Verfügung steht, ist höchste 
Zuverlässigkeit auch in dieser neuen Auflage gewähr- 
leistet. Man kann nur wünschen, daß das Werk dieses 
Mal in kürzerer Zeit abgeschlossen werden möchte als 
seiner Zeit die erste Auflage. 


Dresden. Willers. 


W, Blaschke, Einführungin die Diffe- 
rentialgeometrie. (Grundlehren der Math. 
Wissenschaften in Einzeldarstellungen mit beson- 
derer Berücksichtigung der Anwendungsgebiete. Bd. 
LVIIL) VII + 146 S. mit 57 Abb. Berlin-Göttingen- 
Heidelberg 1950. Springer-Verlag. Preis brosch. 
16,— DM., geb, 18,60 DM. 

Verf. gibt auf knappstem Raum eine umfassende 
Darstellung der metrischen Differentialgeometrie des 
Raumes. Er beginnt nieht mit der Kurventheorie, 
sondern behandelt die Kurven von vornherein als Trä- 
ger von Streifen. Krümmung, Windung, Frenetsche 
Formeln erscheinen als Spezialfälle von Streifengrößen 
und -formeln. In der Flächentheorie stellt er die in- 
neren Eigenschaften und die geodätischen Linien voran 
und bringt dann erst die äußere Flächenlehre. — Als 
Hilfsmittel benutzt er nicht, wie man sonst gewohnt 
ist quadratische, sondern im Anschluß an Cartan 
lineare Pfaffsche 'Differentialformen. Dadurch und 
durch ausgiebige Benutzung symbolischer Schreib- 
weisen gewinnen die Formeln eine überraschende 
Kürze und Durchsichtigkeit, selbst solche wie die für 
die Gaußsche Krümmung, die von Gauß-Codazzi 
usw. 

So ist das Buch für den reifen Mathematiker, der etwa 
den Inhalt des 1. Bandes der berühmten Blaschke- 
schen Werke über Differentialgeometrie beherrscht, 
ein Wegweiser zu einer wesentlichen Vertiefung seiner 
Einsichten und ein Genuß. 

Noch eine Frage: Auf S. 104 wird der ‚hervorragende 
jüdise Zahlentheoretiker und Geometer Min- 
kowski‘ erwähnt. Muß manin Deutschland wirklich 
schon wieder nach der Abstammung eines Gelehrten 
fragen ? 

Dresden. Ott-Heinrich Keller. 

6. Bol, Projektive Differentialgeo- 
metriel. (Math. Lehrbücher, herausgeg. vom Math. 
Forschungsinstitut Oberwolfach unter Leitung von 
Prof. Dr. W. Süss Bd. IV). VII + 365 S. m. 15 Abb. 
Göttingen 1950. Vandenhoeck u. Ruprecht. Preis 
brosch. 17,80 DM, geb. 20,— DM. 

Die projektive Ditferentialgeometrie behandelt die- 
jenigen differentialgeometrischen Eigenschaften der 
Figuren, die gegenüber projektiven Abbildungen in- 
variant sind. Sie wurde hauptsächlich in den letzten 
50 Jahren in Amerika, Italien, Japan gepflegt. Das 
Buch ist das erste Lehrbuch in deutscher Sprache. 
Dieser I. Bd. bringt die projektive Differentialgeo- 
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metrie der ebenen Kurven, dann, nach einer Einfüh- 
rung in die räumliche projektive Geometrie, die Li- 
niengeometrie und in die Lehre von den Asymptoten- 
linien die Theorie der Raumkurven und Flächenstrei- 
fen. Der II. Bd. soll dann die Flächentheorie behan- 
deln. Mehrdimensionale Fragen liegen bis auf ge- 
legentliche Hinweise außerhalb des Blickfeldes. Beim 
Aufbau umgeht Verf. die Normierungsfrage, die sich 
offenbar nicht ungekünstelt behandeln läßt. Dafür 
richtet er seine Definitionen so ein, daß sie sich bei 
Parameter-Transformationen möglichst einfach ver- 
halten. In einem Anhang sind die benötigten Sätze 
und Formeln aus anderen Gebieten, z. B. Differential- 
gleichungen, zusammengestellt. Einem vertieften 
Studium werden durch ein umfassendes Literatur- 
verzeichnis der Arbeiten von 1931—1950 (seit dem 
Lehrbuch von Fubini-Cech) die Wege geebnet. 

Die Darstellung ist so gehalten, daß ein Student der 
Mathematik der mittleren Semester ohne große An- 
strengung folgen kann. Kenntnisse der klassischen 
Differentialgeometrie und der projektiven Geometrie 
werden nicht vorausgesetzt, wohl aber ein gewisses 
Vertrautsein mit und ein’Geschmack an deren Frage- 
stellungen. Es ist zu hoffen, daß das Buch die Kennt- 
nis seines Gegenstandes auch in Deutschland allgemein 
verbreiten und für seine eigentümlichen Schönheiten 
Freunde und Mitarbeiter werben wird. 


Dresden. Ott-Heinrich Keller. 


Dr. Wilhelm Maack (apl. Prof.a.d. Univ. Hamburg), . 


Differential-und Integralrechnung. 
Bücher der Mathematik und Naturwissenschaften, 
herausgeg. v. Dr. Henry Poltz). 2358. m. 17 Abb. 
Wolfenbüttel und Hannover 1949. Wolfenbütteler 
Verlagsanstalt G. m. b. H. 


Der Verf. beabsichtigt mit diesem Buch, ‚‚den Leser 
in großen Zügen mit den Methoden und Fragestel- 
lungen der Differential- und Integralrechnung bekannt 
zu machen. Es wird genau so viel geboten, wie un- 
erläßlich scheint, dem Lernenden Grund und Boden 
unter die Füße zu geben“. Man darf sagen, daß.das 
hierdurch gesteckte Ziel durchaus erreicht wurde. Der 
Leser wird in lebendiger Weise zum abstrakten ana- 
lytischen Denken geführt, ohne ihn etwa durch den 
Zwang, alles in möglichst weitgehender Allgemeinheit 
zu betrachten, zu ermüden,. — Der erste Teil behandelt 
die Funktionen einer Variablen. Es wird von der 
reellen Zahi als etwas Gegebenem ausgegangen; nach 
Einführung des Funktionsbegriffs und der Behandlung 
des Grenzwertes von Zahlenfolgen und Funktionen 
folgen Stetigkeit und differenzierbare Funktionen. 


Hieran schließt sich, nach einem Exkurs über unend-, 
liche Reihen das Riemannsche bestimmte Integral 
“und der Fundamentalsatz der Differential- und Inte- 


gralrechnung, wobei man sich bei der Berechnung un- 
bestimmter Integrale auf solche rationaler Funktionen 
beschränkt. Der zweite Teil behandelt Differentiäl- 
und Integralreehnung von Funktionen zweier Ver- 
änderlichen bis zum Gaußschen Integralsatz unter 
den üblichen vereinfachenden Annahmen. — Eine 
Reihe netter Kleinigkeiten in der Darstellung und den 
Beweisen wird auch dem fortgeschrittenen Leser Neues 
bieten. 


Dresden. ‚Maruhn. 

Wilhelm Maack (Prof. a.d. Univ. Hamburg), Fast - 
periodische Funktionen. (Die Grund- 
lehren der Mathematischen Wissenschaften in Einzel- 
darstellung mit bes. Berücksichtigung der Anwen- 
dungsgebiete, Bd. LXI.) VIII + 214 S. Berlin/Göt- 
tingen/Heidelberg 1950. Springer-Verlag. Preis brosch. 
21,60 DM., geb. 24,60 DM. 

In den Jahren 24 bis 26 schuf H. Bohr seine 
Theorie der fastperiodischen Funktionen. Es sind 
dies solche reelle oder komplexe Funktionen der 
reellen Variabeln x, die erstens einmal zu jedem vor- 
geschriebenen e>0 eine Fastperiode 7 besitzen, für 


/ 
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die also und zwar für jedes x die Ungleichheit 
Ka@+rD)—f(e))<e gilt, für dieesaber zweitens, weil 
das sonst zu vage wäre, zu jedem &>0 eine Zahl 
L(e) gibt, so daß in jedem Intervall der Länge Z eine 
Fastperiode enthalten ist. Diese Funktionen lassen 
sich durch endliche trigonometrische Reihen beliebig 
genau gleichmäßig approximieren. Die. Bedeutung 
der Theorie für die Praxis liegt damit auf der Hand. 
Es ist daher dem Verf. zu danken, daß er die Theorie 
dieser Funktionen, die er jetzt eigentlich fastperiodisch 
nennt, besonders und unabhängig im vierten Kapitel 
dargestellt hat. Die Theorie ist nämlich seit Bohr 
erheblich erweitert worden. Ihre Grundlage ist der 
Begriff der Gruppe in der heutigen Fassung als einer 
Menge von Operationen, für die eine assoziative Multi- h 
plikation existiert. Daherbeginnt das vorliegende Buch 
mit den Darstellungen endlicher Gruppen, um dann im 
zweiten Kapitel die abstrakte Theorie der fast perio- 
dischen Funktionen auf Gruppen anzuschließen. Von 
dieser allgemeinen Theorie aus ist dann der Übergang 
zu Sonderfällen verhältnismäßig einfach, so zu den 
periodischen Funktionen auf der ebenen Drehungs- | 
gruppe, zu den eigentlichen fastperiodischen Funk- 
tionen auf der Gruppe der reellen Zahlen und endlich 
zu den Kugelfunktionen auf der räumlichen Drehungs- ’ 
gruppe. Dieses siebente und letzte Kapitel bedarf d 
im vorhergehenden Kapitel noch einer tiefgründigen 
Untersuchung über kompakte Gruppen. Im Anhang | 
finden sich Literaturhinweise mit einer historischen 
Darstellung der Entwicklung der Theorie. Dieser An- 
hang wie auch die Einleitungen zu den einzelnen 

Kapiteln sind in einer besonders angenehm. lesbaren 

Form geschrieben, was das Verständnis der recht 
schwierigen Theorie erheblich erleichtert. Ref. er- 

freut es zu sehen, daß seine Untersuchung über eine 

Basis aller Zahlen usw. schon so historisch geworden 

ist, daß man sie ohne Nennung seines Namens vor- 

tragen kann, was $8. 87 und 88 geschieht. 


Landshut. Hamel. 


W., Lietzmann (Prof. a. d. Univ. Göttingen), ‚Wo 
steckt der Fehler?“ Mathematische 
Trugschlüsse und Warnzeichen. 183 S, 
mit 120 Abb. Leipzig 1950.: B. G. Teubner Verlags- _ 
gesellschaft. Preis geb. 5,40 DM. 


Das kleine Bändchen erschien zum ersten Mal im 
Jahre 1913 in der von Lietzmann und Wit- 
ting herausgegebenen Mathematischen Bibliothek 
in einem Umfang von 57 Seiten. Von Auflage zu Auf- Fe 
lage nahm der Umfang, besonders dank der Mitarbeit 
der Leser, zu. In der vorliegenden ist den beiden Ab- i [ 
schnitten ‚‚Täuschungen und Fehlschlüsse‘“ und 
‚„‚Trugschlüsse“ ein neuer ‚Warnzeichen aus der Ana- 
lysis des Unendlichen‘‘ hinzugefügt worden. Der In- 
halt des Bändchens ist Unterhaltungsmathematik im 
besten Sinne des Wortes. Natürlich sind nicht alle _ 
Beispiele gleichwertig; aber es findet sich viel vor- 
treffliches darunter. Durch Verwendung manches 
dieser Beispiele wird man Unterricht oder Vorlesung _ 
anregender gestalten können. Die sorgfältige Durch- 
arbeitung ist eine ausgezeichnete Schulung des logi- 
schen und mathematischen Denkens. Möchte sin 
recht viele Hände kommen. PER 


Dresden. Wi 1 le TS x Dr 
Jahresbericht der Deutschen Mathematiker-Ver- 
einigung. Herausgegeben, von Prof. Dr. H. Brandt 
(Halle) und Prof. Dr. E. Kamke (Tübingen). Bd.54 
Heft 1. 54 u. 18 S. Bielefeld 1950. Verlag für Wissen- 
schaft und Fachbuch. Preis: 12,75 DM., für Mit- 
glieder der DMV. 850DM. HE 
Als eine der letzten deutschen mathematischen V, or R 
kriegszeitschriften hat jetzt nach dem Neuaufbau < 
Deutschen Mathematiker-Vereinigung auch der Ji 
berichtder DMV. sein Erscheinen wieder aufgenom 
Er wird jetzt von Prof. Dr. Brandt-Halle u 
Dr. Kamke-Tübingen herausgegeben und s 
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alten Weise fortgeführt werden. Er hat daher wie 
früher zwei Teile: ‚Berichte und Abhandlungen‘ und 
‚Kleinere Mitteilungen“. Bd. 54 wird in zwei Heften 
von insgesamt neun Bogen erscheinen. Das vor- 
liegende erste Heft bringt im ersten Teil eine Gedenk- 
tafel für die seit 1933 verstorbenen deutschen Mathe- 
matiker. Es folgen vier Abhandlungen. Am Schluß 
des zweiten Teiles ist ein Verzeichnis der augenblick- 
lichen Mitglieder gegeben. Hoffen wir, daß die Bände 
bald wieder im alten Umfang erscheinen können. 


Dresden. Willers. 


F. Hund (o. Professor an der Universität Jena), 
Einführungin die theoretische Phy- 
sik. IV. Band: Theorie der Wärme (Meyers 
kleine Handbücher, Bd. 54/55). 330 S. mit 46 Abb. 
Leipzig 1950. Bibliographisches Institut. Preis 
5,80 DM. V. Band: Atom- und Quanten- 
theorie (Meyers kleine Handbücher, Bd. 56/57). 
314 S. mit 48 Abb. Leipzig 1950. Bibliographisches 
Institut. Preis 5,380 DM. | 


Die beiden letzten Bändehen der theoretischen 
Physik von Hund stehen insofern in engerer Be- 
ziehung, als der vierte Band nicht nur die phänomeno- 
logische Thermodynamik enthält, sondern auch schon 
die statistische Theorie der Wärme in den Grundzügen 
behandelt. Ferner wird die Theorie der Wärme- 
strahlung gebracht, wobei die Quantenhypothese ein- 
geführt werden muß. Im fünften Band wird hierauf 
zunächst noch einmal eingegangen, wo das Versagen 
der klassischen Theorie bei Fragen des Atombaues be- 
sprochen wird. Es folgt dann die ältere Quanten- 
theorie, die mit dem Korrespondenzprinzip arbeitet, 
und schließlich die moderne Quantenmechanik. Dem 
Charakter einer Einführung entsprechend sind manche 
Gebiete verhältnismäßig kurz dargestellt. Auffällig 
ist diesin den Abschnitten über die statistische Theorie 
der Wärme, die sehr wenig Anwendungen enthalten. 
Die Darstellung der Atomphysik ist ausgezeichnet 
durch vorbildliche Klarheit. Dieses wohldurchdaehte 
Bändchen zu studieren bereitet Genuß. 


Dresden, Receknagel 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr. J. Horn (ehem. o. Prof. a. d. Techn. Hoch- 
schule Darmstadt). Partielle Differentialglei- 
chungen. (Göschens Lehrbücherei Gr.I, Bd. 14.) 
Vierte, unveränderte Auflage. 228 S. mit 8 Abb. 
Berlin 1949. Walter de Gruyter & Co. Preis geb. 
14,— DM. | 


Suter-Traub, Die Methode der Festpunkte. 
Dritte, neubearbeitete Auflage, XIL + 216 S. mit 
232 Abb. und 7 Tafen. Berlin-Göttingen-Heidelberg 
1951. Springer-Verlag. Preis geb. 21,— DM. 


Dr. H. Küstner, Fünfsteilige Logarithmen der 
natürlichen Zahlen und der Winkelfunktionen bei 
dezimalgeteiltem Altgrad. Anhang physikalisch- 


n 


technische Tabellen. VI-+ 157 
1950. Volk- und Wissen-Verlag. 


S. Berlin-Leipzig 


Dr.-Ing. Hans Mahl und E. Gölz, Elektronen- 
mikroskopie (Die kleinen Studienbücher, Band 
58). 1798. mit 39 Abb. und 5 Tafeln. Leipzig 1951. 


. Bibliographisches Institut. Preis kart. 5,80 DM. 


Dr. A, P. Speiser, Entwurf eines elektro- 
nischen Rechengerätes unter besonderer 
Berücksichtigung der Erfordernis eines minimalen Ma- 
terialaufwandes bei gegebener mathematischer Lei- 
stungsfähigkeit. (Mitteilungen aus dem Institut für 
angew. Math. a. d. E. T.H. in Zürich). 548. mit 
4Abb. und 4 Tabellen. Zürich 1950. 


NACHRICHTEN 


Darmstadt: Dr.rer.nat. W. Schöbe, Prokurist 
bei der Allianz-Lebensversicherung A.G., in, Stuttgart, 
j erhielt die venia legendi für Mathematik. x 
f Dr.-Ing. habil. Theodor Zech, Studienassessor an 
der Odenwaldschule in Hambach wurde nach seiner 
Rückkehr aus der Kriegsgefangenenschaft erneut als 
Privatdozent zugelassen. 
Dr.-Ing. H. Buchholz, früher Dozent an der Techn. 
Hochschule Berlin-Charlottenburg und wissenschaft- 
' licher Mitarbeiter der AEG, jetzt wissenschaftlicher 
' Mitarbeiter am Post- und Fernmeldetechnischen Zen- 
tralamt Darmstadt, erneuerte seine venia legendi jetzt 
an der Techn. Hochschule Darmstadt und wurde kurz 
darauf zum nichtplanmäßigen Professor ernannt. 


) 


% 


Zu H. Holstein. Uber die äußere und innere 
laminarer 


Ei ‚vor kurzem in dieser Zeitschrift veröffentlichter 
tikel von H. H ale tein [1] enthält Einwände 
n einige der Ergebnisse, welche die Verfasser der 

enden "Zuschrift in ihren Arbeiten [2] und [3] 


E in » 
- « n 


Dr. rer. nat. H. Söhngen, früher wissenschaftlicher 
Mitarbeiter an der Deutschen Versuchsanstalt für Luft- 
fahrt Berlin, zuletzt in Frankreich in der Industrie- 
forschung tätig, erhielt die venia legendi für Mathe- 
matik. 


Der nichtplanmäßige Professor für technische Me- 
chanik, Dr.-Ing. H. Dietz ging unter Beibehaltung 
seiner venia legendi als wissenschaftlicher Mitarbeiter 
an die MAN Mainz-Gustavsburg über. 


Freiburg i.Br.: Dr. H. Görtler, o. Prof. für 
angew. Math. a. d. Univ. Freiburg wurde zum Mit- 
glied des Institute of Mathematical Statisties in Ann 
Arbor, Michigan USA ernannt. 


 ZUSCHRIFTEN AN DEN HERAUSGEBER 


veröffentlicht haben. Es scheint den Verfassern, daß 
diese Kritik im wesentlichen auf einem Mißverständnis 
beruht. Diese Zuschrift soll dem Zweck dienen, dieses 
Mißverständnis aufzuklären. 


Wir rekapitulieren zunächst ganz kurz diejenigen 


-Punkte der klassischen Stabilitätstheorie, in bezug 


auf welche Holstein und die Verfasser anscheinend 
übereinstimmen. 


Zuschriften an 
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Die grundlegende Differentialgleichung für kleine 
periodische Störungen einer ebenen, parallelen 
Strömung, 


DA) — 2020" +at!® ; \ 
—iaR{(U(y) —e) (8" — 020) _ ur —0) U) 


kann für große Werte von & R auf zwei wesentlich ver- 
schiedene Weisen gelöst werden. Es existiert erstens 
ein Fundamentalsystem ®;(y,aR), (i=1,..., 4), 
das asymptotisch für «R — 00 die Form 


PD, (y,aRk) 9,(Y) 
D,(yaR)9,(y) 

D, (war) "RAW a(y) 
D,(y,ak) 


hat. Hier sind gy,(y) und p,(y) Lösungen der reibungs- 
losen Gleichung 


U) — ep" — ap) —- U'g=0... (3) 


. (2), 


— lo 


nl 
m 


eV# RQ () o(y) ) 


und 
ET, 

Ay) = SViTW) eo) dy, mit Uly)—c=0 (4). 
Yo 


Y,(y) und p,(y) können so gewählt werden, daß p,(y) 
bei y=y, regulär ist, während 9,(y) dort eine Sin- 
gularität vom Typus (y— yo) log (y—y,) hat. Die 
Eigenschaften der Funktion o(y) sind unwichtig für 
die vorliegende Diskussion. 

Zweitens kann man die Differentialgleichung (1) 
der Transformation 


n= (y — Yo) (a R)Y? 


unterwerfen. Die transformierte Differentialgleichung 
besitzt dann ein Fundamentalsystem von Lösungen, 
die konvergente Potenzreihen in (a R)-1/3 sind. 

Die Störungen der parallelen Strömung müssen 
außer der Differentialgleichung (1) noch vier homo- 
genen Randbedingungen genügen. Eine Störung ist 
gedämpft, neutral oder angefacht, je nachdem ob 
c;—= Im c negativ, null oder positiv ist. 

Es besteht auch noch allgemeine Übereinstimmung 
darüber, daß im Innern der Strömung ein Gebiet 
existiert, das eine besondere Untersuchung erfordert, 
weil dort die reibungslose Gleichung möglicherweise 
keine ausreichende Beschreibung der wirklichen Strö- 
mung liefert. Die Meinungsverschiedenheiten be- 
treffen den mathematischen Ursprung und die Natur 
dieses „inneren Reibungsgebiets‘“. 

a) Holsteins Behandlung. Die Funktion 9,(y) 
hat eine Singularität beiy=y,. Da alle Lösungen der 
Differentialgleichung (1) bei y= y, regulär sind, kon- 
struiert Holstein, wie vor ihm Tollmien, eine Funktion 
9)(y,&R), die in der Umgebung von y=y, eine 
bessere: Annäherung an eine Lösung von (1) liefert. 
Diese Korrektur stellt nach Holstein den wesentlichen 
Einfluß der Viskosität in der inneren Reibungs- 
schicht dar. Nach [1], 8.,45 muß diese Korrektur 
nurinnerhalb einer Umgebung von y,, die für 
aR — © auf einen Punkt zusammenschrumpft 
berücksichtigt werden. 

b) Die Behandlung von Lin und Wasow. Im 
Gegensatz zu dem von Holstein untersuchten Phä- 
nomen, hat das innere Reibungsgebiet, das in [2], 
S. 131 und in [3], S. 870 diskutiert wird, mit dem 
Verhalten der Lösung ®,(y,aR) außerhalb einer 
Umgebung von y=y, Zu tun. Diese Umgebung 
kann beliebig klein sein, ist aber unabhängig von aR. 
Um dies klarer auszuführen, bezeichnen wir, wie 


in [3], die drei krummlinigen ‚Sektoren, in welche die. 


Kurve Re Q@(y)=0 die komplexe y-Ebene zerlegt, 
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mit S;f=1, 2,3). In [3], Satz 2 und 3 wurde 
bewiesen, daß die asymptotische Beziehung ®,(y,aR) 
9,(y) nicht in einem ganzen, den Punkt y=y, 
umgebenden Ring der y-Ebene gilt, sondern nur in 
zwei von den drei Sektoren S;. Im dritten Sektor 
nähert sich D,(y,‚«R) für «R— oo überhaupt 
keinem endlichen Grenzwert, vorausgesetzt, 
daß y außerhalb einer Umgebung von y, bleibt. 

Aus dieser Tatsache folgt leicht (s. [3], $. 870), daß 
es für gedämpfte Störungen ein inneres Reibungs- 
gebiet von endlicher Breite gibt (zum mindesten, 
wenn die Dämpfung nicht allzu klein ist), während für 
angefachte Störungen keine Reibungsschicht existiert. 

c) Holsteins Einwände. Die drei Kurven 
Re Q= 0 sind auch die Linien, über die hinaus die 
asymptotischen Darstellungen der Lösungen ®,(y,&R) 
und ®,(y,a&R) in (2) in unstetiger Weise fortgesetzt 
werden müssen (s. [3], Satz 1). Aber da der Beitrag 
dieser Lösungen zu einer wirklichen Störung, wie man 
leicht zeigt, für a R — oo nach null strebt, kann diese 
Eigenschaft von ®, und ®, nicht unmittelbar für das 
Auftreten des inneren Reibungsgebiets verantwortlich 
sein. Wenn wir richtig verstehen, ist dies der Kern 
von Holsteins Kritik auf $. 49 von [1]. Die Antwort 
ist, daß es in der Tat nicht das asymptotische 
Verhalten von ®, und ®,, sondern das von ®, 
ist, welches unser inneres Reibungsgebiet verursacht. 
Es ist allerdings richtig, daß die asymptotischen Eigen- 
schaften von ®, und ©, eine entscheidende Rolle in 
der Untersuchung des asymptotischen Verhaltens 
von ®, spielen (s. [3], S. 856—864). Der Zusammen- 
hang der Eigenschaften der Lösungen ®, und ®, mit 
denen von ®, und ®, hat auch Lin zu der Entdeckung 
des exzeptionellen dritten Sektors geführt. In Hin- 
blick auf Satz 3 in [3] ist es jetzt klar, daß die zwei 
inneren Reibungsschichten, die Lin. in [2], S. 133 er- 


- wähnt, in Wirklichkeit die Grenzen eines Reibungs- 


gebietes endlicher Breite sind. 

Holsteins Arbeit enthält noch einige andere Be- 
merkungen, mit denen die Verfasser nicht völlig über- 
einstimmen. Aber es ist vielleicht besser, die Dis- 
kussion auf die oben erwähnten Punkte zu beschrän- 
ken, in denen die Meinungsverschiedenheiten am klar- 
sten zum Ausdruck kommen. 
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Erwiderung auch, 

C. C. Lin kam in der oben unter [2] zitierten Arbeit 
zu dem Ergebnis, daß es bei laminaren Strömungen 
mit überlagerten Störschwingungen für den gedämpf- 
ten, neutralen und angefachten Fall 2 bzw. 1 bzw. 0 
innere Reibungsschichten gibt. Demgegenüber stellte 


‘ich in [1] fest, daß dieses Zahlentripel gleich 1,1, 1 


ist. Gemäß vorstehender Zuschrift von Lin und 
Wasow ist nunmehr Lins Ergebnis in 1, 1, 0 um- 


geändert. Hierdurch ist größere, wenn auch nicht 
bereinstimmung mit meinem Ergebnis _ 
Für eine Revisiön meiner kritischn B- 


vollständige 
erzielt. { 
merkungen in [1] sehe ich keinen Anlaß. 


Southsea, Hants, England. . H.Holstein, 


erkannt werden, anderenfalls wird das folgende Quartal noch geliefert. Verantwortlich für den Anzeigenteil: Curt F. W. Sc 
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